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influence the structural and control design of the aircraft (Houbolt, 1973; 
Etkin, 1972; Blsplinghoff and Ashley, 1957). Techniques for computing rolling 
and pitching moments and other aerodynamic forces, which are influenced by 
spatial turbulence, have been developed theoretically and, in general, utilize 
isotropic homogeneous turbulence (Diederich and Drischler, 1957; Eichenbaum, 
1972; Eggleston and Diederich, 1956; Houbolt, 1973; Lichtenstein, 1978; Kordes 
and Houbolt, 1953; Houbolt and Sen, 1972; Pastel, et al.. 1981; Akkarl and 
Frost, 1982; Diederich, 1957; Ringnes and Frost, 1985; Frost and Lin, 1983). 
It is normally accepted, however, that the turbulence in the atmospheric 
boundary layer close to the earth's surface, which is encountered by an 
aircraft during approach and takeoff, and turbulence associated with 
thunderstorms and clear-air roll waves is generally not isotropic. 
Additionally, turbulence shed by orographic features can also create 
relatively large-scale turbulence that is typically not Isotropic nor 

homogeneous. 

Spatial turbulence statistics have been computed from data measured with 
single tower to heights not exceeding much more than 100 m. Towers, however, 
provide spatial turbulence information only in the vertical (Davenport, 1961; 
Brook, 1975), which is uninteresting to aircraft design. Some studies have 
been carried out with tower arrays based on two or three towers located at 
various horizontal separation distances. The data normally reported from 
these studies is the coherence function for longitudinal velocities (Panofsky 
and Mizuno, 1975; Panofsky, et al., 1974; Krlstensen and Jensen, 1979; Pielke 
and Panofsky, 1970; Frost and Lin, 1983). These towers are normally less than 


1 


: - v '\ ;Sdm In height. Due to the fact that turbulence Information required 
for aircraft design Is at much higher altitudes (even in the terminal area 
data to heights of roughly 500 m are required), tower data are somewhat 
limited In their application. 

For this reason a NASA program has been underway to determine time 
histories and various statistical characteristics of three components of gust 
velocity measured simultaneously at the wing tips and at the nose - of a - 
specially instrumented B-57B airplane. The instrumentation system has been 
designed and installed on the airplane and several flights have been carried 
out (see Table 1.1). The flights include turbulence samples taken near storms 
in the Denver-Boulder, Colorado, area. Results from Flights 21, 22, and 26* 
are reported in considerable detail in Frost, et al. (1985a), Camp, et al. 
(1984), Frost (1983), Campbell, et al. (1983), and Chang, et al. (1986). 
Turbulence measurements with the aircraft during Flights 40, 44, 64, and 65 
have been compared with data obtained using remote radar sensing techniques 
(Frost, et al., 1985b; Huang, et al., 1985; Frost and Huang, 1983). Also, 
measurements of turbulent fluxes of momentum, heat, and moisture relative to 
orographic features were made during Flights 60, 61, 63, and 66. Analyses are 
presented in Chang and Frost (1985), Theon, et al. (1986), and Frost, et al. 
(1985c). 

The purpose of the present study is threefold: 

1. Perform statistical analyses of the acquired flight data with 
- emphasis on long data runs in continuous turbulence and glide 
slope runs for simulated takeoffs and landing approaches. 

Flight 31 flown at NASA Dryden was carried out specifically for 
this purpose. 


♦Flights 21, 22, and 26 were originally numbered 6, 7, and 10 

respectively, and are so referred to in the references cited. 




*Tor more information, please contact the Aeroservoelasticity Branch, 
NASA Research Center, 
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Study, develop, and/or modify statistical models, as necessary 
from the standpoint of providing an analytical expression for 
use in design analysis. 

3. Analyze effects of instrumentation characteristics and data 
processing effects on reduced gust velocity data. 

Flight 31 contained runs over sufficiently long distances at level 

flight to provide turbulence time histories long enough to assure a high 

statistical degree of freedom and to determine spectral characteristics for - 

wavelengths as great as 10,000 ft (3000 m). The meteorological correlations 

and statistical analyses of data from Flight 31 is the primary thrust of this 

report. Plans for Flight 31 were also to Include takeoff and touch-and-go 

runs to investigate non-stationary turbulence along the glide slope associated 

with the vertical variation of the horizontal wind In the atmospheric boundary 

layer (Panofsky and Dutton, 1984; Haugen, 1973). Non-stationary data calls 

for statistical ensemble analysis techniques (see Frost and Moulden, 1977; 

Wang and Frost, 1932). Ensemble statistics requires a collection, or 

ensemble, of sample records of the turbulence process. Wang and Frost (1982) 

have shown that a minimum of six flights down the glide slope under similar 

prevailing meteorological conditions are required to assure meaningful 

results. Unfortunately, ensemble analyses could not be carried out because 

during Flight 31 only one touch-and-go and one takeoff run were recorded. 

Thus, insufficient approaches or takeoffs were made under similar prevailing 

meteorological conditions to permit ensemble averaging. 

Another problem associated with the turbulence measurements carried out 
along non-level flight paths Is that the system of equations presently 
utilized to remove the airplane motions from the recorded data are based on a 
linearized model which assumes only small perturbations about wing level, 
horizontal flight. Analysis using the full non-linear system of equations 
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^^>s£n carried out for typical data and compared with the computations from 
the linearized system (Frost, et al., 1983). In general, the effects of a 
glide slope or climb-out angle less than 10* are negligible on the computed 
turbulence gust velocities. 


Data for all runs in Flight 31 including the touch and go are provided 
in Appendix A. Details of the flight path, the time histories, and selected 
statistical analyses including probability distributions, correlations, 
spectra, etc. are given In this appendix. The statistical analyses described 
In Section 2 were applied to the turbulence measurements for all runs. 

Although these analyses are strictly applicable to statistically stationary J 

data, little evidence of non-stationary effects Is observed in any of the data 
except for Run 10 as described later. This observation is true for the 


touch-and-go runs and takeoff run, also. 


The philosophy associated with Appendix A Is to provide the data after 
applying sufficient statistical analysis to allow the reader to distinguish 
data sets which are of interest to his specific application. The complete 
data can then be obtained on magnetic tape from NASA Langley Research Center 
(LaRC) for conducting the reader's own analysis. With this in mind, the 
complete data from Runs 1 through 16 of Flight 31 have been given In the 
appendix. Selected runs, however, are analyzed in more detail throughout 
Section 2 and compared with theoretical and empirical models currently 
available for correlating turbulence data. In general. It was not necessary 
to develop new theoretical models because the data fit existing models quite 
well, as also described in Section 2. There were, however, a few exceptions 
where modifications to two-point correlation and spectrum models were 
required. These are also described in Section 2. 


\ 
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his study, a distinction Is made between a two-point spatial 
correlation and the commonly used term "cross-correlation." The terminology 
cross-correlation in this report Is reserved for a correlation between 
different velocity components; for example, between the lateral and 
longitudinal components or the vertical and longitudinal components, etc. The 
terminology "two-point" spatial correlation refers to a correlation between 
velocity components measured at stations separated In space (e.g., at 
different probe locations on the aircraft). This can be either a two-point 
cross-correlation between dissimilar velocity components separated spatially 
or It can be a two-point correlation between like velocity components 
separated spatially.* A single-point spatial correlation termed an 
auto-correlation Is defined as a correlation between like velocity components 
measured at the same location (e.g., with the same wing tip probe) but 
separated in time. Note that a two-point spatial correlation can also be 
separated in time (i.e., time dependent or lagged in time). 

Procedures for estimating two-point spectra from finite digitized time 
histories are not straigntforward. Considerable Insight to these procedures 
which Is not readily accessible In the literature was gained during this 
study. This insight is described in detail In Section 3. 

The von Karman analytical correlation and spectrum model for atmospheric 
turbulence frequently referred to in the literature (Hinze, 1975; Houbolt, 


*The authors prefer the terminology two-point auto-correlation for a 
correlation between like velocity components measured at spatially separated 
positions. In deference to the reviewers, however, who found this terminology 
confusing and concluded that an auto-correlation must be a correlation of a 
signal with Itself, the correlation between the same velocity components 
measured at different positions in space Is called a two-point common 
component correlation or, where no confusion exists, simply a two-point 
correlation in this report. 
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' ~^>l964; Etkin, 1972; Panofsky and Dutton, 1984) Is used extensively In 


»is study. In Section 2, comparison of this model is made with the 


experimental data. 


In general, agreement with the von Karman auto- 


correlation and the one-point auto-spectrum is good. This Is surprising in 
view of the fact that von Karman’s model is generally assumed valid only for 
isotropic homogeneous turbulence. It should be noted in this regard that for 
the long duration, level runs of Flight 31 analyzed in the current study, the 
turbulence Is not expected to be Isotropic nor even homogeneous. Most runs 
were carried out over both flat desert and mountainous regions (peaks of 6000 
ft ASL) at low altitude (2500 to 10,000 ft ASL) or during touch-and-go flights 
at the airport. It was therefore surprising that a model based on the 
assumption of homogeneous isotropic turbulence correlated the data well. 

Although one-point spectra are addressed In some detail, emphasis in this 
report Is on the two-point spectra and correlation functions. Comparison of 
the twc-point correlations and spectra with a theoretical model (based on the 
von Karman model) originally proposed by Houbolt and Sen (1972) Is made In 
some detail. Correction to this model was required and made as described In 
Section 2. In general, the experimental data agree with the theoretical model 
after the corrections. Appropriate care must be exercised, however, in 
computing the two-point spectra from the digitized data. This Issue Is 
described in-depth in Section 3. 

During analysis of several of the flights, a number of instrumentation 
characteristics were uncovered which influenced the accuracy of the data. 
Although significant effort In the past has been devoted to evaluating effects 
of instrumentation characteristic and of data reduction procedures on the 
accuracy of the measured turbulence data, some additional factors were 


V 











^< v ' ' ‘-^^yne, 1980; Murrow and Rhyne, 1981; Meissner, 1976; Crooks, et 
kTi 967). The addition of the wing tip booms and the interest In measuring 
wind speeds as contrasted to fluctuations about the mean wind were, in 
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general, responsible for the new instrumentation problems. These problems 
along with recommended correction or removal procedures are discussed in 
Section 4. In particular, the INS Schuler position and velocity drift errors 
and the suspected flow vane sensor misalignment problems are addressed. 

Data from Flights 63. 66, 73, and 74 were used for analyzing the 
instrumentation errors. Data from these flights are used only for purposes of 
analyzing errors in this study. Also, a discussion of the Influence of 
departure from straight and level flight on the computed turbulence when the 
data are reduced using the linearized equations which are strictly valid only 


for level flight is given in Section 4 and in Appendix B. 
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2,0 STATISTICAL ANALYSIS OF DATA 


A statistical analysis of the B-57B aircraft data for 16 runs (Runs 1 
through 16) of Flight 31 on November 29, 1982, Is described In this section. 
The procedures for analyzing the turbulence data as well as interpreting the 
analyzed results are strongly influenced by the stationarlty of the data. 
Non-stationary or non-homogenous data represent all classes of data whose 
statistical properties change with time or with position. Figure 2.1 
illustrates three different examples of non-stationary data; these include 
data with a time-varying mean, data with a time-varying mean square, and data 
with a time-varying frequency structure (Bendat and Piersol (1971)). The vast 
majority of physical data actually fall into the former category. 

The theoretical ideas and processing techniques for the stationary data 
do not, for the most part, apply to data which are non-stationary. A totally 
adequate methodology does not exist yet for the analysis of all types Oi 
non-stationary data. In general, an ,'nsemble-averaging technique (Bendat and 
Piersol, 1971) provides a method to analyze the statistical properties of the 
non-stationary data (see application of this by Frost and Huang (1983)). By 
Inspecting the time histories of the aircraft-measured turbulence data (shown 
in Section 2.2), one can easily see considerable patchiness and non- 
stationarity in these data sets. However, only one sample record of 
turbulence over common terrain - and similar prevailing meteorological 
conditions Is available for analysis from each run of the B-57B aircraft data. 
Therefore, through necessity the statistical properties of the data presented 
in this report are calculated by assuming that the measured data are 


stationary. 



Figure 2. 1 




Gerieral information and statistical values for each of the 16 runs are 
£lven in Appendix A. The analysis for each run presented In the appendix 
consists of seven parts as follows: 


1. Flight altitude and horizontal wind velocity along the flight 
path. 

2. Time histories of gust velocities, gust velocity differences 
between wing tips, and the aircraft's normal acceleration... 

3. Average turbulence parameters, integral length scales, and 

correlation coefficients of gust velocities. 

4. Probability density functions for gust velocities and gust 

velocity differences. 

5. Normalized one- and two-point correlation functions of gust 

velocities. 

6. Normalized one- and two-point spectral density functions of 

gust velocities. 

7. List of all parameters measured and the range of their extreme 
and average values. 

A map illustrating all ground tracks for Flight 31 over terrain as 
recorded by the INS during the flight is provided in Figure 2.2. The cross 
section of the vertical profile of the terrain beneath the flight path is 
given fo” each run in Appendix A. 

The atmospheric stability is of importance in turbulence considerations. 
The temperature gradient in the atmospheric boundary layer is a measure of the 
stability of the atmosphere. Figure 2.3 shows the temperature recorded for 
all runs of Flight 31 superimposed on the temperature profile measured by the 
weather balloon. Each represents a 5-second averaged temperature. A 
scattering of the averaged temperature at different altitudes is seen in the 
figure. This scattering is believed to represent the spatial temperature 
variations along a flight path which usually covers more than 10 miles 
horizontally. Temperature profiles measured during the takeoff run and during 
the touch-and-go run. Run 2, were converted to potential temperature profiles 
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and are^plotted in Figure 2.4. The arrows indicate climb or descent. The 
strong negative temperature - gradient near 3700 ft, characteristic of an 
unstable boundary layer, as shown in the second half of Run 2 is believed to 


be associated with wake flow generated by a mountain peak up-wind of the 
flight path (see the terrain contours in Figure 2.2). 

The methods and results of the statistical analyses given in Appendix A 
are described and discussed in detail in the following subsections. The 


experimental data are also compared with theoretical models. 


2 . l Flight Altitude and Horizontal Wind Velocity » ~ong the Flight Path 

The first part of the analysis for each run in Appendix A Includes 
flight altitude (ASL), the corresponding terrain height (ASL), the flight 
direction, and five-second averaged horizontal wind vectors recorded along the 
flight path. The terrain height Is obtained from digitizing a large-scale 
contour map along each ground track of the flight as shown in Figure 2.2. 
Also, tabulated are the date, the time (PST) at which the run began, and the 
duration of the run in seconds. 

Run 1 Is the takeoff leg of the flight. Run 2 started with an approach 
and then made a go-around at approximately 100 ft above the ground. The 
approach and go-around flight path were at a glide .lope angle of 
approximately three degrees. The terrain features over which the majority of 
the B-57B Flight 31 experiment was flown are characterized by regions of low 
and high mountainous terrain (Runs 1 and 2 were over flat terrain). 

Data from these two runs, however, cannot be expected to be 
statistically stationary since they represent ascent and descent through 
variable wind conditions associated with the atmospheric boundary layer. 
There is no clear evidence in the results of the statistical analysis of these 
data sets, however, that suggest non-stationary or even non-isotropic effects. 
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Note also'ln Runs 3. 4, 8, 9, and 11 there are what appear to be relatively 
large excursions from level flight. This Is primarily an Illusion due to the 
exaggeration of the vertical scale In the plot. The apparent climb and 
descent paths are In all cases less than seven degrees, elevational angle. 

Inspection of the flight paths given In Figure 2.2 shows that they may 
be categorized as occurring over a flat region, a low mountainous regions, and 
a high mountainous regions. The low mountainous region Is further subdivided 
into two subregions. The first, subregion is that for which the underlying 
terrain includes gradual and monotonic increase or decrease in elevation. The 
second subregion is one which includes terrain having more than one peak or 
valley. Table 2.1 categorizes the type of terrain associated with each run. 


TABLE 2.1. Terrain Category for Flight 31. 


Terrain Cateaory 

Run Number 

Flat Region 

Low Mountain Region 

1, 2 

Single Peak 

10, 15, 16 

Multi-Peak 

5, 6, 9, 11, 
12, 13, 14 

High Mountain Region 

3, 4, 7, 8 


The flight paths plotted from the INS data for Runs 4, 8, 11, 13, and 14 
shown in Figures A.16, A.36, A.51, A.61 and A.66 suggest that the aircraft 
flew through mountain peaks. The cause of this obviously, impossible result is 
associated with an INS drift problem which Is discussed in Section 4. Errors 
in the recorded longitude and latitude measurements result In an Incorrect 
aircraft position relative to the fixed terrain features (see Figure 2.2). 
However, only when the influence of terrain on the turbulence Is to be 
assessed does the error Influence the data analysis. 
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Table 2.2 shows the time duration of each run and the effective mean 
wind direction relative to the airplane. In. most runs the effective mean wind 
direction Is nearly perpendicular to the flight path. The measured mean wind 
speed in Flight 31 ranges from 5 m/s to over 20 m/s. Run 3 is the longest of 
the 17 runs making up Flight 31. The landing leg of the flight was recorded 
as Run 17. This run was only 47 seconds, which is not statistically 
meaningful and, therefore, is not analyzed in this report. 

TABLE 2.2. Time Duration and Mean Wind Direction. 


Run 

lumber 

Length of 
Record 
(sec) _ 

1 

135 

2 

213 

3 

694 

4 

283 

5 

144 

6 

63 

7 

203 

8 

226 

9 

334 

10 

172 

11 

333 

12 

138 

13 

270 

14 

209 

15 

233 

16 

100 

17 

47 


Wind Direction Observed 
bv the Airplane 

Head wind . , _ 

Cross wind (left to right) 
Cross wind (left to right) 
Cross wind (right to left) 

Cross wind (left to right) 
Cross wind (right to left) 
Cross wind (left to right) 
Cross wind (right to left) 

Cross wind (right to left) 
Cross wind (left to right) 
Cross wind (left to right) 
Tail wind 

Cross wind (left to right) 
Cross wind (right to left) 
Head wind 
Tail wind 
Head wind 


1 

\ 


i 

l 



■ ■ ■ u, ! - [ ’ 'US 


2.2 Time Histories of Gust Velocities. Gust Velocity Differences Between Wing 
Tips, and the Aircraft's Normal Acceleration 

The second part of the analysis for each run recorded in Appendix A 
shows the gust velocity time histories for the three probes located at the 
aircraft's nose and wing tips. The left, center, and right probes are 
designated with subscripts L, C, and R, respectively. The time histories of 
the spatial velocity differences between the right and left probes are plotted 
for the longitudinal (u), lateral (v), and vertical (w) velocity components. 
The definition of longitudinal, lateral, and vertical are along and 
perpendicular to the mean flight path, respectively. Also plotted with the 
vertical velocity time histories is the time history of the aircraft's normal 
acceleration along Its flight path. The sampling rate Is 40 samples per 
second. 

Figure 2.5 shows a plot of the time histories for Run 3 of Flight 31. 
Run 3 was the longest record lasting 694 seconds. The velocity fluctuations 
are typical of the measured data. One observes from the data that there are 
no significant variations between velocities measured at the three probes. 
Therefore, one can surmise that length scales associated with these turbulence 
data are typically larger than the wing span (19.5 m). In Run 3, the aircraft 
encountered significantly more intense turbulence from 560 seconds to 640 
seconds. The turbulence was encountered at approximately 11:00 a.m. (PST) 
just after the aircraft had climbed from 5600 ft to 6600 ft at an approximate 
7* climb angle. During the climb, which started at 530 seconds, visual 
Inspection shows no dlscernable change In the turbulence during the climb and 
for roughly a half mile after leveling off at 547 seconds. During the period 
from 560 seconds to 640 seconds, however, the aircraft encountered much 
stronger turbulence as it flew over a high mountainous area with well 
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the aircraft's normal acceleration for Run 
November 29, 1982. 
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Figured. 5. (continued). 
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Figure 2.5. (continued). 






£c) Vertical component and aircraft's normal acceleration. 



pronounced peaks reaching up to 6500 ft. The wind was blowing over the 
mountains perpendicular to the aircraft flight path (see Figure 2.2). It is 
believed that the Intense turbulence encountered by the aircraft during this 
time period was associated with the disturbed air flow from the nearby 
mountain peaks. The strong vertical turbulence Induced an increased 
fluctuation to the aircraft's normal accelerations. 

2.3 Average Turbulence Parameters. Integral Length Scales, and Correlation 
Coefficient of Gust Velocities 

The third part of the analysis for each run in Appendix A is a tabulated 
listing of the average values of several important turbulence parameters for 
the left, center, and right probes. The statistical parameters include mean 
airspeed, standard deviation of gust velocity, standard deviation of gust 
velocity difference. Integral length scale, and the correlation coefficient of 
the gust velocity. The mean airspeed and the standard deviations of the gust 
velocities and their differences are calculated on the basis of the total time 
history. However, In analyzing the data to obtain the correlation coefficient 
and the integral length scale, the total time history is segmented such that 
the total record Is a multiple of segments of 1024 datum points. In computing 
the correlation coefficients and the length scales, (modify: .... scales, 

the spectrum was first computed by a technique which applied the Fourier 
transform directly to the original digitalized data. The correlation Is then 
....). The approach first computes the spectrum directly from the 
turbulence time history. The correlation Is then computed from the Inverse 
Fourier transform of the spectrum. Finally, the length scale Is computed by 
Integrating the normalized correlation function as described later. 

Table 2.3 lists the mean airspeed for all 16 runs In Flight 31. The 

average mean airspeed for all runs Is 102 m/s. The mean airspeeds at the 

25 


TABLE 2.3. Mean Airspeed (m/s) for Flight 31, November 29, 1982. 


Run 

Number 

V L 

v c 

% 

1 

81.13 

78.92 

81.21 

2 

87.82 

85.79 

87.51 

3 

104.21 

102.52 

104.60 

4 

104.78 

102.62 

104.32 

5 * 

105.79 

103.53 

105.33 

6 

104.31 

102.19 

104.01 

7 

101.47 

99.23 

100.93 

8 

103.22 

101.05 

102.86 


Run 

Number 

\ 

v c 

Vr 

9 

103.15 

100.84 

102.84 

10 

117.27 

115.20 

116.70 

11 

107.01 

104.47 

106.49 

12 

101.03 

98.55 

100.56 

13 

103.30 

101.40 

103.30 

14 

103.38 

101.07 

102.99 

15 

107.74 

105.40 

107.23 

16 

109.41 

107.07 

108.82 


individual right, center, and left probes are 102.8 m/s, 100.6 m/s, and 102.5 
m/s. The mean airspeeds measured at the right and left probes are larger than 
that at the nose by about 2 m/s. This difference can possibly be due to flow 
deceleration In front of the aircraft nose and/or flow acceleration over the 
wlngtlps. Approximate potential flow analysis for a Ranklne body (Karmacheti, 
1966) suggests the former mechanisms. Similar velocity differences were also 
found by Frost, et al. (1985a) In the analysis of Flight 21 for the same 
experimental aircraft. 

Table 2.4 lists the standard deviation of gust velocities for all 16 
runs of Flight 31. The standard deviation of the gust velocities varies from 
1.68 to 7.46 m/s for the longitudinal component, and from 1.42 to 5.57 m/s for 
the lateral component. The vertical gust component standard deviation ranges 
from 1.08 to 3.45 m/s. Table 2.4 also lists the standard deviations of the 
gust velocity differences between the right and left probes. The standard 
deviation of the gust velocity differences has a high of 1.75 m/s for the 
longitudinal component, 1.68 m/s for the lateral component, and 2.00 m/s for 
the vertical component. The standard deviation of the gust velocity. Itself, 
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TABLE 2.4. Standard Deviation (m/s) of Gust Velocity and Gu*t Velocltv 
Difference for Flight 31. * 3 


Run 

Number 


1 

2.12 

2 

3.23 

3 

3.17 

4 

3.73 

5 

2.49 

6 

3.64 

7 

3.03 

8 

3.93 

9 

4.10 

10 

2.04 

11 

3.74 

12 

1.68 

13 

2.49 

14 

2.51 

15 

7.46 

16 

5.68 


Run 


Number 

°uL 

1 

1.74 

2 

3.20 

3 

3.19 

4 

3.75 

5 

2.56 

6 

3.54 

7 

3.07 

8 

3.89 

9 

4.18 

10 

2.02 

11 

3.77 

12 

1.70 

13 

2.59 

14 

2.52 

15 

7.32 

16 

5.74 


°vR 


2.11 

2.31 

2.03 

1.16 

5.25 

2.36 

4.07 

2.80 

4.06 

2.76 

3.67 

3.41 

3.00 

2.23 

5.17 

2.52 

5.12 

2.^0 

4.57 

2.40 

2.10 

2.25 

1.43 

1.66 

5.57 

2.43 

3.54 

2.37 

2.84 

3.45 

3.21 

3.21 


°vL 

°wL 

2.05 

2.58 

2.09 

1.17 

5.31 

2.31 

4.04 

2.79 

4.10 

2.85 

3.42 

3.12 

3.06 

2.22 

5.20 

2.42 

5.12 

2.34 

4.61 

2.33 

2.18 

2.13 

1.42 

1.59 

5.56 

2.41 

3.42 

2.28 

2.87 

3.35 

3.29 

3.14 


a uC 

®vC 

°wC 

1.69 

1.99 

2.33 

3.20 

2.01 

1.08 

3.15 

5.20 

2.18 

3.73 

4.09 

2.61 

2.47 

4.04 

2.66 

3.54 

3.65 

3.00 

3.03 

3.03 

2.15 

3.89 

5.18 

2.36 

4.10 

5.10 

2.21 

1.99 

4.58 

2.34 

3.76 

2.15 

1.99 

1.68 

1.47 

1.54 

2.48 

5.57 

2.29 

2.47 

3.50 

2.12 

7.31 

2.89 

3.29 

5.59 

3.44 

3.02 

°AuRL 

a AvRL 

°AwRL 

1.20 

1.10 

0.77 

0.94 

0.77 

0.87 

1.29 

1.23 

1.37 

1.59 

1.39 

1.62 

1.41 

1.38 

1.42 

1.74 

1.68 

1.92 

0.85 

0.80 

0.89 

1.22 

1.08 

1.31 

1.35 

1.16 

1.45 

0.41 

0.31 

0.38 

1.12 

1.01 

1.26 

0.90 

0.74 

0.91 

1.53 

1.39 

1.59 

1.29 

1.12 

1.37 

1.45 

1.24 

1.49 

1.75 

1.61 

2.00 


t 
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1$ always larger than the standard deviation of the gust velocity difference 
between probes. 

Table 2.5 lists the integral scales L for 16 runs of Flight 31. 
The turbulence Integral length scale is usually estimated by Integrating a 
normalized one-point auto-correlation function from zero to infinity with 
respect to temporal or spatial lag. The normalized correlation function (also 
called the correlation coefficient), B x (t), Is given by 

B x (t) - x(t)x(t+x)/a x o x 

where x Is any of the velocity components. Due to noise In the measured data, 
the auto-correlation coefficient, however, nearly always oscillates about zero 
due to either real physical effects but most probably due to aliasing and 
other digitizing effects. Therefore, In this report the Integral length scale 
Is obtained by Integrating the normalized auto-correlation function to the 
point where It first crosses zero (2. = S * VT or t = T): 

S T 

L = J B x (*)dl = V J B x (x)dx (2-1) 

0 0 

Detailed study of different definitions of the integral length scales Is 
given in Frost and Lin (1983) (also see Houbolt, et al. (1964)). Frost and 
Lin (1983) suggest that using the L as defined In Equation 2.1 In theoretical 
models gives best agreement with experimental results. This length scale Is 
therefore used throughout this report. 

In addition to the Integral scales calculated from the normalized values 
of one-point auto-correlation functions, UrUr, VrVr, and WrWr, the Integral 
length scales are also estimated with Equation 2.1 using the normalized values 
of the two-point correlation functions, UrOl, VrVl, and WRWj;. These two 
integral scales have the same order of magnitude for each corresponding 
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TABLE 2.5. Turbulence Length Scales for Flight 31. 


Integral Length Scale (m) 


Run 

Number 

l uR 

l vR 

1-wR 

l uRL 

l vRL 

l wRL 

i 

297.7 

149.7 

255.1 

248.4 

35.3 

254.5 

2 

325.9 

250.1 

79.1 

322.4 

251.8 

89.3 

3 

234.0 

425.6 

116.9 

258.4 

422.8 

115.3 

4 

419.8 

350.8 

66.9 

408.0 

344.7 

61.9 

5 

333.9 

168.6 

189.7 

317.5 

173.5 

204.0 

6 

364.7 

92.0 

51.7 

344.5 

104.2 

47.5 

7 

562.8 

249.6 

287.6 

532.2 

242.9 

283.5 

8 

306.7 

364.4 

232.9 

302.5 

380.7 

249.6 

9 

327.8 

338.0 

93.9 

341.5 

338.0 

83.9 

10 

641.8 

729.7 

832.2 

638.3 

742.9 

863.8 

11 

370.0 

246.1 

203.3 

375.6 

241.7 

193.1 

12 

127.7 

252.6 

202.4 

137.5 

250.1 

190.8 

13 

156.0 

428.8 

83.7 

148.6 

424.4 

82.6 

14 

174.9 

204.4 

66.8 

161.3 

205.4 

64.5 

15 

540.0 

225.8 

526.1 

526.5 

225.3 

494.0 

16 

348.1 

362.2 

95.0 

347.5 

336.5 

115.3 


turbulence velocity component (see Table 2.5). The Individual velocity 
component characteristics do not vary appreciably across the wing span which 
Is In agreement with the fact that the calculated length scales are much 
larger than the 19.5 m wing span of the aircraft. This Implies that the 

energy-containing turbulence fluctuations essentially engulf the total 
airfoil. 

Finally, Table 2.6 shows the two-point correlation coefficients of the 

gust velocities computed for 16 runs of Flight 31. The symbols PrU L , vrVl, 

and wrw l represent the two-point common velocity component correlation 

functions for longitudinal, lateral, and vertical components, respectively, 

whereas UrVr, and wrIIr represent the one-point cross-correlation 

functions, and urVl, VrWl, and represent the two-point cross -correlation 

functions. Although several other correlations of the gust velocities could 
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TABLE 2.6. Two-Point Correlation Coefficient of Gust Velocity for Flight 31. 


Run 

U R U L 

VrVl 

WR»L 

UrVr 

VrWr 

*R°R 

UrVl 

VrWl 

w R u L 

Number 

0 VR° V L 

°WR a WL 

a UR°Vl_ 

° v R 0w R 

OwrOur 

d UR a v R 

^ v R a WL 

°WR a ut 

a WR°UL 

1 

0.75 

0.34 

0.80 

0.08 

0.14 

0.52 

-0.40 

0.11 

0.45 

2 

0.80 

0.81 

0.82 

0.00 

-0.05 

0.11 

-0.02 

-0.03 

0.10 

3 

0.80 

0.91 

0.75 

0.09 

-0.19 

0.06 

0.04 

-0.19 

0.05 

4 

0.88 

0.91 

0.80 

-0.19 

0.20 

0.09 

-0.19 

0.20 

0.06 

5 

0.87 

0.90 

0.90 

-0.09 

-0.10 

-0.17 

-0.08 

-0.09 

-0.20 

6 

0.82 

0.90 

0.81 

-0.18 

0.60 

-0.10 

-0.18 

0.61 

0.00 

7 

0.92 

0.90 

0.90 

0.02 

-0.21 

0.39 

0.01 

-0.21 

0.32 

8 

0.89 

0.81 

0.80 

-0.20 

0.15 

0.00 

-0.19 

0.10 

0.03 

9 

0.80 

0.90 

0.80 

0.30 

0.20 

0.19 

0.30 

0.20 

0.18 

10 

0.98 

0.99 

0.98 

0.08 

0.00 

-0.47 

0.09 

-0.01 

-0.45 

11 

0.83 

0.81 

0.78 

-0.21 

-0.09 

0.48 

-0.28 

-0.10 

0.40 

12 

0.66 

0.81 

0.78 

0.00 

0.30 

-0.22 

0.01 

0.31 

-0.20 

13 

0.80 

0.91 

0.79 

-0.18 

-0.32 

0.25 

-0.19 

-0.32 

0.22 

14 

C.79 

0.90 

0.77 

0.18 

0.19 

0.10 

0.13 

0.27 

0.07 

15 

0.90 

0.88 

0.88 

0.32 

0.06 

0.02 

0.30 

0.01 

0.00 

16 

0.85 

0.86 

C.85 

0.49 

0.05 

-0.10 

0.49 

0.02 

-0.10 


be estimated, the combinations shown in Table 2.6 are sufficient to detect any 
trends or physical effects associated wlthr the normalized spatial correlation 
computed from these data. Note that the appreciable difference In value 
between the one-point auto-correlation evaluated at zero lag, (shown In Table 
2.4) and the two-point correlations evaluated at zero lag (shown in Table 2.6) 
is that correlation coefficients (i.e., normalized values) are tabulated in 
Table 2.6 whereas non-normal 1 zed correlations are given In Table 2.4. It is 
clear from Inspection of Table 2.6 that the correlation between like 
components of turbulence has a roughly uniform decrease In value of 20 percent 
over the wing span of the aircraft. 

All of the two-point correlation coefficients between common velocity 
components are larger than 0.75 except the value of VRVjVo VR a VL for Run 1 and 
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■DrDl/^ur^ul for Run 12m The fonTier "^y be associated with Run 1 being a take- 
off flight path (see Figure A.l). No explanation is evident for Run 12. The 
one-point cross -correlation coefficients are the Reynolds stresses (Frost and 
Moulden, 1977; Hinze, 1975) and are thus a measure of momentum transfer. For 
Isotropic turbulence the cross-correlation terms theoretically are zero. The 
very low values shown In Table 2.6 suggest that the atmospheric turbulence Is 
Indeed nearly Isotropic. It Is believed that the large value of vrWr /a v 
o w and 

R R 

VRW[/®VR a WL f or Run 6 is caused by the very short averaging time of 63 seconds 
associated with the run. It therefore does not represent a meaningful 
statistical average. Similar arguments can be made for other unjustifiably 
large values of the cross-correlation coefficient. 

2.4 Probability Density Function for Gust Velocities and Gust Velocity 

Differences 

The fourth part of the analysis for each run in Appendix A contains the 
probability density function of the turbulent wind velocities. Data measured 
by the B-57B aircraft for all three different probe positions and for all 
three velocity components (longitudinal, lateral, and vertical) are plotted. 
The degree of or lack of normality of the turbulent wind velocities is 
Illustrated by comparing the experimental probability density distributions 
with the theoretical normal distribution and the theoretical non-Gausslan 
probability density model (modified Besse 1 function distribution see Reeves, 
et al. (1974)). 

The probability density function for the turbulence wind velocities is 
defined as: 

n / y \ _ 11m Probrx < xft) < x + Axl 
p ' ' Ax-0 Ax 


( 2 . 2 ) 


where x(t) may be u, v, w, Au, Av, or Aw and where Prob[x < x(t) < x + Ax] Is 
the probability that the turbulence wind velocity at time t lies within a 
specified speed Interval. The Gaussian probability density function Is given 
by: 

_ (x-x) 2 

P(x) - -±z e 2 ° 2 (2.3) 

Cv£x 

where x Is the mean value of x(t) and o Is the standard deviation of x(t). In 
calculating the probability density distributions In Appendix A, the gust 
velocity and the gust velocity differences are normalized with their standard 
deviations. 

The non-Gausslan probability density distribution Is given by Reeves, et 
al. (1974) as: 



1 + 2t 2 r 2 


exp 




+ 


(1 + 2i 2 r 2 ) 




(2.4) 


where r Is an adjustable parameter which Is a measure of the degree to which 
the distribution Is non-Gauss ian and l is the dummy variable of Integration. 
If r » 0, the function Is exactly the Gaussian function; however, as r 
Increases, the distribution departs from the Gaussian probability density 
function and approaches a modified Bessel function distribution as r-*-. 

Figure 2.6 shows typical probability density distributions for Run 3 of 
Flight 31. The upper half of the figure shows the probability density 
distributions of the three individual gust velocity components and the bottom 
half of the figure shows the probability density distributions of the gust 
velocity differences between the right and left probes. These probability 
density calculations for the measured turbulence do not fit the normalized 
Gaussian distribution very well. Fitting the Individual probability 


32 



Probability 


'•ii . 


Gaussian 
Non-Gauss ian 


R = Right 


R = Right 


R = Right 
C = Cente 
L = Left 


S 0.6 
3 

es 

2 0.4 


Gust Velocity Distributions 


r = 1 .5 

* i < i i i rrr 

+ Measured 
Values 
“ U R-L_ 



r = 1.5 

1 1 l, i i f i 

+ ileasured 
Values 

av r-l _ 



+ ileasured 
Values 
R-L_ 


v 


Gust Velocity Difference Distributions 

Figure 2.6. Comparison of probability density functions for gust 

velocities and gust velocity differences (normalized with 
the standard deviation; with theoretical models, Fliqht 
31, Run 3 (r = degree of non-Gausslan). a 







i 


distributions to Equation 2.4 by adjusting r, a non-Gausslan form which 
provides a closer fit to the experimental data Is found. Inspection of Figure 
2.6 and similar figures In Appendix A shows the non-Gausslan distribution 
gives a very good fit of the gust velocity difference probability 
distribution. The Individual gust velocity probability distribution however. 
In many cases, appears to be blmodal. This Is associated with trends In the 
mean velocity that have not been removed from the data. 

Runs 1 and 9 are clear examples of the effects of trends In the mean 
wind on the probability distribution. For Run 1 during clinfc-out, the 
longitudinal mean wind which Is essentially a headwind (see Figure A.l, 
Appendix A) will Increase from zero at the surface to the value aloft. This 
vertical variation In the mean wind will typically vary logarithmically 
(Panofsky and Dutton, 1984). By simply assuming the mean wind speed Is 
uniform and removing a constant value from the data (as was done In this 
study) causes the velocity fluctuations about the mean at low levels to be 
mainly negative and at higher levels mainly positive fluctuations (see Figure 
A.2, Appendix A). Thus, there Is a blmodal distribution In the probability 
density function of the velocity fluctuations. This blmodal effect can 
probably be eliminated by removing a logarithmic velocity profile* trend. 
However, this was not done. 

Now consider Run 9. The quasi-steady horizontal wind speed along the 
flight path Is shown In Figure A.41. For the Initial part of this flight, the 
winds were partially headwinds with a dominate northward direction. During 
the latter part, the winds became partially tailwinds with a westerly 


♦Note vertical variation of horizontal 
logarithmic (Panofsky and Dutton, 1984). 


wind with height Is typically 
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direction. Again, removing a uniform average wind speed from the data results 
In the longitudinal velocity fluctuations being mainly negative during the 
Initial part of the flight and positive during the latter part (see Figure 
A.42). Again, this results In a strong blmodal distribution In both the 

longitudinal and lateral wind speed gust distributions as shown In Figure 
A. 43. 

Returning to the discussion of the analytical models which best fit the 
data. It Is clear that Equation 2.4 fits the experimental data considerably 
better than the Gaussian distribution. The value of r which gives the best 
fit of the data changes for different velocity components and from run to run. 
Table 2.7 lists values of r determined from "eye-ball" fits of the data for 
the three components of the gust velocity and gust velocity differences. The 
variation In r might be expected to be a result of the underlying surface 
roughness. Inspecting of the terrain features beneath each flight path, 

however, suggests ro apparent relationship between surface roughness and the 
value of r, nor Is there an obvious correlation between r and altitude. 
Further work Is required to associate the degree of non-Gauss lanness of the 
atmospheric turbulence with physical causes. 

2.5 Normalized One- and Two-Point Corre lation Functions of Gust VelnrH-W 
The fifth part of the analysis for a given run In Appendix A Is the 
normalized one- and two-point correlation functions of the turbulent wind 
velocities at the right and left wing tips. The correlation function between 
the same velocity components at two different positions separated by a vector 
distance c Is defined as (Panchev, 1571; Hlnze, 1975); 

T 

Rx(c»t) * Y J x(£,t) x'(^+c, t+r)dt 
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T 




where t Is the lag time, x and x' designate any one of the velocity components 

u, v, and w, and £ Is the position vector at which the velocity x Is measured. 

TABLE 2.7. Values of r (Equation 2.4) Which Represent a Measure of the 

Degree of Departure from a Gaussian Probability Distribution of 
the Gust Velocities. 


UK 

/ 


Run 

Number 


u R» u C> u i V R> V C» V L WR.w c ,w|_ Au r _ l Av r _ l Awj^l 


1 

1.5 

2 

2.5 

3 

0.5 

4 

2.5 

5 

3.5 

6 

3.5 

7 

2.5 

8 

4.5 

9 

2.5 

10 

1.5 

11 

2.5 

12 

1.5 

13 

1.5 

14 

3.5 

15 

4.5 

16 

1.5 


1.5 

1.0 

1.5 

1.0 

3.5 

4.5 

2.5 

1.5 

1.0 

3.5 

1.5 

1.0 

2.5 

3.5 

4.5 

1.5 


1.5 

1.0 

1.5 

2.5 

3.5 

4.5 

2.5 

1.5 

1.0 

2.5 

1.5 

1.0 

1.0 

3.5 

1.5 
1.5 


3.5 

1.5 

1.5 

2.5 

4.5 
4.5 

4.5 

2.5 

2.5 

2.5 

2.5 

1.5 

1.5 

4.5 

7.5 

1.5 


3.5 

1.5 

1.5 

2.5 

7.5 

6.5 

6.5 

2.5 

2.5 

4.5 

2.5 

1.5 

1.5 

4.5 

7.5 

1.5 


0.5 

1.5 

2.5 

3.5 

6.5 

7.5 

6.5 

3.5 

2.5 
2.5 

2.5 

1.5 

2.5 

4.5 

7.5 

1.5 
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The normalized correlation function: 


B x (£»t) 


R x (c,t) 

°x®x' 


( 2 . 6 ) 


Is called the correlation coefficient, c x and c x ‘ are the standard deviations 
of x(£,t) and x'(£+i-,t), respectively. If £ is not equal to zero, and x and 
x* are the same velocity component, the correlation function Is called the 
two-point common component correlation function In this report. The absolute 
value of B x (c,t) Is always less than one. At £ « 0, R x (£,t) reduces to the 
one-point auto-correlation function given by: 


(2.7) 


r x( t ) - y / x(t) x(t + x)dt 
0 


The correlation function may be evaluated using a time history sunroation 
technique or by using the direct Fourier transformation method for computing 
spectra. Steely and Frost (1981) and Frost and Lin (1983) have compared the 
direct method with the summation technique and found they give Identical 

results. The direct method Is therefore used throughout this report unless 
otherwise stated. 

Theoretical models of the one-point auto-correlation and of the 
two-point correlation are the von Kanaan model (Hlnze, 1975) and the Houbolt 
and Sen (1972) extension of the von Karmen model, respectively. The von 
Karman theoretical model for the normalized one-point auto-correlation 
functions for longitudinal and transverse velocity components Is expressed as: 

- '[£] 1/3 ' ; >/3(sy 


!l(! > Mff] - (&7] (3.8, 

where c = 2 /r(l/3), a * 1.339, K Is a modified Bessel function of the 

second kind, r Is the gamma function, L is the Integral length scale, and c Is 
the spatial lag distance. The subscripts L and T refer to longitudinal and 
transverse, respectively. The longitudinal and transverse velocity 
correlation are defined as shown In Figure 2.7. 

The von Karaan correlation Is In principle only valid for Isotropic 
turbulence. The more general non-lsotroplc velocity correlation Is a second 
order tensor given by (Hinze, 1975) as: 

R lj(5.0 = ITtfU Uj u + O ei(0 ejU + O (2#g) 

3? 


I 







where ei (0 represents the direction cosines of the velocity vector at the 
position 5 with respect to the 1 th axis and ej (5 + c) Is similarly defined at 
a distance c from the position £. The symbol u-f(£) is the Instantaneous 
component of the velocity fluctuation with respect to the mean at the position 
5 and ujU + O is similarity defined. The general correlation, RijU.O, is 
thus described Jn_terms of nine components. When the turbulence is isotropic 
and homogeneous it can be shown that the correlation can be expressed solely 
in terms of the longitudinal and transverse correlations shown in Figure 2 . 7 . 

In the present investigation the velocity components are expressed 
relative to the axis of the aircraft (the assumption of small angles is evoked 
(see Appendix B)). For the longitudinal and transverse correlations the 
velocity components must be resolved parallel and perpendicular to the line 
between the two measuring points as illustrated in Figure 2 . 8 . Therefore, to 
transform the longitudinal and transverse correlations to the aircraft frame 
of reference, the cosines in Equation 2.9 must be taken into account. 

Frost, et al. (1985a) have shown following Hlnze (1975) that for 
isotropic turbulence (see Figure 2.8): 


R uU) - _ «T 2 Bt(c) + — £ ol 2 B l (c) 

C z + s z + s z 

( 2 . 10 ) 

R v(c) * ? ~ °L 2 b l( 0 + — r ot 2 b t(c) 

5 z + S Z 5 Z + S 2 

( 2 . 11 ) 

The vertical velocity correlation Is, of course: 
R w(c) - oj 2 Bj(c) 

( 2 . 12 ) 


This model is referred to as the Houbolt and Sen model since Houbolt and Sen 
(1972) utilized it with Equations 2.7 and 2.8 early on to develop a two-point 
spectrum for use in design analyses. (It should be noted that in actual fact, 
Houbolt and Sen did not account for the direction cosines and hence their 
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longitudinal spectrum Is Incorrect.) Fcr Isotropic turbulence, ay * * a 

which Is not the case for the experimental data as Is apparent from Table 2.4. 

Figure 2.9 shows normalized one- and two-point correlation functions. 
The correlations are from Run 3, Flight 31. All correlation coefficients are 
plotted versus the spatial lag distance, c * Vt, In the direction of flight. 
The normalized one-point auto-correlation functions are plotted In the upper 
part of the figure, and the normalized two-point correlation functions In the 
lower part. The two-point correlations have both negative and positive time 
lags. Only the positive lag Is given In the figure. Negative lags behave 
similarly but are not synmetrlc. The Influence of negative lag appear In the 
phase angle of the two-point spectrum which Is discussed In a later section of 
this report. The area obtained by Integrating the one-point auto-correlation 
coefficient from zero spatial lag to the point where the correlation 
coefficient first crosses zero Is defined as the Integral length scale (see 
Table 2.5). Comparisons of the experimental data with the von Karman 
theoretical one-point auto-correlation coefficient and with the Houbolt and 
Sen (1972) theoretical two-point correlation function are shown In the figure. 
The Integral length scale, L, used In the theoretical models was that 
determined as described above. Using length scales determined from other 

definitions (see Frost and Lin, 1983) gave no better and. In most cases, 
poorer agreement with the experimental data. 

In general, auto-correlations are expected to decay faster for the 

i 

vertical and lateral components than for the longitudinal component. Results 
of Run 3 shown In Figure 2.9 appear to be an exception to this rule since 
similar plots of correlations for other runs given In Appendix A behave as 
expected. However, employing the length scales, computed as described In the 
theoretical model, results In the one-point auto-correlation coefficient 


Aircraft data 
von Karma n 



gust velocities with theoretical models. Flight 31, Run 





fitting the measured data quite well (see for example Runs 6, 14). The 
experimental correlation coefficient does depart, however, from the theory. 
It Is higher than the value predicted by the von Karman model at the larger 
spatial lags. 

It Is interesting that poorer agreement with von Karman's theoretical 
models occurs for the high altitude flight Run 10 than for the others. This 
Is surprising since it is generally assumed that turbulence at higher 
altitudes Is Isotropic. The wind at the higher level, however, may have been 
stratified with embedded gravity waves. This Is suggested by the high degree 
of correlation shown in Figure A. 49 and inspection of the time history In 
Figure A.47 which suggests Run 10 encountered a wave pattern. 

Consideration of Figures 2.7 and 2.8 shows that R u (c) defined by 
Equation 2.10 converges to ot 2 Bx( 0 at K » 0. Inspection of the correlation 
coefficients plotted In Appendix A shows this to be approximately true in most 
cases. In turn, as £ becomes large R u (c) approaches the longitudinal 
correlation aL 2 B|_(c) - This Is also approximately true based on Inspection of 
the experimental results. The above observation suggests that the turbulence 
Is reasonably isotropic for all runs except Run 10. Run 10 at high altitude 
as noted appears to be associated with wave motion. This Is even more 
apparent In the cross-correlation coefficients described next. 

The cross-correlation function of two sets of random data describes a 
general dependence between the variations of the sets. The two-point 
cross-correlation function Is given by: 

T 

Rxy(C') - T J x(Ct) y ({+?, twjdt (2.13) 

0 

where x and y are time histories of any two of the turbulence velocity compo- 
nents u, v, and w, t represents the lag time, and £ Indicates the position 
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vector. For a given c. the function R xy (c,x) Is always a real-valued function 
which may be either positive or negative. Furthermore, R xy (c,x) does not 
necessarily have a maximum at t = 0, nor Is R xy (c,x) an even function as was 
true for the one-point auto-correlation functions. However, R xy (c,x) does 
display the symmetric relation (Bendat and Plersol, 1971): 

Rxy(C.-r) ■ R yx (c,T) (2.14) 

where x and y are Interchanged. 

The normalized cross-correlation function Is then defined as: 


B xy (c,T) 


Rxy(c»t) 

a x a y 


(2.15) 


where a x and a y are the standard deviations of x(j*,t) and y(£+c,t), 
respectively. At c = 0, Rxy(0 4nd B xy (x) are called the one-point cross- 
correlation and the normalized one-point cross-correlation functions, 
respectively. 

Figure 2.10 shows typical normalized one- and two-point cross- 
correlation functions for Run 3 In Flight 31. The upper half of the figure 
shows the one-point cross -correlation coefficients for three combinations of 
the turbulent velocity components measured with respect to the right wing tip 
of the aircraft. The lower half of the figure shows the two-point cross- 
correlation coefficients for three corresponding combinations of the turbulent 
velocity components measured from the right and left wing tips. Since the 
wing span Is much smaller than the characteristic length scale of the 
turbulence, the two-point cross-correlation coefficients are quite similar to 
those of the one-point cross-correlation for all runs. 

The cross-correlation coefficients, shown In Appendix A, are, with the 
exception of Run 10, generally small and almost constant with spatial lags. 
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Figure 2.10. Normalized one- and two-point cross-correlation functions of 
gust velocities. Flight 31, Run 3. 








The cross-correlations for Run 10 Increases with spatial lag having a maximum 
at c * 1500 m. Inspection of the time histories suggests a wave phenomenon at 
the high altitude at which Run 10 was made. The first wave occurs at roughly 
17 seconds which corresponds to c * Vt of approximately 1700 m. The u and v 
components are approximately 180* out of phase resulting in a strong 
cross-correlation at q = 1500 m. 

In contrast to Run 10,. the cross-correlctions for the other runs have 
values on the order of 0.5 or less but show no pronounced peak. The value of 
these correlations (i.e., approximately 0.5) are higher than exp-, ted but the 
high values may be due to the short time records. For best results, the 
cross-correlation function of Equation 2.13, requires the sample record length 
to approach infinity. However, this Is not well approximated for several of 
the runs. The cross-correlation coefficient for Run 3, which has the longest 
sample record is very small and Is expected to be the best representative of 
the true cross-correlation coefficients. 

2.6 Normalized One- and Two-Point Spectral Density Functi ons of Gust 
Velocities 

The sixth part of the analysis for each run In Appendix A Is the 
spectral analyses of the turbulence velocity components. The spectral 
analysis Includes the normalized one- and two-point spectral density functions 
of the gust velocity components measured with respect to the right and left 
wing tips of the aircraft. In addition, the auto-spectrum and the one-point 
cross-spectrum are compared with predictions from theoretical models. The 
spectra presented in this report are one-sided spectra (see Bendat and 
Plersol, 1971). 

The definition of spectral density functions In terms of Fourier 

transforms of the correlation functions yields two-sided spectral density 
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functions which are defined for both positive and negative frequencies (—,-) 
and are denoted by S(f). Assume that the auto- and cross-correlation 
functions R x (c,t) and R X y(c,x) exist, as defined In Equations 2.5 and 2.13. 
At c - zero, R x (t) and R xy (x) represent the one-point auto- and cross- 
correlation functions, respectively. The two-sided auto- and cross-spectral 
density functions are given by: 


Sx(c.f) - J R x (?.T)e-J 2 * fT dx 


(2.16) 


s xy (c,f) = J R X y(c,x)e~ J2 * fT dx 


(2.17) 


respectively. The letter j denotes the Imaginary number j » and f Is the 


frequency. 


From the symmetry properties of the correlation functions. It follows 


that: 


s x("0 = S x (f); z = 0 
S x (c»-f) * S x (c,f); c * 0 


(2.18) 


S X y(C,— f) = sjy(c,f) = Sy X (c,f) 

where designates the complex conjugate. These equations state that the 
two-sided one-point auto-spectral density functions are real, non-negative, 
and even functions of f, whereas the two-sided two-point spectral density 
functions and the one- and two-point cross-spectral density functions are 
complex-valued functions of f. 

The one-sided spectral density functions, * x (c,f) and * xy (c,f) where f 
varies only over the frequency range (0,-) are defined by: 
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* 2S x (c,f) 0 < f < - otherwise zero 

*xy(c.f) « 2S xy (c.f) 0 < f < - otherwise zero (2.19) 

In terms of the correlation function, the one-sided one-point auto-spectral 
density function becomes: 

* x (f) * 4 J r x( t ) cos 2xfx dx 0 < f < - (2.20) 

0 .. .. . .. .. 

The one-sided two-point spectral density function Is: 

* x (c,f) = 2 J R x (c.T)e _j2xfT dx 

3 C x (c,f) - JQ x (c.f) (2.21) 

and the one-sided one- and two-point cross-spectral density function Is: 

«• 

*xy(H.f) ■ 2 J R xy (c,x)e" J2xfT dx 

3 C X y(c,f) - JQ X y(c*f) (2.22) 

where C x (c,f) and C xy (£,f) are called the coincident spectral density 
functions (co-spectrum) and are even functions of f, and where Q x (c,f) and 
Qxy(c.f) are called the quadrature spectral density functions (quad-spectrum) 
and are odd functions of f. An alternative way to describe the complex- 
valued spectral density functions Is with the polar form, *(c,f) * 
i#(c.f)ie 0 (£»O, defined In terms of an absolute magnitude and a phase angle: 

|*x(c.O| = Jc x 2(c,f) + Q x 2 (c,f) , ©x(c,f) » tan-1 — - (2.23) 

C x (c,f) 

and 
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(2.24) 


- Jc xy 2(c,f) + Q xy 2(c.f) , 0 xy (c,f) « tan-1 Qxy(C,f) _ 

C X y(c*f) 

The magnitude of the comp lex- valued spectral density functions 
represents the energy associated with fluctuations at specific frequencies 
within the turbulent flow. For each run of Flight 31. only the magnitude of 
the normalized spectral density functjons . are presented In Appendix A. A 
segment-averaging technique was used to compute and smooth the raw spectral 
estimates obtained from the direct Fourier transform of the individual data 

segments. Data smoothing procedures for two-point spectra are discussed in 
Section 3. 

Figure 2.11 shows a typical plot of the normalized spectral density 
functions for Run 3 In Flight 31. The upper half of Figure 2.11 shows the 
normalized one-sided one-point auto-spectral density functions for the three 
respective turbulence velocity components measured at the right wing tip. The 
theoretical von Karman spectral density functions are also plotted for 
comparisons. The comparisons show good agreement between the experimental 
results and those predicted by the theoretical models. The Integral length 
scales, L u , L v , and L w , which were computed from the longitudinal, lateral, 
and vertical correlation functions, respectively, were used In the theoretical 
models. The experimental data Is higher than the von Karman predictions at 
high frequencies. The spectra have been corrected for variance error but not 
for aliasing nor bias error. These effects are discussed In Section 3. 

The lower half of Figure 2.11 shows the normalized one-sided two-point 

spectral density functions for three respective turbulence velocity components 

measured at the right and left wing tips of Run 3. Also plotted for 

comparison are the Houbolt and Sen tneoretlcal models. This two-point 

spectrum model Is derived from the Fourier transform of Equations 2.10 and 
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Figure 2.11. Comparison of normalized one- and two-point spectral density 

functions for gust velocities with theoretical models. Flight 31, 



2.11 where Bi_(c) and By(c) are the von Karman longitudinal and transverse 
correlations (Equation 2.8), respectively. The form of the theoretical 
spectra is: 




Co u z 



K 5 /6(Z) - (s/Lu)5/3 Z l/6 k 1/6 (Z) 


(2.25) 


for the longitudinal component, and 


*v(s,f) “ Ca v 2 


s_ 

IlJ 


5/3 


(H 


11/3 


|f^ K5/ 6 (Z) ' a 2 z ll/6 Kll/6(Z) + ^ s/L v) 5/3 Z 1/6 Kl/6(Z)j 

(2.26) 

for the lateral component. For the vertical component, the Spectrum Is given 
by: 

,5/3 f si 11/3 

U 


*w( s » f ) = Co w“ 


l-W 

V 


8 iy 


Z 5/6 


k 5/6( z ) " 


a 2 zll/6 


k 11/6( z ) 


(2.27) 


where 



s 

1 + 

2*fL 
a 

La J 


V J 


a = 1.329 


where L is any Integral length scale of L u , L v , and L w , V is the mean 
airspeed, and K Is a modified Bessel function of the second kind. 

Equations 2.25, 2.26, and 2.27 are plotted In Figures 2.12, 2.13, and 2.14 
for different s/L values, respectively. The termination of the curves In 
Figure 2.12 for the two-point spectra for the longitudinal velocity component 
are not arbitrary. At this point, the spectra based on Equation 2.25 takes on 
negative values. Negative values occur when the last term In the brackets of 
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Equation 2.25 becomes negat1ve._ This corresponds to the ratio 2K 5/6 (I,/( 2 

Kl/5(Z)) becoming unity. A piot of this ratio versus wave ouster for 

different values of s/l is shown in figure 2.15. The value of wave nuster for 

which the ratio becomes unity is indicated by the solid circles, inspection 

of this figure shows that these values nf ws.,^ 

values of wave nuiroer correspond to 

wavelengths of approximately one half the separation distance s. - " 

The idea of a negative spectrum is initially inconsistent with one's 
noraal thinking. However, after further consideration, it is totally 
consistent with physical reasoning that the energy contained in fluctuations 
of wavelengths smaller than one-half the separation distance s would be zero 
or even negative, it is also very likel, that values of the lateral and 
vertical spectrum shown figures 2.13 ,„ d 2.14 should be truncated at 
corresponding values of reduced frequency for which the longitudinal spectrum 
is truncated. The energy contained in eddies of size smaller than the 

separation distance predicted by the model for the lateral and vertical 
spectra Is not likely meaningful. 

Equations 2.25 through 2.27 were derived analytically in this study. 
Campbell (1984) obtained results consistent with those shown in Figures 2.12 
2.13. and 2.14 by numerically integrating equations similar to Equations 2.10 
and 2.11. Campbell obtained negative values for the longitudinal spectrum, 
however, he erroneously contributed them to round-off errors in his numerical’ 
integration (Campbell. 1986). further work is needed to fully resolve the 
meaning of the two-point spectrum at high frequencies for which It becomes 
negative. However, it is believed that it is consistent with physical 
principles to simply truncate the two-point spectrum at these negative values. 

Returning to a consideration of Figure 2.11, the two-point spectral 
density functions, calculated from the experimental data, are consistent with 
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qua Lion 2.25) becomes negative. 





the theoretical model until a mid-range frequency value. Above that 
frequency, the theoretical model drops off rapidly compared with the 
experimental data. Since the spectral estimate is calculated from a digitized 
turbulence time history with a finite record length, the departure of the data 
from the theory is due to aliasing and truncation error. This will be 
discussed further In Section 3. 

The normalized one- and two-point cross-spectral density functions for 
three combinations of the three respective turbulence velocity components 
measured at the right and left wing tips in Run 3 are shown In Figure 2.16. 
As mentioned earlier, the spectral density functions plotted in this report 
are magnitude only. The shape of the normalized magnitude of the spectral 
density functions are observed to be very similar for both the auto-spectral 
density functions and the cross-spectral density functions. This would not 
intuitively be expected because the normalized auto-correlation and cross- 
correlation functions, Figures 2.9 and 2.10 respectively, are very different. 
However, the non-normal ized or absolute value of the auto-spectral density 
function (i.e., $ as contrasted to <j>/o 2 ) Is about one order of magnitude 
larger than the value of the cross-spectral density function. The above 
observations suggest that the eddies of a given frequency contain 
cross-component energy proportionate to the distribution of common component 
energy; however, the cross-components are out of phase and have little 
correlation. Phase relationships for the spectra are dictated by the 
magnitude or quadspectra (see Equations 2.23 and 2.24). The quadspectra for 
the two-point common component spectra are very small; however, for the 
cross-spectra, both one-point and two-point, the quadspectra are of the same 
magnitude as the co-spectra. The former result indicates little phase shift 
between common-components displaced spatially whereas the latter result 
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indicates significant phase shift between uncommon components regardless of 
spatial displacement. - - 

Similar analyses on the turbulence velocities gathered by the NASA B-57B 
aircraft (Frost, et al., 1985a) for the two-point common component spectra 
have shown that the quadrature spectra have values near zero (<0.1) for all 
frequencies. Emphasizing that the phase shift between the same turbulence 
components measured at the different probes are negligible. Again, “this is ’ 
probably because the wing span of the airplane is much smaller than a 
characteristic length scale and significant phase shift would not occur for 
most of the turbulent eddy sizes involved. Therefore, the shapes of the 

one-point and two-point correlation functions will be similar to each other. 
However, as noted for the one- and two-point cross-spectra, the quadrature 
spectra are comparable with the corresponding coincident spectra, which means 
the phase angles are significant. Thus, when utilizing turbulence cross- 

spectral functions to assess the influence of the gust gradient on an 
aircraft's response, the phase angle of the spectral function is an Important 
parameter. 

There is little information on theoretical or empirical models for 
one-point cross-spectra and virtually no information on two-point cross- 
spectra in the literature. Reeves, et al. (1974) suggests the following 
two-sided one-point cross-spectral density function to relate the u and w gust 
components at low altitudes: 
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where r. A, and B are arbitrary parameters satisfying the inequalities: 
r > 0, 
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In this equation, r is the parameter defined in Equation 2.4. The model 
developed by Reeves, et al. (1974) is basically for low-altitude applications. 
Therefore, the lowest flight level was chosen for comparison of the experiment 
with the tneory. Inspection of the flight altitude of all runs of Flight 31 
shows that Run 12 was flown at the lowest average altitude of approximately 
400 ft above the ground. The normalized one-point cross-spectral density 
function, $ uw (f)/a u a w , for the turbulence velocity components u and w at the 
right wing tip in Run 12 was calculated from Equation 2.28 and the results are 


shown in Figure 2.17. 

The parameters A and B for Run 12 are chosen to best fit the 
experimental data. The value of r is taken as 1.25 the average value for u 
and w tabulated in Table 2.7. Figure 2.17 shows good agreement between the 
experimental results (symbol "x") and the results predicted by the Reeves 
model (solid line) except In the high-frequency regions. As discussed earlier 
for the one-point auto-spectrum, the higher experimental values in the 
high-frequency region are probably due to aliasing and truncation error. 
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2. 7 List of All Parameters Measured and the Range of Their Extreme and 

Average Values - — ~ ~ 

Finally in Appendix A, a table is presented for each" run which lists all 
parameters recorded during a flight: the units, the maximum and minimum 

values, the mean value, the root mean square value, and the number of data 
points for each parameter. These parameters are stored on magnetic tapes In 
the following order (the symbols used to represent the variable are appended 
in brackets): 

1. Mountain Daylight Time (MDT) In seconds for each record [t]. 

2. Roll rate measured by body-mounted roll-rate transducer (positive 
with right wing going down), rad/s [d<|>/dt]. 

3. Normal acceleration at c.g. (positive up) g units [a n ]. 

4. Pitch rate measured by body-mounted pitch-rate transducer (positive 
with nose going up), rad/s [de/dt]. 

5. Pitch attitude measured In the vertical plane (positive with nose 

up), rad [0]. - 

6. Roll attitude of airplane with reference to horizontal (positive 
with right wing down), rad [$J. 

7. Airplane heading measured In a horizontal plane clockwise from true 
north (always positive), 0 to 360* range H»j]. 

8. Sensitive airplane heading obtained from with arbitrary zero at 
the instant the data switch is turned on (positive with nose rightl 

±15* range 3 ' 

9. Airplane heading measured in a horizontal plane clockwise from true 
north (always positive) 240* to 600* range [^J. 

10. Sensitive airplane- heading obtained from i|^, with arbitrary zero at 
the Instant the data switch Is turned on (positive with nose riahtl 
±15* range [A^]. 

11. Normal acceleration at the left wing tip (positive up), g units 
UnlJ' 

12. Normal acceleration at the right wing tip (positive up), g units 
l a nRJ* 
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13. Longitudinal acceleration at the c.g. (positive forward), g units 
[a x ]. 

14. Lateral acceleration at the c.g. (positive toward right wing), g 
units [ay]. 

15. Angle of attack measured at the airplane nose boom (positive with 
flow vane trailing edge up), rad [oq]. 

16. Angle of sideslip measured at the airplane nose boca (positive with 
flow vane trailing edge toward right as viewed frca the aircraft), 
rad [Bel- 

17. Temperature of the INS pallet, *F [Tj], 

18. Temperature of the instrument pallet, *F [Tp], 

19. Vertical acceleration of the INS stable element (positive up), g 
units [a z ]. 

20. Angle of attack measured at the right wing tip boca (positive with 
flow vane trailing edge up), rad [or]. 

21. Angle of sideslip measured at right wing tip boom (positive with 
flow vane trailing edge toward right as viewed frea the aircraft), 
rad [8r]. 

22. Angle of attack measured at the left wing tip booa (positive with 
flow vane trailing edge up), rad [cqj. 

23. Angle of sideslip measured at the left wing tip boon (positive with 
flow vane trailing edge toward right as viewed from the aircraft), 
rad [Bj* 

24. Yaw rate measured by a body-mounted yaw-rate transducer (positive 
with nose going right), rad/s [dtp/dt]. 

25. Total temperature, *C [T 0 ]. 

26. Impact pressure measured at the left wing tip boom, P a [qcL^ - 

27. Impact pressure measured at the airplane nose boom, P a [qccl* 

28. Impact pressure measured at the right wing tip boom, P a [qcR]. 

29. Free-stream static pressure measured at the airplane nose boom, KP a 

[p]. 

30. Temperature of IRT, *C [Tjrj]. 

31. Wet-bulb temperature, *C [T wb ]. 


63 



32. Turbulent fluctuation of Impact pressure at the left wing tip boca. 

Pa CflctL-]- - 

33. T’lrbulent fluctuatlor of Impact pressure at the airplane nose boom, 
f-a Cflctcl. 

34. Turbulent fluctuation of Impact pressure at the right wing tip 
boom, P a [QctR]« 

35. Deflection angle of the right aileron, deg [5 a R]. 

36. Deflection angle of the left aileron, deg [4 a [_]. - 

37. Deflection angle of the elevator, deg [6 e ]. 

38. Deflection angle of the stabilizer, deg [5 S ], 

39. Deflection angle of the rudder, deg [6 r ]. 

40. Thrust ratio of the right engine to the maximum thrust, percent 
CTr]. 

41. Thrust ratio of the left engine to the maximum thrust, percent 
[T L J. 

42. Deflection position of the flap system [6f]. 

43. Deflection position of the speed brake system [6 S bJ. 

44. Distance to go from the present position of the aircraft to the 
next waypoint set on the INS (always positive), m [At]. 

45. Bearing to destination, i.e., bearing from the aircraft's present 
position to the next waypoint set on the INS, measured In a 
horizontal plane clockwise from true north (always positive), deg 
C YB i * 

46. Longitude of aircraft as measured by INS, deg [LONG]. 

47. Latitude of aircraft as measured by INS, deg [LAT]. 

48. Track angle of airplane measured In a horizontal plane clockwise 
from true north (always positive), deg [yr). 

49. Airplane heading, measured in a horizontal plane clockwise from 
true north, rad [j»], 

50. East-west component of the airplane Inertial velocity as measured 
by INS (positive toward east), m/s [Vg]. 

51. North-south component of the airplane inertial velocity as measured 
by the INS (positive toward north), m/s [V^]. 
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52. 


53 . 

54 . 

55 . 

56 . 
57 - 60 . 

61 . 


62 . 


63 . 

64 . 


65 . 


66 . 

67 . 

68 . 

69 . 

70 . 

71 . 

72 . 

73 . 


Pressure-derived altitude based on standard atmosphere tables, km 
[h p ]. _ 

Computed free-stream temperature, *C [T c ]. 

Computed east-west wind component (positive toward east), knots 

[W E 1. 

Computed north-south wind component (positive toward north), knots 

[W N ]. 

Computed magnitude of wind vector, knots [WJ. 


Computed direction of wind vector, deg [*w3* 

True airspeed computed from Impact pressure measurement at right 
wind tip boom, m/s [Vr]. 

True airspeed computed from impact pressure measurement at the 
airplane nose boom, m/s [Vq]. 

True airspeed computed from the Impact pressure measurement at left 
wing tip boom, m/s [Vl). 

Incremental pressure-derived altitude with reference to value at 
beginning of run (positive when altitude Increases), m [Ah p j. 


Computed corrected inertlaldlsplacement, m [Ah c ]. 


Computed longitudinal component of gust velocity at right wing tip 
boom (positive In direction of flight path), ra/s [ur]. 


Computed longitudinal component of gust velocity at airplane 
centerline nose boom (positive In direction o. flight path), m/- 

Cur]. 


Computed longitudinal component of gust velocity at left wing tip 
boom (positive in direction of flight path), m/s [U[_]. 


Computed lateral component of gust velocity at right wing tip 
(positive toward right), m/s [vr). 

Computed lateral component of jjust velocity at airplane centerline 
nose boom (positive toward right), m/s [v c ]. 


Computed lateral component of gust velocity at left wing tip boom 
(positive toward right), m/s [vl). 

Computed vertical component of gust velocity at right wing tip boom 
(positive up), m/s [wr]. 

Computed vertical component of gust velocity at airplane centerline 
nose boom (positive up), m/s [wq]. 
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74. Computed vertical component of gust velocity at left wing tip boom 
(positive up), m/s [w|_]. 
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3.0 SPECTRAL ESTIMATION 


Spectral analysis of atmospheric turbulence generally Involves the Fourier 
transform of a digitized finite-duration velocity fluctuation time history. 
The digitization process and the truncation associated with the finite time 
increment of the time history result in both aliasing and bias errors; while 
the random nature of turbulence results in variance errors. Although these 
errors cannot be totally eliminated, appropriate filtering will reduce their 
magnitudes. It was found that the magnitudes of the respective errors and the 
effects of the filtering process are quite different for two-point spectra 
than they are for one-point spectra. Because of the significance of the 

difference, a detailed discussion of the effect of aliasing, bias, and 
variance errors on the one-point and two-point spectra is given In this 
section. The magnitude of the errors are also estimated for typical data such 
as that reported in Appendix A. 

To illustrate the magnitude and nature of the errors, an analytical von 
Karman one-point correlation and the Houbolt and Sen (1972) two-point 
correlation are used to investigate aliasing and the bias error (also called 
truncation error or spectral leakage) associated with discrete Fourier 
transforms. Spectrum for each of these correlations has been computed 
analytically. The analytical spectra Is then compared graphically with the 
spectrum estimate calculated from the discrete Fourier transform (DPT) of the 
digitized analytical correlation functions. This comparison, based on 
analytical models, illustrates the aliasing and bias errors occurring simply 
from the digitization and truncation process. 

The DFT of the actual turbulence data is then computed. The resulting 
spectra which are now based on a random signal are then compared with the 

spectra calculated analytically from the digitized continuous correlations 
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functions. It is shown that the spectra calculated front the random data 
• contain not only aliasing and bias errors but also variance errors. The use 

of segment- averaging (Bendat and Piersol, 1971) to reduce the variance error 
° f the spectra computed from the random turbulence data time histories Is then 

♦ 

discussed. 

3.1 Graphical Illustration of the Discrete Fourier Transform 

The generation of errors associated with the Fourier transform of a 
digitized function can be conceptually explained by a graphical illustration. 
Following closely the development of Brigham (1974), consider some function 
h(x) and its Fourier transform H(f) illustrated in Figure 3.1a. To determine 
the Fourier transform of h(x) by means of digital analysis techniques, it Is 
necessary to digitize h(-r) at discrete increments in time. 

Digitizing h(x) in Increments of Ax Is equivalent to multiplying it by the 
comb function shown in Figure 3.1b. The comb function, A 0 (x), has the Fourier 
transform, A 0 (f), shown in the corresponding figure. The Fourier transform of 
the product h(x)A 0 (x) Is given by the convolution integral of H(f) and A 0 (f) 
designated by H(f)*A 0 (f), i.e. 

vo 

H(f)*A 0 (f) = J H(f’)Ao(f - f')df’ ( 3 .!) 

Figure 3.1c Illustrates H(f)*A 0 (f). 

Note that the transform H(f)*A 0 (f) differs from the analytical transform 
by the appearance of images of the analytical spectrum H(f) displaced along 
the frequency axis at a spacing of til/Axi. Each of the images contribute 
some energy to the true spectrum centered about f = 0. This effect Is called 
aliasing which occurs due to working with a digitized function. 
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The Fourier transform pair in Figure 3.1c is still not suitable for 
machine computation, however, because an Infinity of digitized values of h(x) 
Is considered; It is necessary to truncate the sampled function h(x) so that 
only a finite number of points, say N, are considered. The rectangular or 
truncation function, w(x), and its Fourier transform, W(f), are illustrated in 
Figure 3. Id. The product of the infinite sequence of impulse functions 
representing h(x) and the truncation function (i.e. , h(x)A 0 (x)w(x)) yields the 
finite length time function Illustrated in Figure 3.1e. The Fourier transform 
of the truncated, digitized function is given by the convolution of H(f)*A 0 (f) 
with W(f) or [H(f )*A 0 (f )]*W(f ). 

As illustrated in Figure 3.1e, the frequency transform now has a ripple to 
it; this effect has been accentuated in the illustration for emphasis. The 
form of W(f) for the rectangular function of unit magnitude is: 

W(f) = 2T sin(2*Tf)/2*Tf (3.2) 
Hence, if the truncation (rectangular) function is Increased in length, then 
the s1n(f)/f function will approach an Impulse; the more closely the s1n(f)/f 
function approximates an impulse, the less ripple or error due to truncation 
will be Introduced by the convolution. Therefore, It is desirable to choose 
the length of the truncation function as long as possible. 

The effect of digitization and truncation on typical turbulence 
correlation/spectrum Fourier transform pairs is discussed in the following 
sections. The discussion Is presented in terms of an example for the one- and 
two-point correlations and spectra, respectively. For the one-point 
correlation, assume the function h(t) in Figure 3.1 represents the theoretical 
von Kantian transverse correlation function given by Equation 2.8 and repeated 
here in lag time domain (c = Vx) for convenience: 
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and the function H(f) corresponds to the spectrum given by an analytical 
Fourier transform of Equation 3.3, i.e.: 
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For the theoretical two-point correlation, the Houbolt and Sen (1972) 
model is used, I.e.: 

i /*} . . 

(3.5) 
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, a = 1.339 


where 


c = Js 2 + (Vx)2 

and the theoretical spectrum Is given by: 
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3.2 Aliasing 

Consider first the problem of aliasing. — Aliasing can be described by 
considering a continuous record which Is sampled such that the time interval 
between sample values Is At seconds. The sampling rate is then 1/At samples 
per second. However, at least two samples per cycle are required to define a 
frequency component in the original data as illustrated by the sketch in 






/ 


Figure 3.2. Hence, the highest frequency which can be defined by sampling at 
a rate of 1/At samples per second Is 1/2At Hz. Frequencies In the original 
data above 1/2At Hz will be folded back into the frequency range from 0 to 
1/2 At Hz, and be confused >/ith data in this lower range, as illustrated In 
Figure 3.1. The cutoff frequency f c = 1/2At Is called the Nyquist frequency 
or folding frequency. For any frequency f in the range 0 < f * f c , the higher 
frequencies which are aliased with f are given by (see Bendat and Piersol, 

1971): 

(2f c ± f), (4f c ± f), ... , (2nf c ± f), ... ( 3 - 7 ) 

To demonstrate the magnitude of aliasing, we have digitized Equation 3.3 
representing h(r) in Figure 3.1 for values of L = 300 m and V = 100 m/s and we 
have used a discrete Fourier transform (DFT) to compute H(f). The resulting 
value, H(f)0FT, is plotted along with the analytical H(f) In Figure 3.3. The 
analytical function can be computed to as high a value of frequency as desired 
by integrating Equation 3.3 mathematically to give Equation 3.4 (i.e., 
mathematically f can approach infinity). However, to employ a DFT method to 
compute the spectra from Equation 3.3 a finite record length, T, must be used. 
Spectra computed for two different values of T = 12.8 sec and T = 51.2 sec as 
contrasted to infinity for the theoretical model are shown In Figure 3.3. 
Considerable departure of the theoretical spectrum curve from the DFT 
determined spectrum curve Is observed. The reason for this departure is 
associated with both aliasing and truncation errors as discussed in the 
following. 

Since the theoretical value of H(f) is known, the turbulence energy 
aliased into a given frequency f can be computed by Inputting the values from 
Equation 3.7 into Equation 3.4, i.e.: 




— — -- 7 * 
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(3.8) 


Aliased energy at frequency f « E H(2nf c ± f) 

n=0 

Table 3.1 Illustrates the magnitude of aliasing for the case f * 15 Hz and 20 

Hz, respectively. The values In the table are obtained by summing Equation 

3.8 for 0 < n < 20 with H(f) given by $ w (f) from Equation 3.4. 

TABLE 3.1. Comparison of Aliased Spectrum Values with True Analytical Values 
at f - 15 Hz and 20 Hz (see Figure 3.4). 



True theoretical value 0.8075 x 10"3 0.5000 x 10~3 

Aliased val ue 1.579 x 10~3 1.424 x 10'3 


Figure 3.4 is an enlargement of Figure 3.3 for the 1 to 100 Hz frequency 
range. (The variance error bars shown on the figure are described later.) 
The aliased values for the 15 and 20 Hz frequencies In Table 3.1 are plotted 
on the figure. The plot clearly shows that the major portion of the 
difference between the theoretical curve and the DFT curves is due to 
aliasing. This Is evident from the fact that when the aliased energy Is added 
to the theoretical curve, the results almost coincided with the H(f) DFT 
functions. There remains a small difference which is attributed to bias 

error. Note that this small difference decreases with increasing record 
length. 

Table 3.2 shows a similar aliasing calculation for the two-point spectra 
(Equation 3.6). The energy aliased Into the 15 Hz frequency is almost zero 
and In the 20 Hz frequency Is only doubled. 

Figure 3.5 is a comparison of the two-point spectrum computed from a DFT 


of Equation 3.5 digitized at At = 0.025 sec with the analytical value from 
Equation 3.6. (The solid circles and variance bars are described later.) 





Frequency (Hz) 


Figure 3.3. The theoretical von Karman one-point auto spectrum compared 
with the DFT computed values (triangular spectrum window). 






(sec) 



Figure 3.4. Comparison of aliased one-point auto-spectrum with true 
analytical values. 



Figure 3.5. 


Comparison of theoretical two-point spectrum (Equation 3.6) 
with the OFT of the correlation function (Equation 3.5). 
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Significant difference between the analytical H(f) and the DFT, H(f)oFT* is 
observed. The difference In this case, however. Is not due to aliasing as Is 
evident from the aliased values plotted on the figure. The significant 
departure of the DFT value from the analytical value is due to bias error. 

TABLE 3.2. Comparison of Aliased Values for Two-Point Spectra with True 
Analytical Value. 

f « 15 Hz f » 20 Hz 

True analytical value 0.1764 x 10-1° 0.2239 x 10“13 

Aliased value 0.1764 x 10-1° 0 .4478 x 10-13 




3.3 Truncation Error 

g-fas error occurs due to truncation of the time history and appears as 
ripples In the DFT curve as Illustrated in Figure 3.1e. Bias error and 
truncation error are in effect the same thing and are referred to 
Interchangeably throughout the remainder of this Chapter. Bias error Is 
Influenced by the use of the lag window w(t). In Figure 3. Id the function 
w(t) Is a rectangular lag window. 

The spectra shown In Figures 3.3 and 3.5 have beer corrected for bias 
errors with a triangular or Bartlett lag window, w(-r), defined as: 


w(t) * 1 * * ,T ' < 21 (3 ‘ 9) 

The Fourier transform of w(x) Is called the spectral window, W(w). The 
Bartlett spectral window Is given by: 
sln^ Tu> 


W(«) 


xTuj2 


( 3 . 10 ) 


Bias or truncation errors occur when attempting to compute the power spectrum 
S(w) of a real, stationary process, x(t), from a single realization of x(t) 
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available only over a finite Interval (-T,T). To demonstrate the generation 
of bias error when computing turbulence spectra, consider first the 
calculation of the spectrum from the Fourier transform of the R(x) function. 
This development closely follows the excellent presentation In Papoulls, 1977. 
We will use as an estimate of R(x) the time average: 

T- i x 1 /2 

R(x)-|y J x(t + J] x(t - f]dt (3.11) 

-T+ixi/2 

where the function Is defined as above for ixi < 2T; and, for ixi > 2T, It Is 
assumed to be zero. The estimate of the spectruin, S(w), Is then given by the 
transform of R(x): 

2T _ 

i(o>) - { R(x)e" JuT dx (3.12) 

-2T 

It follows from Equation 3.11 that the expected value of R(x) In terms of the 
true R(x) Is: 

E{R(x)} * R(x) (l - ^f-)p2T(x) = R(x)q 2T (x) (3.13) 

where P 2 t(t) and q2T( T ) are a pulse and a triangle lag wlno.w, respectively 
(see Figure 3.6). Thus the estimated autocorrelation function Is biased, 
because Its mean Is not the true autocorrelation function at lag x. We can 
say, however, that this R(x) is asymptotically unbiased because the ixi/2T 
term vanishes at T •» -. We could easily get unbiased estimates by using 2T - 
ixi to divide the Integral In Equation 3.11 rather than 2T. The form used 
here Is preferred for the reason that R(x) In this form can be expressed In 
terms of the given x(t) as a convolution. 

R(X) = 2f x(t) * x( ' t) 
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Figure 3.6. Rectangular and triangular lag window and their respective 
Fourier transform. 




From Equation 3.12 the expected value of the spectrum Is: 

2T 

E{S(«)} - J R(T)02T(T)e' jt0T dx - S(«) * - sin2 & (3.14) 
-2T * TuZ 

where designates the convolution Integration as defined In Equation 3.1. 
Thus, tne estimate S(w) equals the convolution of S(u») with the kernel 
s1n2(T(i))/xT<i)2 (l.e., Bartlett window). The estimator Is therefore biased. 
However, 

E{S(u>)} S(u) (3.15) 

The effect of a rectangular versus triangular lag window Is Illustrated In 
the following. Figure 3.7 shows the spectrum computed from a discrete Fourier 
transform of Equation 3.3 without a lag window (or In other words with a 
rectangular lag window, see Figure 3.6). Figure 3.7 should be compared with 
the spectrum In Figure 3.3 which was computed with a triangular lag window. 
The two curves are not appreciably different In terms of departure from the 
analytical spectrum except they do show some dissimilarity In shape near the 
cutoff frequency. 

It should be observed that the magnitude of aliasing for the spectrum In 
Figure 3.7 Is of the same magnitude as that given In Table 3.1, for Figure 
3.3. One may conclude then f or the one-point spectrum the primary source of 
error Is due to aliasing and that the bias errors are negligible. 

For the two-point spectrum, however, one cannot draw this conclusion. 
Figure 3.8 shows the spectrum computed from Equation 3.6 by the DFT using a 
rectangular lag window. This figure should be compared with Figure 3.5 which 
shows the two-point spectrum computed with a triangular window. Notice that 
the curves behave considerably different In Figure 3.8 than they do In Figure 
3.5. A ripple In the spectrum curve at high frequencies appears In Figure 
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Frequency (Hz) 


Figure 3.7. Comparison of the theoretical von Karman one-point auto- 
spectrum with the DFT computed values (square lag window) 




Figure 3.8. Comparison of the analytical two-point spectrum (Equation 3.6) 
with the DFT values computed from truncating Equation 3 5 with 
a rectangular lag window. 




3.8. This Is the same ripple as Indicated In Figure 3.1e due to the 

truncation or application of the rectangular lag window to h(t). Thus the 
graphical Illustration In Figure 3.1e clearly reveals how the rectangular 
spectrum window W(u) convolved with the theoretical spectrum results In a 
truncation or bias error shown In Figure 3.8. Employing the triangular lag 
window results in a much smoother spectral curve (Figure 3.5) but there is 
still appreciable error between the DFT and the theoretically computed 
spectra. 

The order of magnitude of aliasing associated with digitization of 
Equation 3.5 is the same as given In Table 3.2. The aliasing error is very 
small. It can, therefore, be concluded that the bias error contributes 
significantly to the error produced from computing two-point spectra with 
discrete Fourier transforms. The fact that the bias error for the one-point 
spectrum Is small compared to aliasing whereas the bias error for the 
two-point spectra is very large compared to aliasing Is a very striking 
observation when one considers that the correlations functions from which 
these two spectrum are computed differ only by just a miniscule amount near 
the zero lag value (see Figure 3.9). 

Figure 3.9 Illustrates the very small difference between the one-point 
correlation function and the two-point correlation function from which the 
spectra In Figures 3.3, 3.5, 3.7, and 3.8 were computed. A very small change 
In the correlation near zero lag causes the very large differences observed in 
the spectra. The effects of. bias error generally not significant for 
autospectra must, however, be carefully considered when computing or 
Interpreting two-point spectra. 

In addition to abiaslng and bias error, there Is a third form of error 
called variance error. Variance error occurs when computing spectra and 
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correlations from random data. The analyses described to this point have 
dealt with deterministic correlation functions which do not cause any variance 
error In the spectral calculation. The aliasing and truncation error 
discussed occur because of the digitization of the function and because of the 
finite length of the record required for application of the DPT. The variance 
error occur because of the random nature of the gust velocities. In 
attempting to minimize both the bias and variance error, there are conflicting 
requirements on record length. These requirements result in accepting certain 
tradeoffs between resolution and accuracy. A more detailed look at the bias 
and variance errors will shed some light on these tradeoffs. 

Consider initially the method of computing the spectrum by first 
calculating the correlation R(t) and then Fourier transforming R(t) to obtain 
*(f) as we have done for the deterministic functions in the previous 
discussions. It follows from Equation 3.14 and 3.15 that If T Is sufficiently 
large, then: 

i{S(u)} - S(u>) ( 3 - 16 ) 
for analytical functions. 

However, even for analytical functions with large T the variance of S(w) will 
not be small for random data. In fact, for any T, 

var[S(u>)] > e2{S(oj)} ( 3 - 17 ) 

Therefore, S(n>) Is not a good estimator of S(u»), no matter how large T Is. 
The reason Is that with random data the values of_the Integrand, R(t), In 
Equation 3.12 are not reliable (have large variance) for t close to ±2T. 
Thus, the power spectrum S(w) of a process x(t) cannot be determined from a 
single sample, no matter how large the sample Is. To reduce the variance of 
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“* eSt,K,te> " e ”' ust «"* a smoothed version of S(„); that Is. we must 

sacrifice resolution. 

3.4 Smoothing the Spectrum 

The variance of the Integral In Equation 3.12 can be reduced by deempha- 

altmg the contribution of R(,) for t near =2T. For this purpose, the 
estimator: 


2T 


Sw(u>) = J R(T)w(x)e~^“ T dx 

-21 


(3.18) 


Is formed where w(,) Is a lag window vanishing for i„ > 21 . In this section. 

we examine the properties of and the factors affecting the selection of 

the window w(x). 

The estimator S w (u) Is the Fourier transform of the product R(x)w(x); 
hence. 


3w(w) = ^ _J S(u - y)W(y)dy = i(ic) * W(w) 

From the above and Equation 3.14, It follows that: 

E{Sw(“>)) = ^ E{S(o,)} * W ( U ) = L. S(w) * slniTa, # ^ 

xTw? 

For a reliable estimation, the duration of W(«) must be large compared to 1/T. 
This leads to the approximation: 
sin? T<a 


(3.19) 


(3.20) 


xTw? 


* W(u>) * W(u) 


(3.21) 


Inserting this approximate expression Into Equation 3.20 gives: 

E{SwWJ * 4r S(w) * W(w) 

2 * (3.22) 

It can be shown that under certain general conditions (Papoulls, 1977) the 
variance of ^(w) Is given by: 
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a 


varfS^u))] =. ^ s2( u ) 
where 


cj * 0 


2T 


E w - J w2( T )d T - J W2(w)du) 

-2T 


(3.23) 


(3.24) 


Equations 3.22 and 3.23 dictate the factors affecting the selection of the 
window pair w(t) and W (w). For the bias error, 

b = E{S*( tt )} - S(m) {3 25) 

to be small. W(u>) must be of short duration. For the variance to be small, E w 
must be small. We shall presently see that if T is sufficiently large, then 
both requirements can be reasonably satisfied. 

3.5 Window Selection 

For a satisfactory estimation of S(w). the variance of the estimator S^u) 
must be small compared to $2( u ) or, equivalently, the variance ratio 

3 =, varfSw(u>)l 
S2(«) 


(3.26) 


must be very small compared to 1: 

S <<C 1 (3.27) 

This Is the case If E w (see Equation 3.23) Is very small compared to 2T: 

E w » 2TS « 2T 

The above requirement leads to the conclusion that w( T ) must take on 
significant values only In an Interval (-M.M) such thatM « 2T. We shall 

assume that iw( T )i < 1 for all t and that, for i t i > M, It Is not just small 
but it vanishes: 

" (t) * 0 for 111 > M (3.28) 

From these assumptions, it follows that 
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E w < 2M so that 8 <, M/T (3.2S) 

Thus, to satisfy the variance requirement (Equation 3.27), we must choose 
M such that 

M « T (3.30) 

With M so determined, the shape of the window Is selected so as to minimize 
the bias 


b = J SC 01 * - y)W(y)dy - S(u>) (3.31) 

The bias b depends not only on W(u>) but also on he shape of S(oj). Therefore, 
there Is no well-defined optimum window. However, If T Is sufficiently large 
and W(w) > 0, It can be shown (Papoulls, 1977) that: 

co 

b „ J ^(ajjdw (3.32) 

— an 

The problem Is to find a positive function W(w) for a specified value E w 
(Equation 3.24) which minimizes the integral In Equation 3.31 or 3.32. We 
will now consider some estimates of the bias and variance errors which may be 
expected to occur In the computation of turbulence spectra. 


3.6 Variance Error 

The variance error associated with atmospheric turbulence spectra which 
follow Equations 3.4 and 3.6 can be estimated from Equations 3.23 and 3.24. 
The lag window for these figures Is w(t) - 1 - iTi/2T ra where T m - T/M. Hence, 




(3.33) 


Substituting into Equation 3.23 gives: 
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(3.34) 


vartSv/w)] » S 2 (u) 

The magnitude of the variance errors estimated from Equation 3.34 is Indicated 
on Figures 3.4 and 3.5 by the variance error bars. The variance Is computed 
from Equation 3.34 at f = 15 Hz and 20 Hz and at f = 10 Hz and 15 Hz, for 
Figures 3.4 and 3.5, respectively, with M = 5 which is the typical number of 
segments used In computing the spectra shown in Appendix A. It should be 
noted that the spectra shown In Figures 3.3 and 3.5 have no actual variance 
error since they are computed from deterministic functions. However, the 
error bars do Indicate the one standard deviation error which can be expected 
In analyses of random turbulence signals which physically obey the analytical 
Equations 3.4 and 3.6. 

Now consider the magnitude of the bias error. Equation 3.31 can be used 
to compute the bias error. Values for the convolution of S(w) (given by 
Equation 3.4 for the one-point spectrum and by Equation 3.6 for the two-point 
spectrum, respectively) with W(w) for a triangular spectral window (see Figure 
3.6) are given In Table 3.3 for frequency values of 10, 15, and 20 Hz. The 
values of the convolution Integral which Includes bias error given in the 
table for the two-point spectrum are plotted on Figure 3.5 (as marked by the 
solid circles). They coincide very closely with the DFT curves clearly 
illustrating that the departure of the DFT curves from the theoretical curve 
for the two-point spectral is due almost entirely to bias error. 

The Important conclusion from Table 3.3 Is that a two-point spectrum 
computed with truncated digitized turbulence data with no prior knowledge of 
the actual form of the spectrum may show appreciable energy In the high 
frequency range whereas In reality there is no energy at those frequencies. 
The bias error in the two-point spectrum is approximately 385,400 percent at 
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TABLE 3.3. Values of l/2x _J[ S(u - y)W(y)dy for S(u>) Given by Equations 

3.4 and 3.6, Respectively, and a Triangular Spectral Window 
(N = 512, Ax = 0.025). 



One-Point Spectrum 


Two-Point SDectrum 

f 

Bias Error 

True 


Bias Error 

True 

10 

0.8000 x 10-3 

1.5731 x 

10*3 

0.6286 x 10~ 4 

0.1631 x 10' 7 

15 

0.4011 x 10-3 

0.8075 x 

10-3 

0.1789 x 10~ 4 

0.1764 x 10' 10 

20 

0.2460 x 10-3 

0.5000 x 

10-3 

0.1568 x IQ’ 4 

0.2218 x 10-13 


10 Hz as contrasted to only about 50 percent for one-point spectrum. The 
remarkable phenomenon, however, is that the two spectra are computed from 
correlations which are almost identical except for a very small difference at 
zero lag (see Figure 3.9). These factors have strong implication when 
computing spectra from turbulence time histories where the true spectrum is 
not known a priori. With this in mind, improved windows for reducing bias 
errors in two-point spectra were investigated. 

Considerable study of computing spectra for random data has been carried 
out in the communication engineering field. Several alternate windows for 
smoothing the spectra estimates have been proposed. Papoulls (1977) gives the 
expression: 

w(t) = | [sin y t| + [l - -y^cos f t , 0 < i t i < T (3.35) 


and 


W(w) 


- 4t,2 (1 + COS Tu) 

(*2 - t 2 * 2)2 


(3.36) 


This spectral window is called the minimum-bias window because it minimizes 
the value of the integral in Equation 3.32. When this window is applied to 
the two-point correlation, the results shown in Figure 3.10 are achieved. The 
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= Analytical solution 

= Discrete Fourier transform 
(-12.8 sec <_ t < 12.8 sec, 
Lz 3 0.025 sec) 

= Discrete Fourier transform 
(-51.2 sec <_ t <51.2 sec, 
L*z = 0.025 sec) 
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Figure 3.10. Results of applying minimum-bias window (Equation 3.34) 
to the DFT computation of the two-point spectrum from 
Equation 3.5 as compared with the analytical value. 


bias error Is almost completely eliminated to frequencies as high as 15 Hz. 
The fluctuation in the curve at higher frequencies Is believed to be a 
computer-generated numerical roundoff error. Although this window has not 
been used In computing the turbulence spectra in this report further 
consideration of Its use Is recommended. 

Attention Is now directed toward calculation of spectra from the 
turbulence time histories. The spectra discussed will now not only contain 
aliasing and bias error but also variance error. 

3,7 Spectrum Calculati on from a Finite Turbulence Time History 

Consider the variable x(t) as a velocity fluctuation In the interval 
(~T»T). The correlation estimator is then given by: 

T-iti/2 

rm - J 


x [t + JJxjt - jjdt , i x i < 2T 

-T+iti/2 

The Integral in Equation 3.37 Is the convolution 

R(t) x'(t) * x'(-t) 

where 


(3.37) 


(3.38) 


x'(t) = x(t)p2t(t) = 


iTl < T 
iTl > T 


(3.39) 

equals the truncated time record of x(t). 

It follows from Equation 3.39 and the convolution theorem that the Fourier 
transform of R(t) gives the estimate of the spectrum: 

.2 


2T - T 

S(u>) - J R(x)e ^“ T dx = J x(t)e~^ u,t dt 

-2T -T 


(3.40) 


Thus. S( u ) can be determined directly from the given sample or turbulence time 
history x(t). This approach was used In computing all spectra In Appendix A 
and those discussed In Section 2. 
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It is interesting to consider now three numerical methods for estimating 
the power spectrum S(w) In terms of the given segment, x'(x) = x(x)p 2 t(t). 
All three methods are statistically Identical but differ only in the 
computational procedures. 

1. Determine the sample auto-correlation R(x) by convolving x'(x) with 
x'(-x) as In Equation 3.36: 

R(t) * x '(x) * x’(-x) (3.41) 

Multiply R(x) by the window w(x) and compute the Fourier transform of the 
product as in Equation 3.18 to get S^u). The required operations are one 
convolution, one multiplication, and one Fourier transform. 

2. Compute the Fourier transform of x'(x) to get X'(u). Multiply X'(u) 
by its conjugate and form the sample spectrum S(w): 

S(u.) = y |X'(c)| (3.42) 

Convolve S(u) with the window W'w) to get: 

Sw(uj) = S(u) * W(ui) (3.43) 

The required operations are one Fourier transform, one multiplication, and one 
convolution. 

3. Compute X'(uj) and S(w) as In method 2. Find the Inverse Fourier 
transform R(x) of S(w). Form the product R(x)w(x) and compute Its Fourier 
transform S^w) as in method 1. The required operations are two multiplica- 
tions, two Fourier transforms, and one inverse Fourier transform. 

Method 3 has been used for computing spectra from the turbulence time 
histories discussed in this report. In most cases, however, the lag window 
operation has not been performed for reasons described later. The variance 
error for spectrum estimates obtained by a direct Fourier transform operation 
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on the digitized time history (Method 3) may be determined differently from 
that given by Equation 3.23. 

Consider the spectrum function of a stationary (ergodic) Gaussian random 
process x'(t). An estimate of S(u>) can be obtained from Equation 3.42. The 
narrowest possible bandwidth resolution from Equation 3.42 is B e - (1/T). To 
determine the variance of the estimate of S (w), observe that the Fourier — 
transform X(io) Is defined by a series of components at frequencies f - k/T; k 
* 1»2,3, etc. Further observe that X'(<i)) Is a complex number where the real 
and imaginary parts, X' R (u>) and X'i(id), can be shown to be uncorrelated random 
variables with zero means and equal variances (Bendat and Piersol, 1971). 
Since a Fourier transformation is a linear operation, X' R (<u) and X'jtw) will 
be Gaussian random variables If x(t) Is Gaussian. The random variable x(t) is 
not strictly Gaussian as discussed In Section 2; however, the variance error 
to be described is expected to be representative of the error associated with 
atmospheric turbulence. It follows then that the quantity 

i X ' (to) 1 2 = X' R (u>) + X'i(w) (3.44) 

Is the sum of the squares of two independent Gaussian variables. It can be 
shown that each frequency component of the estimate S(u>) will have a sampling 
distribution given by 

SM x 2 2 

S(«) * ~T (3.45) 

where X 2 2 Is the chi-square variable with n = 2 degrees of freedom. 

Note that the result In Equation 3.45 Is Independent of the record length 
T, that is. Increasing the record length does not alter the distribution 
function defining the random error of the estimate. It only Increases the 
number of spectral components In the estimate. The variance error of the 
estimate Is substantial. The mean and variance of the chi-square variable are 

95 


) 




\ 



n and 2n, respectively. Thus the normalized standard error, which defines the 
variance error of the estimate Is: 

g _ var^rSfmll _ 2 

S2(«) ' n (3.46) 

For the case at hand, n = 2 so 8 = 1, which means that the standard deviation 
of the estimate Is as great as the quantity being estimated. This is an 
unacceptable error for most applications. 

In practice, the variance error of an estimate , luced by Equation 3.42 
Is reduced by smoothing the estimate in one of two ways. The first way Is to 
smooth over an ensemble of estimates. This can be done by computing 
Individual estimates from M independent sample records, x-j(t); 1 = 1,2,3, ...M 
and then averaging the M estimates at each frequency of a spectral component 
as illustrated in Figure 3.11a from Bendat and Plersol (1971). The second way 
is to smooth over frequency. This can be done by averaging together the 
results for Z contiguous spectral components In the estimate from a single 
sample record as Illustrated In Figure 3.11b. In either case, the smoothing 
technique approximates the expectation operation in Equation 3.42. 

For smoothing the spectra presented In Appendix A, ensemble averaging has 
been used throughout this study. Care must be used when ensemble averaging, 
however. Each spectrum from the M segments of the total time history Is 
generally complex. To estimate the mean (or magnitude) of the spectrum from 
the segments some authors Imply ensemble-averaging of the absolute values: 

j*xy( s »f)| * ^ I^xy,l(s,f) - JQxy,l(s,f)| (3.47) 

where N ra Is the number of segments. However, few authors deal with two-point 
spectra. Jenkins and Watts (1969) correctly defined the mean of the complex 
two-point spectral function as: 


\ 


' *•* JKT+ . 
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(3.48) 


r 1 5 ^ — MSB >*** T^ fS ^ wsr ^^ T ~rrsg^/~v ^s^-t 


l®xy( s *0 I 


1 N m i N ra 

"■ w Vl (s '° ' 3 (5 ^ Vi' 5 ' 0 


wherein the real and Imaginary parts of the spectral function are segment 
averaged separately. 


Figures 3.12 and 3.13 demonstrate the difference between the two-point 
spectrum estimations from Equations 3.47 and 3.48, respectively. Both cases 
are compared with the Houbolt and Sen theoretical model. Equation 3.48 gives 
better agreement with the theoretical model. It Is also Important to note 
that when smoothed by the Equation 3.47 technique, the two-point spectrum In 
Figure 3.12 Is almost Indistinguishable from the one-point spectra as shown In 
Figure 2.11. Notice the variance Is larger for a two-point spectrum than for 
a one-point spectrum. Also, one- and two-point cross-spectra have apparent 
higher variance levels. The Increase In variance Is due to the fact that 
variability Is introduced by two separate processes rather than one (l.e., 
averaging C X y and Q X y separately as contrasted to averaging C X y Itself which 
Is the only contribution to the one-point auto-spectrum). 


3.8 Lag Windows for Reducing Bias Error 

A number of data windows for reducing bias errors are described In the 
literature. No single window has been Identified as most appropriate for an 
atmospheric turbulence signal. A cosine tapered data window to smooth the 
data at each end of the record Is commonly used In the literature but was 
found to have no effect on spectra calculated and was not used In this report. 

Lag windows (as contrasted to data windows) are applied to the correlation 
function as defined In Case 3, page 94. Figure 3.14 shows the two-point 
spectrum for a turbulence measurement and a digitized deterministic model with 
a rectangular lag window. The shape of the spectrum from the data deviates 
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0.01 0. 1 C 1.00 10.00 100.00 

Frequency (Kz) 

Figure 3.12. Segment-averaged two-point auto-spectrum using Equation 3.47. 

(Heavy solid line designates Houbolt and Sen's theoretical 
model . ) 
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Figure 3.13. Segment-averaged two-point auto-spectrum using Equation 3.48. 

(Heavy solid line designates Houbolt and Sen's theoretical _ 
model . ) qq 




Flight 31, Run 3 



Figure 3.14. Two-point spectrum computed directly from data with 
rectangular lag window. 




significantly from the values predicted by the model even at low frequencies, 
t The minimum amplitude .window Equation 3.35 was also used with the direct FFT 

of the turbulence time history as shown In Figure 3.15. The window does not 
c correct the random data Input to the same degree It corrects the deterministic 

Input. It does, however, give a better correction of the bias error than any 
of the other windows used in this study. Further Investigation of the effect 

of lag windows on two-point spectra computed directly from the turbulence Is 
recommended . 




FI ight 31 , Run 3 



Frequency, f (Hz) 


Figure 3.15. Two-point spectrum with minimum amplitude window 
Equation 3.35. 



4.0 INSTRUMENTATION ERROR ANALYSIS 


4.1 Instrumentation Problems 

The Instrumentation platform Is the B-57B aircraft. This aircraft is a 
U.S. Air Force (USAF) version of* the English Electric Canberra and was built 
under license by the Martin Company. The B-57B, designed as a tactical 
bomber, first flew In 1954 but is no longer In use by the USAF. NASA uses the 
aircraft as a flight research tool to measure wind velocity, turbulence, 
temperature, and other properties of the atmosphere. The aircraft is equipped 
to gather data on gust gradients across the 19.5 m (64 ft) wing span. 
Characteristics and dimensions of the B-57B are given in Figure 4.1. 
Additional information about aerodynamic coefficients and stability 
derivatives can be obtained from Ringnes and Frost (1985). The 
instrumentation on the B-57B include three airspeed probes located at the nose 
section and at each wing tip. The flight angles, sideslip angle, and 
angle-of-attack are measured at the same three locations. Also, 
accelerometers are placed at both wing tips and at the center of gravity 
(c.g.) for use In studying wing vibrations. Ground speed, Euler angles and 
angular rates, acceleration components, and geographical location are provided 
by the Inertial Navigation System (INS).* Details on the Instrumentation and 
its accuracy is given in Meissner (1976), Rhyne (1980), and Murrow and Rhyne 
(1981).** 


♦Carousel IV made by AC Electronics, Division of General Motors 
Corporation. 

**The INS and pressure transducers used in the flights are different 
from reported in these references. 
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Figure 4.1. Characteristics of the B-57B 









During previous research (Chang and Frost, 1985; Frost, et a!., 1985a; 
Ringnes and Frost, 1985) using data gathered with the B-57B aircraft, various 
uncertainties in the measured wind velocities and turbulence measurements were 
traced to Instrumentation characteristics. Frost, et al. (1985a) have pointed 
out Irregularities In the total pressure measurements and postulated that 
water droplets may have occasionally been Ingested in the pitot tubes. These 
caused spikes In the turbulence spectrum at approximately 15 Hz. No similar 
spikes were observed for the data from Run 31 which Is analyzed in this 
report. Chang and Frost (1985) and Frost and Huang (1983) also noted that 
there are, in some cases, discrepancies In the calculation of the total wind 
vectors. These were attributed to problems with boom alignment and with the 
INS. In this section, the data reduction procedures of the quasi-steady wind 
vector and of the turbulence are reviewed In detail to pinpoint how 
Instrumentation errors might affect the wind measurements. The magnitude of 
the errors are estimated and methods of correcting for them suggested. 


4.2 Wind and Gust Velocity Equations 

The velocity of a moving airmass with respect to earth, in this study. 
Is obtained by vectorial ly subtracting aircraft velocity with respect to the 
air mass from aircraft velocity with respect to earth. These velocities are 
referred to as airspeed and ground speed, respectively. Since airspeed Is 
measured In a body-axis (airplane fixed) reference system. It Is necessary to 
rotate the airspeed vector Into the Inertial (earth fixed) frame of reference. 
The governing equations are derived In detail In Appendix B (see also Frost 
(1981); Crooks, et al. (1967); Houbolt, et al. (1964); Lenschow (1972); and 
Axford (1968)). The present assumptions used In the equations for removing 
the aircraft motions from the wind vector are straight and level flight 
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without large perturbations. Therefore, small angle assumptions are made for 
the roll, pitch and yaw angles, and for the angle of attack and sideslip 
angle. Furthermore, it Is assumed that the product of sines of any of the 
small angles mentioned above vanishes and the cosines of small angles are 
unity. The application of these linearized equations to computing gust 
velocities for touch-and-go flights and during excursions from level flight 
during a run (e.g.. Run 9 at 7 to 11 miles, see Figure A.41) Is discussed 
later. Based on the small angle, level flight assumptions, the following 
expressions are used for computing the horizontal wind velocity components 
from the measured parameters: 


We = Vg - Vc sin|* - Be - 
W N = V N - V C cos* - Be - 


*xC+ 


VC 

S-xC^ 


(4.1) 

(4.2) 


where * Is aircraft heading and Be ^ sideslip angle. Vc Is the true airspeed 
of the aircraft, V E and are east-west and north-south components of the 
airplane inertial velocity, and 4 x c Is the longitudinal distance measured 
parallel to the x-axIs of the airplane from the INS to the centerline 
measuring station. The higher order term containing * arise because the air- 
speed and ground speed are measured at different locations. 

Wind speed and direction are derived directly from two Independent 
components and are given by: 

W - (W E 2 + W N 2)l/2 (4.3) 




tan"l 



+ x 


for x/2 > tan" 



> - x/2 


(4.4) 
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Positive wind Is defined as a wind blowing towards the east. W E . and north, 
W N , . "hlch In meteorology Is referred to as west and south winds, respectively. 

The turbulence components are calculated In the aircraft-fixed 
coordinate system. A complete derivation of the equations has been carried 
out both by NASA Langley Research Center and by FWG Associates In the past. 
The FWG derivation Is also restated in Appendix B. The linearized equations 
for the center probe are: 

A A A 

uc * V E sin * + V N cos * - V c , . 

^ 4 . 5 ; 

vc = V C B C - Vq* + V E cos * - V N sin * + i xC * + (4>6) 

w C = V C o C - V C 9 + V az + 5. xC 9 - V C B C * (4 ?) 

where*. 9. and * are roll, pitch, and yaw angles of the aircraft. Those 
for the wing tip probes are straightforward modifications of those for the 
center probe. It Is assumed that the average pitch angle of the average 
pitch angle of the average flight path. 9 . Is zero. The caret f) symbol 

Indicates deviation from the mean value and the overbar (-) Indicates average 
value. 

4.3 Sources o f Inaccuracy In Data Reduction 

Instrumentation errors influence the quantities appearing on the 
right-hand side of Equations 4.1 and 4.2 and thus the accuracy of the computed 
wind velocities. Instrumentation errors In the INS ground speed components, 
the airspeed, and the sideslip angle have been Identified and studied. Errors 
In the yaw rate are negligible, and the yaw angle Is believed to be accurate. 

A test to verify yaw angle accuracy Is suggested since yaw angle errors could 
significantly contribute to errors In the calculation of horizontal wind. 
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Of these sources of Instrumentation errors, the most difficult to 
correct is the dynamic error in the velocity Inherent In the INS, termed the 
Schuler error to which aircraft motions contribute. All other errors can be 
removed by careful calibration. The effects on the magnitude of the measured 
wind and also turbulence calculations due to the sources of error in the 
Instrumentation are presented next. 

4.4 Inertial Velocity and Position Errors 

The accuracy of the calculations of horizontal winds depends upon the 
performance of the INS and its capability to provide correct measurements of 
the Inertial (ground) speed of the aircraft. In recent years mechanical and 
electronic advances have greatly Improved INS accuracy. However, a cumulative 
oscillation in the INS stable platform element called the Schuler drift 
effect, first pointed out in the famous paper by Schuler (1923), can be quite 
significant. Inertial navigation theory including derivation of the Schuler 
pendulum effects Is explained in many textbooks (see for example, Boxmeyer, 
1964). The Schuler error is essentially periodic with a period near that of 
an earth radius pendulum, 84.4 minutes.* The error behaves sinusoidally and 

*Huber and Bogers (1983) point out that a platform used In an airplane 
cannot strictly be kept tuned to T 0 * 84.4 minutes after takeoff since R 
(distance between the airplane and center of the earth) and g (gravitational 
acceleration) change with altitude. They propose to define T« = 84.4 minutes 
as the Schuler constant (for the earth). The actual period of oscillation 
proposed by these authors for a specific Schuler-adjusted system takes Into 
account the gravity gradient, the mass distribution In the system, and the 
centrifugal forces due to the velocity of the carrying vehicle. This Is 
called the actual oscillation period. The actual oscillation period of a 
specific Schuler-adjusted system (acceleration Insensitive system) under 
specific circumstances Is given by them as: 

T = k*2x Jrfg 

where k will always have a value between 0.5 and -. 


will thus change polarity. The error caused by a slow oscillation of the INS 
stable platform causes the two horizontal accelerometers to detect a part of 
the gravity vector. This false Indication of acceleration Is carried through 
the Integration for velocity and produces errors In the and W N values. 
Distance traveled or geographical position Is obtained from a second 
Integration of the measured accelerations. Thus the Schuler oscillations will 
create errors' In acceleration, velocity, and position/ The following 
procedures were used to estimate the velocity errors associated with Schuler 
drift. 

Position error can be computed from aircraft data during overflight of 
landmarks where exact geographical locations are known. Since acceleration, 
velocity, and position errors are all interrelated the Schuler error can 
experimentally be investigated by obtaining data on either one of the three 
parameters having a Schuler oscillation induced error. The velocity error Is 
generally small but increases with time, e.g., after several hours of 
operation It can be on the order of 3 to 5 m/s (Rhyne, 1980; Lenschow, 1983). 
The magnitude of the position errors for the IV INS used in B-747 alracraft 
reported by Weber (1975) normally are on the order of 10 nautical miles or 
less even after transatlantic flights. These errors are not critical for pure 
navigation purposes. But, when the objective Is to calculate wind velocity, 
the Schuler error can be quite Important. 

In an attempt to model In-flight Schuler error, data from Flight 63 have 
been analyzed. Specifics about the flight can be found In Table 1.1. A box 
pattern flight plan as shown In Figure 4.2 was flown sequentially at 1000 ft 
levels over Boulder, Colorado, In February 1984. Details of the flight and 
results are given In Chang and Frost (19e5). Each time the B-57B flew the leg 
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Figure 4.2, Event marker location of Bt 57R on box pattern flights. 



heading east, an event marker on the ground was activated to record tne moment 
a north-south running road lined up perpendicular to the flight path (see 
Figure 4.2). INS recorded longitude at the time of the event marker can thus 
be compared with the known longitude of the road to construct the Schuler 
position error (see Figure 4.3a). The exact latitude of the aircraft at the 
time of the event markers Is less certain. In fact. It depends upon the 
ability of the pilot to fly the Intended flight path. But, since the flight 
paths were, flown toward a fixed landmark, only small deviations In the 
latitude position of the east test runs are expected. A similar Indication of 
position errors has also been plotted for the latitude. Figure 4.3b. In both 
cases, the error appears to have a sinusoidal behavior. Curve fits by simple 
trial and error techniques are also plotted on the two figures. The curve 
fits suggest the latitude error has a 77-minute period of oscillation, and the 
longitude has an 111-minute period. The latitudinal period Is reasonably 
close to the Schuler constant of 84 minutes, but the longitudinal period does 
not conform to that for the latitude. Since longitude and latitude errors are 
two components derived from the same stable platform oscillation, equal period 
lengths differing only by a phase angle would be expected. Thus, additional 
Investigation of the discrepancy Is needed. 

Flight 66 (see Table 1.1) followed the same flight pattern as Flight 63, 
and the same technique for marking geographical position by event markers was 
used. Figure 4.4 has been constructed similarly to Figure 4.3. The dashed 
lines outline sinusoidal trends but are not represented by mathematical 
equations. The latitude oscillation In Flight 66 seems to have a period of 
approximately 110 minutes which Is similar to the longitude oscillation of 
Flight 63. The longitude error of Flight 66 contains more scatter In the 
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data, although the period seems to be of roughly the same length as the 
latitude oscillation on this flight. 

The magnitude INS position errors Identified are within a range of less 
than 15 km or 10 nautical miles. From a commercial aircraft operation 
standpoint, these errors are not a large problem, particularly in the 
proximity of an airport where other means of navigation are available in the 
proximity of an airport. However, Schuler position errors are of significance 
for wind measurements. Exact ground tracks are needed to determine terrain 
effects on turbulence such as wake regions behind mountains, etc. An error on 
the order of several kilometers can drastically distort the picture. 

The INS velocity errors are especially important in the wind 
measurements. Horizontal wind components are calculated based on Equations 
4.1 and 4.2. As will be demonstrated, the velocity errors can be of the same 
Grder of magnitude as the wind speed, which will greatly alter the calculation 
of the wind vector. An estimate of these errors is presented in Figure 4.5. 
The velocity error curves are calculated by taking the derivative of the 
position error curve fits illustrated in Figure 4.3. The magnitude of the 
velocity errors determined is within the range of that quoted in the 
literature (Rhyne, 1980; Lenschow, 1983). The influence of these errors is 
demonstrated in Sections 4.7 and 4.8. 

To further investigate the Schuler error Flights 73 and 74 were carried 
out where the aircraft was tracked by the NASA EPS-16 #34 tracking radar. The 
radar track provided the location and the ground speed of the aircraft 
throughout the flight. The investigation of Schuler velocity errors for 
Flight 73 and 74 has not been completed due to the late reception of flight 
data for Flight 73 and of the need to correct the radar tracking. However, 
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data on post-flight Schuler velocity errors recorded on the ground have been 
received from NASA/laRC along with data from Flight 74. The north-south and 
east-west velocity errors of Flight 74 and the ensuing post-flight velocity 
measurements are plotted In Figures 4.6 and 4.7. The in-flight velocity 
errors are obtained by comparing aircraft and radar data assuming the radar 
Indications are free of error. The data recorded on the ground is a direct 
measure of the Indicated velocity from the INS while the aircraft was parked 
and hence not moving. This velocity fluctuation Is attributed to the Schuler 
error. The INS was left on during the entire time span covered In the plots. 
The magnitude of the errors are within expected limits. Both figures show one 
complete cycle of a near perfect 84-minute Schuler oscillation in the 
post-flight data In the latter part of the test period. This is in keeping 
with Huber and Bogers (1983) who noted that near the ground without 
accelerations Involved the Schuler oscillations will have an 84.4-minute 
period. But, In the first half of Flight 74 the errors are more random in 
their behavior and the oscillation Is an Irregular period. This complicates 
attempts to model or predict the error In advance. Lenschow (1972) suggests 
that post-flight data recorded with a stationary aircraft be used to back out 
the error. He proposed to simply trace back a recorded post-flight error 
oscillation with an 84-minute period constant amplitude sinusoidal curve. The 
present study shows, however, that both the period and the amplitude of the 
velocity error are altered substantially during flight and thus the Lenschow 
(1972) approach would not be successful here. This observation is in keeping 
with Huber and Bogers (1983) physical description of the Schuler effects. 
Additional investigation of Flight 74 Is needed to determine if the INS errors 
are accurately described in Figures 4.6 and 4.7. While the Inertial velocity 
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rror In north Inertial speed on Flight 74. 


measurement errors strongly influence the horizontal wind vector calculations, 
they generally have little effect on the gust velocity computations because 
the effect of the slow variations in velocity is greatly diminished or 
eliminated when the average velocity is removed. 

4.5 Flow Vane Errors 

Ringnes and Frost (1985) observed In analyzing the 8-57B data that 
constant differences existed between the angles of attack measured at the 
three different stations along the wing. The constant offset from the true 
value again has little Influence on the computed turbulence since the mean 
value is removed during the computation. The angle of attack terms have 
negligibly small effect on the computed values and therefore the inaccuracies 
cause no problems of the total horizontal wind vector. The cause of the angle 
of attack difference, however, were attributed to misalignment of the wing tip 
booms. 

The average sideslip angles were also found to be different from the 
expected value. All aircraft are designed directionally stable and will fly 
with zero average sideslip angle unless forcefully kept in a sideslip flight 
condition. The average sideslip angle of 2.23 degrees, for example, recorded 
at the centerboom on Flight 63 is therefore attributed to error. The source 
of the error is not clear but boom misalignment or problems with the data 
acquisition system are suspected causes. Again, the average sideslip error Is 
removed In the turbulence calculations, but it does affect the computed value 
of the horizontal wind vector noticeably as will be demonstrated. 

4.6 Airspeed Errors 

Frost, et al. (1985a) observed a difference in airspeed measured by the 
three separate wing probes. They compared average values for all runs on 
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Flight 21 and reported an average difference between the right and center boon 
measured velocltes of 1.82 m/s. The difference between the right and left 
wing tip measured airspeed was 0.79 m/s. The overall averaged airspeed was 
about 105 m/s. In Flight 31, the airspeed difference between the right and 
left probes Is 0.3 m/s at an average airspeed of 102 m/s. The accuracy of the 
horizontal wind vector calculations depends upon the quality of airspeed 
measurements. Possible Instrument calibration, position errors, or conversion 
from Indicated to true airspeed can cause these Inaccuracies. Also, the lack 
of separate static pressure transducers at the wind tips could have 
contributed to the Inconsl stances. A test flight conducted with the B-57B 
also revealed a value of horizontal wind speed of 2.5 m/s lower at the center 
boom than at the wing tip booms at a relative airspeed of roughly 122 m/s 
(Ehernberger, 1987). An approximate analysis based on a potential flow 
solution for a Ranklne body (Karamchetl, 1966) predicted a 6 percent error. 
This Is expected to be high because the B-573 is a more streamlined body than 
a Ranklne body, but the results do support the hypothesis that the airspeed 
may be retarded sufficiently by the aircraft body to produce the relative 
airspeed difference a boom's length from the nose. This 2 percent error Is 
accounted for In the following Investigation of the influence of 
Instrumentation errors on horizontal wind calculations. 

4.7 Gust Velocity Corrections 

The only instrumentation errors of those reported above which would 
noticeably effects the gust velocity calculation based on Inspection of 
Equations 4.5, 4.6, and 4.7 Is airspeed. The magnitude of correction to the 
gust velocities due to airspeed corrections Is illustrated in Figure 4.8. 
Uncorrected turbulence Is plotted directly from the tapes received from NASA 
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Langley Research Center. The "corrected" turbulence has been computed with 
the predicted Inertial velocity, airspeed, and .sideslip angle errors removed. 
The differences between the two computations are small and only detectable for 
the lateral and vertical components where total airspeed enters the Equations 
4.5 and 4.7. It Is apparent that even the airspeed error Is of little 
significance In gust velocity calculations. 

4.8 Horizontal Wind Vector Correction 

The INS velocity and position Indication, sideslip angle, and airspeed 
errors Identified as described above have been removed from the recorded data 
on some runs of Flight 63. The influence these errors have on the calculation 
of horizontal winds are demonstrated In this subsection. A series of wind 
vectors are plotted before and after corrections have been made along the 
flight path recorded by the INS. Each vector represents a one-second average 
from the 40 samples per second data tapes. 

In Figure 4.9 one of the box patterns flown on Flight 63 Is plotted. In 
this figure, no corrections have been made. There are some obvious 
Inconsistencies In the wind vectors, particularly, at the corners where It is 
expected that the wind would agrt.: closer between the two runs. The aircraft 
made 270-degree turns between runs which take less than two minutes. The wind 
direction is not expected to change significantly during that short of an 
Interval. Instrumentation errors are, therefore, the probable cause for the 
discontinuities In wind direction. Figure 4.10 differs from Figure 4.11 only 
by removal of the 2.23-degree sideslip error In the calculation of the wind 
vectors. It Is debatable whether this correction alone has Improved the wlru. 
vectors but It clearly demonstrates that seemingly small errors have 
significant effect on the wind vectors. In Figure 4.11 corrections have been 
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Fic„.e 4.11. Horizontal wind vectors on Flight 63 after airspeed, beta, 
and inertial velocity and position corrections (Runs 9 
through 12). 
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made for the known errors. The discontinuities In the wind vectors at the 
comers have all but vanished except for the bottom left-hand corner. 
However, as the numerical order of the runs indicates the box pattern was 
flown In a clockwise direction; thus, the beginning of Run 9 and the end of 
Run 12 are separated in time by approximately 15 minutes. Therefore, it is 
conceivable that the wind could have changed In that time span. The position 
errors are not severe for this box pattern but still noticeable. 

Figure '4.12 1's similar to Figure 4.9 except Runs 13 through 16 on Flight 
63 have been plotted. No corrections have been made. Only the disconti- 
nuities In the direction of the wind In the upper left-hand corner and in the 
magnitude of the wind In the upper right-hand corner appear questionable. 
Figure 4.13 illustrates the effect of removing the errors on wind vectors and 
INS Indicated locations. The horizontal wind vectors are more consistent and 
also the location of the runs are In better agreement with the flight plan. 

A third box pattern on Flight 63 (Runs 17 through 20) does not show the 
same Improvement with corrections. Figure 4.14 shows the uncorrected wind 
vector and Figure 4.15 the corrected version. The INS Indicated location Is 
Improved but not the wind vectors. After correction, the wind directions on 
Runs 18 and 20 are In sharp contrast to each other and additional or better 
corrections are needed. 

4.9 Effects of Non-Level Flight 

The algorithms used by NASA Langley computer facility to compute the 
turbulent gust velocities from the measured aircraft data are based on the 
assumption of straight, level flight. The more complete generalized system of 
equations which will allow for departure from level flight are derived In 
Appendix B. Questions arose during the study as to whether those portions of 
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flights for which the aircraft climbed or descended should be removed from the 
data. For example, during Run 3 at approximately 536 seconds (34 miles) Into 
the flight the aircraft climbed approximately 1000 ft (see Figure A.ll). In 
turn, during Flight 9 the aircraft climbed roughly 1000 ft beginning at t * 80 
seconds (7 miles) and descended again at t = 135 seconds (8.4 miles), see 
Figure A. 41. Also, Runs 1 and 2 where the aircraft took off or made 
touch-and-go's. This section shows that algorithms to reduce the data based 
on small angles and perturbations have no significant effect on the computed 
turbulence for runs where departures from straight and level flight occur. 

First It should be noted that because of the exaggerated vertical scale 
in, say for example. Figures A.ll, A.41, and all other plots of this nature 
given In Appendix A, the departure from level flight appears to be severe. It 
should be noted, however, that In no cases Is the climb or descent angle 
greater than 7*. This size angle adequately satisfies the small angle 
requirement defined in the algorithms presently used In the data reduction 
process. However, this statement Is further supported by quantitative 
analyses In the following. 

To Investigate the effects of climb and descent angles on the computed 
gust velocities. Equation B.27, which are used In the NASA Langley algorithms 
and Equation B.15 which FWG has derived and programmed to Investigate the 
effect of "large angles" where compared. The FWG equation still assumes the 
mean roll angles, <fr. Is zero. Equation B.27 was programmed and the turbulence 
time histories at central probe for the descending (e = -2.83*) and dinting 
(6 * 2.9*) segments of Run 2 were computed separately. Figure 4.16 shows the 
comparison of the descending segment. Figures 4.16a. 4.16b, and 4.16c are for 
the longitudinal (u), lateral (v), and vertical (w) components, respectively. 
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The standard deviation of each time history Is also shown In the figure. 
Similar comparison for the climbing segment is shown In Figure 4.17. An 
abnormal spike occurs In the lateral turbulence component at t * 201 seconds. 
It Is believed that the attitude of the aircraft at this moment probably 
deviated from the small angle assumptions significantly. Therefore, the 
maximum difference of the turbulence calculations from Equations B.15 and B.27 
occurs at this point. Although the complete equation (Equation B.l c ) 
calculates turbulence more accurately. Equation 3.27 saves a lot of computer 
time and still holds an acceptable accuracy for small angles considered here. 
For all practical purposes, the two calculations will Introduce only 
negligible difference In the turbulence analyses presented In this report. 

To further Investigate the departure from level flight, the turbulence 
time history for Flight 3 was divided Into two segments. Segment 1 Is from 0 
to 512 seconds (0 to 32 miles) and segment 2 from 512 to 691 seconds (32 to 44 
miles). The turbulence statistics were computed for the total run and for 
each segment Individually. The spectra for the two individual legs of the 
flight are compared with the total run In Figure 4.18; no apparent difference 
Is observed. The turbulence Intensity for each segment of the flight are 
listed In Table 4.1. Difference In turbulence Intensity for each leg of the 
flight are apparent. These differences, however, are not attributable to 
departure from level flight but rather due to patchiness of the turbulence 
associated with terrain features beneath the flight path. Figure 4.19 shows 
the turbulence time history and the approximate location relative to the 
underlying terrain at which the measurement was made. This figure. In view of 
the fact that the mean wind Is essentially out of the plane of the paper at 
approximately 15* , clearly suggests that strong turbulence Is associated with 
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TABLE 4 1. Turbulence Intensities for Total Time History and Each Segment 
Individually for Run 3. 


Comoonent 

Total 

Seament 1 

Seament 2 

ff UR 

3.17 

2.64 

2.87 

0 V r 

5.25 

5.29 

4.27 

*WR 

2.36 

2.10 

4.55 


flow over the mountain peaks. Thus, It Is concluded that the patchiness of 
the turbulence is due to terrain effects and not associated with any departure 
of the aircraft from straight and level flight. 
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5.0 CONCLUSIONS 
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The results of the analysis of Flight 31 coupled with experience from 
previous analysis of flight data from the NASA B-57 aircraft gust gradient 
program, lead to the following conclusions and recommendations: 

1. The probability density distribution of gust velocities In the 
atmosphere are not Gaussian. The distribution of velocity 
differences across the airfoil which filters 4>ut trends- In the 
quasi-steady wind have a definite modified Bessel function type 
distribution, l.e., a higher percentage of small and large 
velocity differences and lower percentage of Intermediate values 
than Is predicted by a Gaussian distribution. The parameter r of 
the modified Bessel function distribution however could not be 
related to the existing meteorology or to specific terrain 
features. It Is recommended that additional work to establish a 
physical meaning of the parameter r be carried out. The 
probability density distribution of the gust components 
themselves, l.e., not the difference, were rather ill behaved In 
this study and in many cases showed blmodal distributions. This 
Is believed to be due to the fact that a trend due to spatially 
varying mean wind along the flight path, caused by terrain 
features or other factors, were not removed from the gust 
velocities when computing the probability density functions. 

2. The theoretical von Karman spectrum fits the turbulence data well 
over the frequency range Investigated in this study (0.04 to 20 
Hz). The theoretical models were computed with length scales 
determined from Integration of the correlation coefficient from 

zero lag to the point where the correlation first becomes zero. 
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The results of the study strongly suggest that the turbulence 
behaves relatively consistent with the assumption of Isotropic, 
homogeneous turbulence despite the fact that flights were made 
over mountainous terrain and during touch-and-go's through the 
atmospheric boundary layer. 

3. The two-point common component theoretical spectra (that Is, the 
spectra for the same velocity component for spatially separated 
positions) proposed by Houbolt and Sen (1972) for the vertical 
fluctuations and the spectra for longitudinal and lateral 
velocity fluctuations derived In this report, agree with the 
experimental flight data provided the care described In Chapter 3 
Is exercised In computing the two-point spectra from the 
truncated, digitized gust velocity time histories. 

4. In calculating one-point auto-spectra and two-point conrnon 
component spectra from a direct Fourier transform of the data can 
result In large errors due to aliasing and truncation estimator 
bias. Aliasing Is the major source of error In the auto-spectrum 
whereas bias Is the major source of error In the two-point common 

component spectra. This Is physically evident since the energy 

* % 

contained at high frequencies In two-point spectra vanishes 
significantly faster than that In one-point spectra. What Is not 
evident Is that a very small departure of two-point correlation 
^coefficient from unity at zero lag can cause high bias errors. 
To remove bias error from two-point spectra It Is recommended 
that the minimum-bias lag window be utilized. 

5. To reduce the variance error associated with a two-point common 
component as well as all cross-spectra it Is Important to carry 
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out segment averaging of the co- and quad-spectra C(s,f) and 
Q(s*f) separately as contrasted to averaging the absolute value, 

l.e., k2(s,f) + Q2(s,f). 

6. The small values (0.07) of the ratio of the spanwldth separation 
distance of 20 m to the typical turbulence Integral length scale 
was found to have a relatively significant effect on the 
two-point spectra In terms of spectrum ..dropoff, at high 
frequencies. However, for other statistical parameters; 1e, 
cross correlations, there were Indications that the separation 
distance of 20 m Is too small to resolve some of the statistical 
Issues of Interest. It Is recommended that experiments be 
carried out with larger separation distances than 20 m for a 

firmer understanding of two-point statistical pamameters. 

7. The one-point and two-point cross-correlations between uncooroon 
velocity components shows, almost zero correlation (further 
supporting the assumption of Isotropic turbulence). However, 
close examination of the complex phase angle associated with the 
cross-spectra Indicate there Is significant phase difference 
between various frequencies. The one-point cross-spectrum appears 
to agree well with the model proposed by Reeves, et al. (1974). 
No empirical expression or analytical model of two-point cross- 
spectra Is available. Further work Is required In this area. 

8. The Instrumentation system and data processing algorithms for the 

NASA B-57B aircraft provides highly accurate measurements of 

turbulent gust velocities. The measurements of the total 

Instantaneous wind speed, however, may contain certain errors 

Induced by some of the present characteristics of the measuring 
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system. In particular the INS Schuler drift problem causes 
significant uncertainty In the position of the flight path and of 
the magnitude of the mean wind speed (l.e., ±2 to 5 m/s). This 
uncertainty In velocity coupled with small Inaccuracies In the 
flow vane measurements (possibly due to boom misalignment or 
other factors), while having Insignificant effect on the gust 
velocity measurement, can result In major errors In the wind 
field. These errors can be corrected If appropriate data other 
than that measured by the on-board Instrumentation system Is 
gathered during the flight. For example, visual observed 
position recorded with a designation marker utilized can be used 
to estimate INS Schuler position drift from post-analysis of the 
data. Since the Schuler drift does not appear to have a constant 
amplitude nor period of oscillation, procedures to correct for 
this error by backing cut the Inertial winds and position from 
measurements made with the aircraft stationary on the runways are 
not feasible. 
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APPENDICES 


APPENDIX A 


RESULTS OF STATISTICAL ANALYSIS OF FLIGHT 31 

General Information, ground track terrain features, and statistical values 
for all runs (except landing operation. Run 17) of Flight 31 on November 29, 
1982, are presented in this appendix. The analysis of each run is given in 
two tables and five figures. The first table shows the turbulence average 
parameters, integral length scales, and correlation coefficients and the 
second one lists all parameters measured and their range of values. Five 
figures show the flight altitude, time history, probability density function, 
normalized correlation function, and normalized spectral density function of 
gust velocities, respectively. 
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Figure A. 2. (continued) 



TABLE A. 1. Average Turbulence Parameters, Integral Length Scales, and 

Correlation Coefficients of Gust Velocities, Flight 31, Runl. 


Mean Airspeed (m/s): 

4. Integral Length Scale (m) 
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l vR L w r 
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Gust Velocity Difference Distributions 


Figure A. 3. Probability density functions for gust velocities and 
gust velocity differences (normalized with the standard 
deviation). Flight 31, Run 1 (r = degree of non-Gaussian) 
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a. One and two-point common component spectra. 

Figure A. 5. Comparison of normalized one- and two-point spectral density functions 
for gust velocities with theoretical models, Flight 31, Run 1. 
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One- and two-point cross-spectra. 
Figure A. 5. (continued). 
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Figure A. 7. Time histories of gust velocities, gust velocity differences, and aircraft 
normal accelerations. Flight 31, Run 2. 
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Figure A. 7. (continued) 



TABLE A. 3. Average Turbulence Parameters, Integral Length Scales, and 

Correlation Coefficients of Gust Velocities, Flight 31, Run 2. 


1. Mean Airspeed (m/s): 4. 
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Figure A. 8. Probability density functions for gust velocities and 
gust velocity differences (normalized with the standard 
deviation), Flight 31, Run 2 (r = degree of non-Gaussian) 
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Figure A. 9. Comparison of normalized one- and two-point correlation functions for 
gust velocities with theoretical models. Flight 31, Run 2. 
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a. One- and two-point common component spectra. 

Figure A. 10. Comparison of normalized one- and two-point spectral density functions 
for gust velocities with theoretical models. Flight 31, Run 2. 



One- and two-point cruss-spectra. 
Figure A. 10. (continued). 
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TABLE A. 4. 


List of All Parameters Measured and Their Range of Values, 
FI i ght 31 , Run 2. 
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Flight 31 , Run 3 
Date: Nov . 29, 1982 

Start Time: 10:44:29 (PST) 

Duration: 694.8 seconds 



(TSV Ti) TM6i3H 


Figure A. 11. Flight altitude and horizontal wind along flight path. Flight 31, Run 
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normal accelerations, Flight 31, Run 
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Figure A. 12. (continued). 
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Time (seconds) 
Figure A. 12. (continued). 
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Time (seconds) 
Figure A. 12. (continued). 



TABLE A. 5. Average Turbulence Parameters, Integral Length Scales, and 

Correlation Coefficients of Gust Velocities, Flight 31, Run 3. 


1. Mean Airspeed (m/s): 

VL Vc Vr 

104.21 102.52 104.60 

2. Standard Deviation of 
Gust Velocities (m/s):. 


a uR 

a vR 

CT wR 

3.17 

5.25 

2.36 

a uC 

°vC 

a wC 

3.15 

5.29 

2.18 

CT uL 

a vL 

ff wL 

3.19 

5.31 

2.31 


3. Standard Deviation of Gust 
Velocity Differences (m/s): 

a AuRL ct avRL ff AwRL 
1.29 1.23 1.37 


Integral 

Length Scale (m’ 

L uR 

l vR 

l wR 

234.0 

425.6 

116.9 

l uRL 

l vRL 

l wRL 

238.4 

422.8 

115.3 


5. Correlation Coefficient of Gust 
Velocities: 
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TTO/^UR^VR 

^R w r/ (J vr ct wr 

w R u r/^wr°ur 

0.09 

-0.19 

0.06 

TJRVL/a UR CT VL 

VRWL/a VR a WL 

WrOl/c Wr c Ul 

0.04 

-0.19 

0.05 
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Gaussian 

Non-Gaussian 



Gust Velocity Distributions 
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Gust Velocity Difference Distributions 


Figure A. 13. Probability density functions for gust velocities and 
gust velocity differences (normalized with the standard 
deviation), Flight 31, Run 3 (r = degree of non-Gaussian) 











a. One- and two-point common component spectra. 

Figure A. 15. Comparison of normalized one- and two-point spectral density functions 
for gust velocities with theoretical models. Flight 31, Run 3. 
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One- and two-point cross-spectra. 
Figure A. 15. (continued). 



ORIGINAL PACE IS 
OF POOR QUALITY 


TABLE A. 6. 


List of All Parameters Measured and Their Range of Values, 
FI i ght 31 , Run 3. 
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(TABLE 

A. 6. continued) 
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normal accelerations. Flight 31, Run 




Figure A. 17. (continued). 



(m/s) 



Time (seconds) 
Figure A. 17. (continued). 



TABLE A. 7. Average Turbulence Parameters, Integral Length Scales, and 

Correlation Coefficients of Gust Velocities, Flight 31, Run 4. 


1. Mean Airspeed (m/s): 

V L V C Vr 

104.78 102.62 104.32 

2. Standard Deviation of 
Gust Velocities (m/s): 


*uR 

°vR 

a wR 

a uC 

CT vC 

°wC 

a uL 

*vL 

c wL 


3. Standard Deviation of Gust 
Velocity Differences (m/s): 

ct auRL ct avRL ct awRL 


4. Integral Length Scale (m): 


*-uR lyR 1-wR 

419.8 350.8 66.9 

l uRL l vRL l wRL 

408.0 344.7 61.9 


5. Correlation Coefficient of Gust 


Velocities: 

u R u L/°u R a u L 

^7RVL/av R CTv L 

WrWl/ct W rCT W |_ 

0.88 

0.91 

0.80 

u R v r/ c ur 0 v r 

v R w r/ <j vr q w R 

W R U R/°WR a UR 

-0.19 

0.20 

0.09 

URVL/a UR av L 

v RWl/°vr^w l 

WRUL/a WR a UL 

-0.19 

0.20 

0.06 



Probability Probability 


Gaussian 

iion-Gaussian 


r = 2.5 r = 1 .0 r = 2.5 



Gust Velocity Distributions 



Gust Velocity Difference Distributions 


Figure A. 18. Probability density functions for gust velocities and 
gust velocity differences (normalized with the standard 
deviation). Flight 31, Run 4 (r = degree of non-Gaussian) . 
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One- and two-point cruss-spectra. 
Figure A. 20- (continued). 


ORIGINAL PAG* B 
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TABLE A. 8. List of All Parameters Measured and Their Range of Values, 
FI i ght 31 , Run 4. 
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Date: Nov. 29, 1932 

Start Time: 11:00:55 (PST) 

Duration: 145.0 seconds 



(1SV Ti) Titian 
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Figure A. 21. Flight altitude and horizontal wind along flight path. Flight 31, Run 
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normal accelerations. Flight 31, Run 






Time (seconds) 


Figure A. 22. (continued). 
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TABLE A. 9. Average Turbulence Parameters, Integral Length Scales, and 

Correlation Coefficients of Gust Velocities, Flight 31, Run 5. 


1. Mean Airspeed (m/s): 

\ Vc Vr 

105.79 103.53 105.33 

2. Standard Deviation of 
Gust Velocities (m/s): 


CT uR 

CT vR 

a wR 

2.49 

4.06 

2.76 

a uC 

ff vC 

°wC 

2.47 

4.04 

2.66 

CT uL 

ff vL 

17 wL 

2.56 

4. 10 

2.85 


3. Standard Deviation of Gust 
Velocity Differences (m/s): 

°auRL p avRL ^awRL 
1.41 1.38 1.42 


4. Integral Length Scale (m): 


LuR lyR 

333.9 168.6 189.7 

luRL ^VRL LwRL 

317.5 173.5 204.0 


Correlation 

Velocities: 

Coefficient of Gust 

u R u L/ a UR CT UL 

VrVl/ctv r ctv l 

WR W L/ a WR a W|_ 

0.87 

0.90 

0.90 

U R V R/°UR a VR 

VrWr/OvrCTwr 

WRUR/a WR a UR 

-0.09 

-0.10 

-0.17 

URVL/a UR cJv L 

VRWL/a VR CT WL 

w R u L/°WR a ui_ 

-0.08 

-0.09 

-0.20 
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Probability Probability 


Gaussian 

Non-Gaussian 


r = 3.5 r = 3.5 r = 3.5 



Gust Velocity Distributions 


r = 4.5 r = 7.5 r = 6 . 5 



Gust Velocity Difference Distributions 


Figure A. 23. Probability density functions for gust velocities and 

gust velocity differences (normalized with the standard 
deviation). Flight 31, Run 5 (r = degree of non-Gaussian) . 
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Aircraft data 
von Karma n 



a. One- and two-point common component spectra. 

Figure A. 25. Comparison of normalized one- and two-point spectral density functions 
for gust velocities with theoretical models, Flight 31, Run 5. 
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One- and two-point cross-spectra. 
Figure A. 25. (continued). 
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A. 10. List of All Parameters Measured and Their Range of Values, 
FI i ght 31 , Run 5. 


CHANNEL 

UNITS 

HIGH 

LOW 

MEAN 

RMS 

STO 

POINTS 

1 TIME 

5EC0NDS 

40279.553 

40134.576 

40207.06540 

40207.08719 

41.86150 

5 80 D 

? PHI DOT 

R AD/SEC 

. 20 8 

-.195 

-.00218 

.04443 

.04430 

5600 

3 AC C L N CG 

G UNITS 

2.327 

.230 

.99249 

1.00515 

.15907 

5800 

< THETA DQT 

RAD/SEC 

.153 

-.151 

.00264 

.01892 

.01873 

5800 

5 THETA 

RAD 

.154 

-.033 

.04342 

.05090 

.02655 

5800 

6 PHI 

RAD 

.133 

-.184 

.00088 

•04024 

.04023 

5800 

7 PSI 1 

RAD 

359. 110 

-.155 

299.64552 

326.65102 

130.06294 

5000 

8 DEL PSI 1 

DEG 

7.929 

-9.556 

-.29595 

2.35956 

2.34113 

5800 

9 PS! 2 

RAD 

363.571 

346.672 

355,57109 

355.57842 

2.28427 

5800 

10 DEL PSI 2 

DEG 

2.13B 

-15.092 

-2.84581 

6.26277 

5.57933 

5800 

11 accl n lt 

G UNITS 

4.691 

-2.620 

1.00522 

1.06613 

.35524 

5 bOO 

12 ACCL N RT 

G UNITS 

4.993 

-2.651 

1.00395 

1.07230 

.37674 

5600 

13 ACCL X CG 

G UNITS 

.202 

-.013 

.04563 

.05809 

.03595 

5B00 

14 ACCL Y CG 

G UNITS 

.143 

-.109 

-.00165 

•02553 

.02548 

5600 

15 ALPHA CTR 

RAD 

.126 

-.106 

-.02076 

.02707 

.01739 

5600 

16 BETA CTR 

RAD 

.178 

-.166 

.00413 

.03557 

.03533 

5800 

I 7 TEMP I 

DEG F 

76.497 

74.076 

75.76997 

75.77070 

,35062 

5800 

18 TEMP P 

DEG F 

61.781 

61.425 

61.61107 

61.61108 

.03676 

5600 

ACCL 2 INS 

G UNITS 

2.272 

.202 

.99627 

1.00934 

.16193 

58 00 

*0 ALPHA RT 

RAD 

.140 

-.105 

-.01016 

.02160 

,01906 

5800 

21 BETA HT 

RAD 

.190 

-.165 

.01321 

.03466 

.03227 

5800 

22 ALPHA IT 

RAD 

.133 

-.086 

.C1B99 

.02656 

.01857 

5800 

23 BETA LT 

RAD 

.142 

-.161 

.00764 

.03341 

.03253 

5600 

2* PSI ODT 

RAD/SEC 

.103 

-.204 

.00271 

.03135 

.03123 

5000 

25 TEHP TOT 

DEG C 

14.259 

8.057 

11.50946 

11.58680 

1,33665 

5800 

26 QC LT 

PSID 

1.121 

.604 

.02579 

.82801 

.06055 

5600 

27 oc CTR 

PSID 

1.003 

.578 

.70997 

.79214 

•05061 

5600 

25 oc RT 

PSID 

1.023 

.560 

•61040 

.82065 

• 0b070 

5600 

29 PS" 

PSIA 

11.581 

11.376 

11.521 20 

11.52138 

.04622 

5800 

70 )TEHP IRT 

VOLTS 

7.371 

2.638 

.6.51844 

6.54426 

.58078 

5800 

31 HVGROM 

DEG C 

-1.547 

-9.759 

-5.86794 

6.15471 

1.65697 

5800 

3? QC2 LT 

PSID 

.070 

.066 

.06965 

•06965 

.00042 

5800 

33 0C2 CTR 

PSID 

.177 

.167 

.17330 

.17333 

•00303 

5800 

34 QC2 RT 

PSID 

.163 

-.035 

,15347 

.13364 

.00727 

5600 

35 PAR 

DEG 

1.169 

-8.916 

-0.76417 

8.76534 

.14316 

5800 

36 DAL 

DEG 

1.513 

-8.596 

-8.47692 

8 .47811 

.14203 

5600 

37 DELEV 

DEG 

3.393 

-15.431 

3.27645 

3.28624 

,25345 

5800 

33 0 STAB 

DEG ' 

.764 

-.26? 

-.25796 

.25836 

,01427 

3800 

39 DRUG 

DEG 

7.612 

-14.981 

7.36917 

7.37639 

.32644 

5800 

40 0 THR R ' ” 

PCT MAX 

61.320 

23.376 

60.04641 

60.84900 

.56216 

5600 

41 DTHRL 

PCT H AX 

60.938 

23.288 

60.44540 

60.44005 

.56558 

5600 

42 DFLP '* " 

POSITION 

1.960 

.559 

.57000 

.57922 

•02204 

5800 

43 D SB 

POSITION 

.174 

.058 

•16219 

.16234 

.00704 

5000 

<4 0 TO 6 

HETER5 

7531608.000 

7525273.037 



1070.13318 

3800 

45, B TO 0 

DEGREES 

72.036 

72.799 

72.61743 

72.61743 

•01070 

3800 

46. LONG 

DEGREES 

-117.953 

-117.976 

-117.96500 

117.96500 

•00601 

3000 

47 LAT 

DEGREES 

35.763 

35.629 

35.69505 

35.69567 

•03898 

5800 

4 ft TRK AN 6 

DEGREES 

10.653 

5.074 

7.73019 

7.07254 

1.44833 

5800 

40 HOG 

RADIANS 

.112 

-.204 

-.03712 

•05631 

.04235 

5600 

60 VE 

H/SEC 

19.460 

8,599 

13.92566 

14.19503 

2.75239 

3000 

51 VN - 

H/SEC 

105. 400 

91.695 

102.33217 

102.37126 

2.62096 

5800 

r ? altitude 

KH 

2.108 

1.964 

Z.006C0 

2.00626 

.03237 

5 800 

t\ TENPC 

DEGREES C 

6,949 

3.293 

6.16740 

6.27054 

1.13273 

5800 

54 F* WND SPD 

KNOTS 

50.684 

5.746 

35.36690 

36.22351 

7.83171 

5800 

55 NS WHO SPD 

KNOTS 

14,845 

-29.404 

-2.01503 

3.26163 

4.86080 

5600 

55 WIND SPEEO 

KNOTS 

51.566 

7.565 

35.81271 

36.60369 

7*56900 

5800 

57 WIND DIREC 

DEGREES 

334.159 

216.077 

273.71739 

273.07670 

9.34098 

5800 

5 B WIND DIAZ 

DEGREES 

154.159 

38.077 

93,71744 

94.10172 

9.34099 

5000 

59 WINO 0 1 R 3 

DEGREES 

334.159 

218.077 

273.71744 

273.87675 

9.34099 

5000 

Rn WlN0 DIR4*“ 

DEGREES 

334.159 

216.077 

273.71744 

273.87675 

9. 34099 

5000 

rj AIRSPEED R 

H/5EC 

117.187 

89.626 

105.33370 

105.40679 

3.92291 

5800 

jjl AIRSPEED C 

H/SEC 

115.996 

89.111 

103.53162 

103.60331 

3.84857 

5600 

AIRSPEED L * 

H/SEC 

122.810 

90.975 

105.79945 

105. 07G52 

3.07892 

5800 

DELTA ALT 

METERS 

59,411 

-04.477 

-42.43732 

53.37152 

32.36932 

5800 

INRTL OISP 

METERS 

37.737 

-82.204 

-41.49497 

52.29241 

31.02511 

5800 

66 UG RIGHT 

M/SEC 

9.156 

- 11.101 

.00000 

2.49780 

2.49601 

5600 

67 UG CENTER 

M/SEC 

fl.546 

-14.061 

.00000 

2.47382 

2.47403 

5800 

59 UG LEFT 

M/SEC 

8,044 

-18.166 

.00000 

2.56170 

2.56200 

5 800 

6 g VG RIGHT 

M/SEC 

6.832 

-17.989 

-.01423 

4.06130 

4.06163 

5800 

70V& CENTER 

H/SEC 

8,216 

-15.520 

-.01413 

4.04630 

4.04662 

560 0 

71 VG LEFT 

M/SEC 

7.272 

-13.953 

-.01740 

4.10711 

4.10730 

5600 

77 W 6 RIGHT ” 

M/SEC 

13.641 

-10.461 

.01471 

2.76231 

2.76238 

5600 

7 1 WG CENTER 

H/SEC 

13.351 

-9.627 

•02009 

2.66070 

2.66093 

5800 

74 W 6 LEFT 

H/SEC 

13.921 

-9.047 

.02031 

2.85220 

2.05245 

5 600 


20S 






(m/s) 
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normal accelerations. Flight 31, Run 



TABLE A. 11. Average Turbulence Parameters, Integral Length Scales, and 

Correlation Coefficients of Gust Velocities, Flight 31, Run 6. 


Mean Airspeed (m/s): 

4. Integral Length Scale (m) 

: 

vl 

v C 

v R 

1-uR L v r L w r 


104.31 

102.19 

104.01 

364.7 92.0 51.7 





LllRI LwRI L w ri 


Standard Deviation of 




Gust Velocities (m/s): 

344.5 104.2 47.5 


*uR 

G vR 

a wR 




3.64 

3.67 

3.41 

5. Correlation 

Coefficient of Gust 




Velocities: 



ff uC 

a vC 

G wC 




3.54 

3.65 

3.00 

WOL/°u R °u L 

VRVL/a VR a VL 

w R w L/ a w R a W[_ 




0.82 

0.90 

0.81 

a uL 

a vL 

° wL 




3.54 

3.42 

3.12 

u R v R/ a u R a v R 

W G VR G WR 

WRTJR/CT WR a U R 




-0.18 

0.60 

-0.10 

Standard Deviation of Gust 

u R v L/°UR a VL 

VRWL/a VR a WL 

w R u l/ g wr°ul 

Velocity Differences (m/s): 

-0.18 

0.61 

0.00 

*auRL 

a AvRL 

a AwRL 




1.74 

1.68 

1.92 
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One- and two-point cross-spectra. 
Figure A. 30. (continued). 
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TABLE A. 12. 


List of All Parameters Measured and Their Range of Values, 
FI ight 31 , Run 6. 


CHANNEL 

UNITS 

I.TJHE 

5ECDNDS 

* PHI DOT 

RAD/SEC 

; ACC l N CG 

G UNITS 

J THETA DOT 

RAD/SEC 

5 THETA 

RAD 

5 PHI 

RAD 

7 PSI 1 

RAD 

B DEL PSI 1 

DEG 

9 PSI 2 

RAD 

™ DU PSI 2 
T1 ACCL N LT 

DEG 

G UNITS 

12 ACCl N RT 

G UNITS 

13 ACCl X CG 

G UNITS 

14 ACCL V CG 

G UNITS 

15 ALPHA CTA 

RAD 

16 BETA CTR 

RAO 

IT TEMP I 

DEG F 

18 TEMP P 

DEG F 

19 ACCL l INS „ 

G UNITS 

70 ALPHA RT 

RAD 

21 BETA RT 

RAD 

22 ALPHA LT 

RAD 

23 BETA LT 

RAD 

24 psi DOT 

RAD/SEC 

26 TEMP TOT 

DEG C 

26 QC LT 

PS 1 D 

2? QC CTR 

PSID 

28 QC RT 

PSID 

29 ps 

PSIA 

3Q T6HP IRT 

VOLTS 

31 MYGROH 

DEG C 

32 QC 2 IT 

PSID 

33 QC 2 CTR 

PSIO 

34 0C2 RT 

PSID 

38"DAR 

DEG 

35 DAL 

DEG 

37 DELEV 

DEG 

38" 0 ST AG 

OEG " 

39 DRUD ”■ 

DEG 

4^DTHRR 

PCT MAX 

4 1 0THRL 

PC T MAX ** 

42 DFLP 

POSITION 

43 DSB ” ' 

POSITION 

44 D TO’G 

METERS 

45ft TO D 

DEGREES 

46 LONG ' 

DEGREES 

47 l AT 

DEGREES 

48 TRK AMG 

DEGREES 

49 HOG 

radians 

50 VE 

H/SEC 

51 VN 

H/SEC 

52 ALTITUDE 

KM 

53 TEMPC 

DEGREES C 

54 EW W NO SPD 

KNOTS 

55 NS WHO SPD 

KNOTS 

56 WIND SPEED 

KNOTS 

57 WIND OIREC 

DEGREES 

56 WIND DIR2 

DEGREES 

59 WIND 0 1 R3 _ 

DEGREES 

60 W I ND DIR4~’ 

DEGREES 

61 AIRSPEED R 

M/SEC 

62 AIRSPEED C 

M/SEC 

63 A I R 5 P E E 0 L 

H/SEC 

64 DELTA ALT 

METERS 

65.1 NR T L DISP_ 

METERS 

66UG _ RIGHT 

H/SEC 

67 UG CENTER 

H/SEC 

68 UG LEFT __ 

H/SEC 

69 VG RIGHT 

" M/SEC 

70 VG CENTER 

H/SEC 

71 VG LEFT 

H/SEC 

72 WG RIGHT 

H/SEC 

73 WG CENTER 

H/SEC 

74 WG LEFT 

H/SEC 


HIGH 

LOW 

MEAN 

MS 

STD 

POINTS 

40514.445 

40451.620 

40403.03270 

40403.03677 

16.14684 

2514 

.4 56 

-.235 

-.00310 

.06217 

.06210 

2314 

2.781 

-.523 

1.01158 

1.04299 

.25412 

2514 

.104 

-.196 

.00335 

•02706 

.02608 

2514 

.146 

,020 

.06762 

.09017 

•02131 

2514 

.263 

-.107 

-.00172 

.05970 

.05969 

2514 

196.300 

179.129 

184.51265 

184.53591 

2.93071 

2514 

14.903 

-2.562 

2.60766 

4.09430 

2.98058 

2514 

550.604 

537.401 

542.62956 

542.63938 

3.26409 

2514 

14.841 

-2.213 

3.10556 

4.24230 

2.8906b 

2514 

5.024 

-2.019 

1.02459 

1.16053 

.54311 

2514 

4.993 

-2.960 

1.01898 

1.17018 

.59156 

2514 

.233 

-.020 

. Q83 53 

.09360 

.04240 

2514 

.135 

-.164 

-.00067 

*03531 

.03530 

2514 

.116 

-.160 

-.01681 

•03140 

.02663 

2514 

,208 

-.165 

.00400 

.04619 

.04803 

2514 

78.476 

77.397 

77.96893 

77.96927 

.22974 

2514 

61.781 

61.601 

61.74117 

61.74122 

.07356 

2514 

3.045 

-.562 

1.01599 

1.04685 

.25236 

2514 

.150 

-.194 

-.00244 

•03131 

.03122 

2314 

.203 

-.167 

.01407 

.04765 

.04553 

2314 

.155 

-.151 

•02282 

.03562 

.02761 

2314 

.174 

-.164 

.00548 

.04307 

.04273 

2514 

• 1 59 

-.126 

.00190 

•04015 

.04011 

2514 

10.715 

7.072 

8,67362 

8.70365 

.72253 

2514 

l.UO 

.641 

.79533 

•80003 

.06656 

2514 

1.037 

.565 

.76269 

.76735 

.06445 

2314 

1.061 

.618 

.79085 

.79568 

.00761 

2514 

11.324 

11.237 

U. 29206 

11.29209 

.02389 

2514 

5.503 

3.121 

4.61436 

4.83667 

.46390 

2514 

-1.942 

-8.194 

-5.40948 

5.63043 

1.56212 

2514 

.072 

.072 

.07180 

• 071 BO 

.00025 

2514 

.150 

,137 

.14859 

.14070 

.00586 

2514 

.137 

.126 

.13074 

.13078 

.00301 

2514 

-9.689 

-9.717 

-9.66945 

9.6694b 

.00261 

2314 

-0.899 

-9.037 

-8.98630 

8.96639 

. 04009 

2514 

4.058 

3.785 

3.92190 

3.92248 

,06743 

2514 

-.283 

-.294 

-.26661 

.26663 

“ .00328 

2514 

6.295 

8.066 

8.15400 

6.15421 

.05977 

2514 

62.590 

62.207 

62.38175 

62.38182 

.09379 

2514 

62.305 

62.012 

62.16127 

62.16136 ' 

* .10322 

2514 

.512 

.496 

.50417 

•S041B 

.00443 

2514 

.207 

.197 

.20250 

•20252 

•00251 

2314 

752505B *456 

7523312.563 



324.55013 

2514 

72.024 

72.790 

72.81079 

72.81079 

.00720 

2514 

-117.997 

-118.000 

-117.99656 

117.99856 

.00061 

2314 

35.686 

, 35,625 

35.6)489 

35,65489 

.01733 

2514 

102.249 

174.899 

177.83402 

177,84529 

2.00259 

2314 

3.447 

3.138, 

3.23368 

3,23411 

.03251 

2514 

9.799 

-4.430 

4.09913 

5,61498 

3.03809 

2514 

-102.169 

-112.961 

-106.9U144 

106.96891 

3.79957 

2514 

2.206 

2.145 

2.16746 

2.16753 

.01690 

2514 

4.602 

2.085 

3.46107 

3.51936 

.63002 

2514 

42.305 

-1,321 

25.54551 

26.42022 

6.74336 

2514 

11.402 

-25.054 

-10,14093 

12.26632 

6.90396 

2514 

42.407 

3.865 

28.29386 

29,12970 

6.92939 

2514 

357.790 

6.301 

290.70819 

291.23767 

17.55703 

2514 

177.798 

-173. $99 

110.70625 

112.09123 

17.55705 

2514 

357.798 

6.301 

290,70025 

291.23773 

17.55705 

2514 

360.755 

101.166 

290.99465 

291.42208 

15.70111 

2314 

121.266 

92.271 

104.01912 

104.16727 

5.55470 

2314 

110.709 

09.929 

102.19461 

102.33983 

5.45103 

2314 

122.780 

93.967 

104.31353 

104.45593 

5.45346 

2314 

1.015 

-60.600 

-37.01934 

41,45465 

16.97942 

2514 

.000 

-52.813 

-30.60195 

34.32676 

15.55452 

2514 

8.054 

-11.346 

.00000 

3,64237 

3.64309 . 

2514 

8.341 

-10.706 

.00000 

3.54397 

3.54467 

2514 

7.622 

-13.375 

.00000 

3,54621 

3.54692 

2514 

17.341 

-9.962 

-.00764 

3.67656 

3.67624 

2514 

15.198 

-9.730 

-.00299 

3,63269 

3.65247 

2514 

16.523 

-9.308 

-.10410 

3,42577 

3.4248T 

2314 

17.793 

-11.422 

-.10907 

3.42038 

3.41952 

2514 

15.330 

-9.243 

-.09051 

3.00378 

3.00277 

2514 

20.597 

-12.616 

-.09179 

3.12148 

3.12075 

2514 
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Time (seconds) 
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Figure A. 32. (continued) 


TABLE A. 13. Average Turbulence Parameters, Integral Length Scales, and 

Correlation Coefficients of Gust Velocities, Flight 31, Run 7. 


1. Mean Airspeed (m/s): 

Vl Vc Vr 

101.47 99.23 100.93 

2. Standard Deviation of 
Gust Velocities (m/s): 


ff uR 

a vR 

°wR 

3.03 

3.00 

2.23 

°uC 

p vC 

ff wC 

3.03 

3.03 

2.15 

a uL 

°vL 

ff wL 

3.07 

3.06 

2.22 


3. Standard Deviation of Gust 
Velocity Differences (m/s): 

ff AuRL °avRL ct awRL 
0.85 0.80 0.89 


4. Integral Length Scale (m): 


LuR *-vR ^-wR 

562.8 249.6 287.6 

LuRL l vRL l wRL 

532.2 242.9 283.5 


5. Correlation Coefficient of Gust 
Velocities: 


u R u L/ a UR a U|_ 

VRVL/a VR a VL 

w R w L/ a WR a W[_ 

0.92 

0.90 

0.90 

urvr/^ur^vr 

v R w R/ 0 vr^wr 

WRUR/a WR o UR 

0.02 

-0.21 

0.39 

TJRVL/a UR a VL 

V R W L/ a VR a WL 

WRUL/a WR au L 

0.01 

-0.21 

0.32 
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One- and two-point cruss-spectra. 
Figure A. 35. (continued). 
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TABLE A. 14. 


List of All Parameters Measured and Their Range of Values, 
FI i ght 31 , Run 7. 


CHANNEL 

UNITS 

I TIME 

SECONDS 

« PHI DOT 

RAO/SEC 

3 ACC L M CG 

G UNITS 

4 THETA DOT 

RAD/SEC 

5 THETA 

RAD 

6 PHI 

RAD 

7 PSI 1 

RAD 

8 DEL P SI 1 

DEG 

9 PSI 2 

RAD 

10 del PSI 2 

DEG 

11 ACCL N LT 

G UNITS 

12 ACCl N RT 

G UNITS 

13 ACCl X CG 

G UNITS 

H ACCL Y CG 

G UNITS 

15 ALPHA CTR 

RAD 

16 BETA CTR 

RAD 

] l TEMP I 

DEG F 

16 TEHP p 

DEG F 

19 ACCL 2 IMS 

G UNITS 

20 ALPHA RT 

RAD 

21 BETA RT 

RAD 

22 ALPHA LT 

RAD 

23 BETA LT 

RAD 

n PSI DOT 

RAD/SEC 

25 TEHP TOT 

DEG C 

26 OC IT 

PSIO 

27 OC CTR 

PSIO 

26 CC RT 

PSID 

29 PS" 

PSIA 

30 TEHP IRT 

VOLTS 

31 HYGROM 

DFC C 

32 0C2 LT 

PSID 

33 0C2 CTR 

PSID 

34 GC2 RT 

PSID 

35 dar 

DEG 

36 dal 

DEG 

37 DELEV 

DEG 

38 f) $ T A 6 

DCG 

39 DRUD 

dtc 

40 n THRR 

PC T HA X 

41 D 1 HR L 

PCT MAX 

42 DHP 

POSITION 

43 DSB 

POSITION 

44 o TO G 

HETERS 

4^ B TO D 

DEGREES 

46 EON G 

DEGREES 

47 LAT 

DEGREES 

48 TP K AMG 

DEGREES 

49 HUG 

RADIANS 

bO VE 

H/SEC 

51 VM 

H/SEC 

52 ALTITUDE 

KH 

53 TEHPC 

DEGREES C 

54 E V VND SPD 

KNOTS 

55 NS WHO SPD 

KNOTS 

56 VINO SPEED 

KNOTS _ 

57 V IN D DIREC 

DEGREES 

58 MIND DIR2 

degrees 

59 WIND DIR3 _ 

DEGREES 

60 WIND DIR4 

DEGREES 

61 AIR SP EE 0 R 

H/5EC 

62 AIRSPEED C 

H/SEC 

63 AIRSPEED l 

* H/SEC 

64 DELTA ALT 

HETERS 

65 INRTL D I SP 

HETERS 

66 UG RIGHT 

h/SEC 

67 UG CENTER 

H/5FC 

6 B UG LEFT 

H/SEC 

69 VG RIGHT 

PI/S EC 

70 VG CENTER 

H/SEC 

71 VG LEFT 

H/SEC 

72 WG RIGHT 

H/SEC 

73 WG CENTER 

H/SEC 

74 WG LEFT 

h/SEC 


HIGH LOW MEAN RMS STD POINTS 


40859. 30Z 

40655.627 * 

.160 

-.143 

2.067 

.159 

.099 

-.110 

.119 

-.009 

.140 

-.004 

334.891 

275.719 

3.057 

-5.942 

353.009 

34 A • 2 06 

3.636 

-5.376 

3.249 

-.600 

3.091 

-1.336 

.230 

.011 

.124 

-.079 

.133 

-.110 

.175 

-.149 

79.016 

77.576 

61.961 

61. 605 

2.062 

.133 

• 160 

-.128 

.174 

-.129 

.217 

-.076 

.151 

-.153 

,005 

-.095 

13.570 

7.170 

.903 

.627 

.846 

. 591 

.670 

.612 

11.361 

11.282 

7.140 

4.794 

-1.936 

-9.759 

.076 

.073 

• 166 

,069 

.140 

,101 

-6.750 

-9.551 

-8.67B 

-9.092 

4,710 

3.306 

.022 

-.350 

9,546 

.705 

64.551 

49.841 

64,346 

49.625 

1.162 

,418 

• 406 

« 006 

7529228.122 

7522746.054 

72.649 

72.706 

-117.978 

-117.966 

35.754 

35.579 

359.961 

• 026 

6.200 

6.039 

7.474 

-.876 

101.162 

07,945 

2.175 

2.103 

0.464 

Z.357 

60.593 

19.506 

20.659 

-21.728 

61.124 

19,063 

303.845 

245.044 

123.645 

65.044 

303.645 

245.044 

303.845 

245.044 

109.065 

91.750 

107.536 

90.222 

111.063 

92.927 

71,366. 

-.779 

56.331 

.000 

12.332 

-9.102 

12.209 

-0.1B4 

11.994 

-8.519 

9.633 

-9.468 

10.724 

-10.3B7 

6 . 489 

-11.962 

12.606 

-6.420 

11.796 

-9.297 

16.661 

-8.046 


, 36410 40737.60660 


-.00271 

.02910 

.99544 

1.00084 

.00306 

,01370 

.06305 

•06569 

.00142 

•02966 

349,36021 

349.39191 

-2.21736 

2.61189 

347.99832 

346.00093 

-1.64475 

2.14194 

1.00664 

1.02644 

1.00569 

1.03030 

•06404 

.07038 

.00096 

.01789 

-.01274 

.01816 

.00669 

,02546 

76.47828 

70.47957 

61.77957 

61.77957 

1.00037 

1.00562 

.00055 

.01438 

.01674 

.02006 

.02615 

.03163 

• 01 U 00 

.02465 

•00266 

.01892 

10.29466 

10.40512 

.74751 

.74855 

,71422 

.71521 

.73949 

.74053 

11.33640 

11.33649 

5,91622 

5.94657 

-7.46468 

7.60921 

,07456 

.07457 

.13343 

.13618 

.11500 

•11646 

-9,17613 

9.17961 

-8.94551 

0.94639 

4.55243 

4.55381 

-.33007 

.33099 

9.16363 

9.16672, 

63.69992 

63.90119 

63.61769 

63.61871 

•44116 

.44144 

.23535 

.23548 



72.61735 

72.81736 

-117.98276 

117.98276 

35.66524 

35,66528 

5.76641 

36.00202 

6,10577 

6,10582 

3.66626 

4,16317 

95.43015 

95.51602 

2.13595 

2,13597 

5,39005 

5.56526 

42.50033 

42.84876 

-4.07377 

7.57962 

43.23770 

43,51399 

275.69601 

276.05272 

95.69606 

96.34450 

275.69606 

276.05277 

275.09606 

276.05277 

100.93810 

100.97547 

99.23665 

99.27276 

301.47274 

101.50925 

32.43672 

33.74127 

34.34776 

35.5U321 

•coooo 

3.03782 

•00000 

3,03101 

.00000 

3,07700 

•00960 

3.00277 

,00890 

3.03892 

.00956 

3.06967 

.00545 

2.23956 

.00060 

2.15312 

.00692 

2.22691 


50.86447 

0156 

.02906 

6156 

,10361 

6156 

.01336 

6156 

.01643 

6156 

.02963 

6156 

1.60740 

6156 

1.30040 

0156 

1.34044 

0156 

1.37220 

0156 

.19034 

0156 

.22365 

8156 

.02920 

6156 

,01787 

6156 

.01297 

0156 

i . .02452 

0156 

.45001 

8156 

.01320 

6156 

.10264 

6156 

.01437 

6156 

.02252 

0156 

.01442 

8156 

.02253 

6156 

.01873 

6156 

1,51067 

8156 

.03939 

8156 

.03759 

8156 

.03936 

8156 

.01311 

6156 

.57999 

0156 

1.47313 

6156 

.00065 

8156 

.02722 

6156 

.01235 

61 56 

.25201 

8156 

.12566 

• 156 

.11232 

0156 

.00895 

6156 

.23797 

6156 

.40372 

• 156 

,36006 

8156 

,01591 

6156 

.00781 

8156 

1883.47666 

8156 

.01034 

6156 

.00249 

0156 

•05000 

0156 

35.53582 

0156 

.02499 

015 6 

1 .97253 

0156 

4.09644 

0156 

.00929 

0156 

1.36553 

0156 

5.39092 

0156 

6.39219 

0 1 5 1 

4 , B95 30 

0156 

9.27708 

0156 

9,27710 

0156 

9.27710 

6156 

9.27710 

0156 

2.74714 

6156 

2.67007 

6156 

2.72263 

6136 

9.2B512 

6156 

8.98441 

6136 

3.03601 

6156 

3.03120 

8156 

3.07718 

0156 

3.00294 

8156 

3.03910 

8X56 

3.06984 

6156 

2.23969 

1136 

2,15323 

«1 56 

2.22704 

• 1 5t 
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Figure A. 37. Time histories of gust velocities, gust velocity differences, and aircraft' 
normal accelerations, Flight 31, Run 8. 




Figure A. 37. (continued) 



(m/s ) 
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Figure A. 37. (continued). 



TABLE A. 15. Average Turbulence Parameters, Integral Length Scales, and 

Correlation Coefficients of Gust Velocities, Flight 31, Run8. 


Mean Airspeed (m/s): 

4. Integral Length Scale (m) 


v L 

VC 

V R 

l uR 

l vR l wR 


103.22 

101.05 

102.86 

306.7 364.4 232.9 


Standard Deviation of 
Gust Velocities (m/s): 

l uRL l vRL l wRL 

302.5 380.7 249.6 


°uR 

a vR 

a wR 




3.93 

5.17 

2.52 

5. Correlation Coefficient of Gust 
Velocities: 

°uC 

a vC 

ff wC 

ti R u l/< j u r ct u l 

W/°v R °v L 


3.89 

5. 18 

2.36 




0.89 

0.81 

0.80 

c uL 

a v L 

CT wL 

II R v R/ a u R a v R 

v R w R/°v R a WR 

w R u R/ a w R CT u R 

3.89 

5.20 

2.42 




-0.20 

0.15 

0.00 

Standard Deviation of Gust 
Velocity Differences (m/s): 

tj R v l/°u R ct v l 

-0.19 

'V R WL/a VR a WL 

0.10 

w R u l/ °w r ct u L 
0.03 


a AuRL 

°AvRL 

ct awRL 

1.22 

1.08 

1.31 
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Probability Probability 


Gaussian 

Non-Gaussian 




Gust Velocity Difference Distributions 


Figure A. 38. Probability density functions for gust velocities and 

gust velocity differences (normalized with the standard 
deviation), Flight 31, Run 8 (r = degree of non-Gaussian) 
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Spatial Lag (m) 

a. One- and two-point common component correlations. 

Figure A. 39. Comparison of normalized one- and two-point correlation functions for 
gust velocities with theoretical models. Flight 31, Run 8. 



Spatial Lag (m) 

One- and two-point cross-correlations 
figure A. 39. (continued). 






a. One- and two-point common component spectra. 

Figure A. 40. Comparison of normalized one- and two-point spectral density functions 
for gust velocities with theoretical models, Flight 31, Run 8. 
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One- and two-point cross-spectra. 
Figure A. 40. (continued). 
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TABLE A.l 6 . 


List of All Parameters Measured and Their Range of Values, 
FI i ght 31 , Run 8 . 


CHANNEL 

UNITS 

1 TINE 

SECONDS 

? PHI DDT 

RAD/SEC 

- ACC L N CG . 

G UNITS 

- THE T A -DOT . — 

RAD 7 S EC 

' THETA. . . 

RAD 

6 PHI _ 

RAD 

' PSI 1 

RAD 

P DEL PSI I 

DEG 

5 PSI 2 

RAD 

ID DEL PSI Z_ 

DEG 

11 ACCl M LT . 

G UNITS 

\2 ACCL N RT 

G UNITS 

13.ACCL X.CG 

G UNITS . 

14 ACCL ¥ CG 

G UNITS 

15 ALPHA C T R 

RAD 

16 BETA CTR 

RAD 

17 TEMP I 

DEG F 

18 TEMP P 

DEG F 

W ACCL l INS 

G UNITS 

70 ALPHA RT 

RAD 

71 BETA RT 

RAD 

22 ALPHA IT 

RAD 

23 BETA LT 

RAD 

24 PSI DOT 

RAD/SEC 

25 T ENP TOT 

DEG C 

26 QC LT 

PS 1 0 

21 QC CTR 

PSIO 

?B0C RT 

P5I0 . _ 

29 PS 

PSIA ... 

30TEHP m 

VOLTS 

31 HYG RDM . 

DEG C. _ _ 

32 QC2 LT 

PS I D 

33 QC2 CTR 

PSID 

34 0 C 2 RT 

PSID 

35 DAR 

DEG 

36 OAL 

DEG 

37 DELEV 

DEG 

38 DSTAB 

DEG 

33 DRUD 

DEG 

40 DTHRR 

PCT MAX 

41 DTHRL 

PCT MAX 

42 OFLP "" 

POSITION 

43 058 

POSITION 

44 0 TO 6 " 

METERS 

458 TO 0 

DEGREES 

46 LONG 

DEGREES 

41 l AT 

DEGREES 

4b TRK ANG 

DEGREES 

4 # HDG „ . 

RADIANS 

bo VE. 

K/SEC . 

51 VN 

fl/SEC 

S? ALTITUDE 

KM 

53 TEMPC _ . 

DEGREES C 

M EW WNO SPO 

KNOTS 

5SNS HNO SPD 

KNDTS 

56 VINO SPEED _ 

KNOTS 

57*VIND OIREC 

“ DEGREES 

AA WIND DIR2 

DEGREES 

Kq WIND 0IR3 DEGREES 

60 WIND DIR4 

DEGREES 

6 ] AIRSPEED R 

M/SEC 

r? AIRSPEED C 

M/SEC 

l\ airspeed l 

M/SEC 

6 J DELTA ALT 

METERS 

6 A INRTL OISP 

METERS 

M UG RIGHT 

H/SEC 

67 U$ CENTER. 

H/SEC 

fcftUG LEFT._ 

M/SEC 

69 V G RIGHT 

H/SEC 

70 VG CENTER 

M/SEC 

VG LEFT 

M/SEC 

?? VG RIGHT 

h/SEC 

73 WG CENTER 

H/SEC 

74 VG LEFT 

H/SEC 


HIGH 

LOW 

MEAN 

RMS 

STD 

POINTS 

41204 .439 

40978.339 

41091.48890 

41091 . 54066 

65.22233 

9037 

.311 

-.284 ; 

-.00291 

.04454 

.04444 

9037 

1.963 

.336 

1.00164 

1.01212 

.14326 

9037 

.103 

-.101 

.00347 

.01625 

.01566 

9037 

.245 

-.059 

.06134 

.08022 

.05170 

9037 

.113 

-.190 

,00356 

.03967 

.03931 

9037 

190.046 

179.129 

183.98462 

183.99642 

2.08406 

9037 

5.931 

-4.779 

.07258 

2.05172 

2.03055 

9037 

547.691 

537. 3B3 

542.02722 

542.03089 

1.996B6 

9037 

6.284 

-4.448 

,40391 

2.09569 

2.05651 

9037 

3.466 

-2.046 

1.C1449 

1.03546 

.29122 

9037 

4.347 

-1.631 

1.01139 

1.05700 

.30717 

9037 

• 234 

-.072 

.06389 

• OB 1 53 

•05066 

9037 

.no 

-.129 

-.00294 

•02241 

.02222 

9037 

.103 

-.106 

-.01853 

,02582 

.01795 

9037 

.107 

-.190 

,00385 

,03193 

.03170 

9037 

76.497 

75.418 

76.06562 

76.06602 

.24603 

9037 

61.961 

61.613 

61.78630 

61.70659 

.02734 

9037 

2.004 

.288 

1,00677 

1 .01739 

.14657 

9037 

.124 

-.106 

-.00610 

•02196 

.02109 

9037 

.192 

-.173 

•01323 

.03237 

.02954 

9037 

.131 

-.064 

« 02 11 B 

•02615 

.01853 

9037 

.161 

-.168 

.00708 

.02982 

.02897 

9037 

.069 

-.103 

•00229 

.02379 

,02368 

9037 

12.192 

8,057 

10.24093 

10.28370 

.93706 

9037 

1.033 

.619 

.78696 

.79234 

.07306 

9037 

.966 

• 601 

.73379 

,75707 

,07038 

9037 

1.019 

.622 

.78172 

.78516 

.07360 

9037 

11.613 

11.297 

11.50133 

11.5C185 

.08570 

9037 

6.9 07 

1.911 

5.78063 

5.86605 

.99750 

9037 

-.178 

-8.976 

-6.30606 

6.46187 

1.41032 

9037 

.071 

.068' 

.06942 

.06942 

,00044 

9037 

.177 

.140 

•16480 

•16511 

.01011 

9037 

.163 

• 132 

.14798 

.14830 

.00973 

9037 

-0.722 

-9.330 

-8.97011 

8.97170 

•16830 

9037 

-8.403 

-0.631 

-Q. 47336 

8.47367 

.07261 

9037 

4.953 

4.682 

4.79274 

4.79341 

.07993 

9037 

-.347 

-.360 

-.33328 

•35330 

•00331 

9037 

10.117 

9,773 

9.94806 

9.94844 

.08661 

9037 

66.016 

63.137 

63.63407 

65.63443 

.22374 

9037 

65.916 

64.844 

65 « 46662 

65.46731 

•29866 

9037 

.395 

.348 

•36829 

• 36853 

.01324 

9037 

.265 

,260 

•27283 

.27293 

.00721 

9037 

7326163.491 

7320103.629 



1731.03911 ’ 

9037 

72.809 

72.797 

72.84398 

72.64399 

.02706 

9037 

-117.961 

-117.992 

-117.96349 

117.98549 

.00369 

9037 

33.700 

3 > • 4 67 

35.59327 

35.59333 

•06252 

9037 

161.966 

174.296 

177.68239 

177.69231 

1,87781 

9037 

3.331 

3.137 

3.22416 

3.22437 

.03732 

9037 

10.753 

-3.436 

4.39653’ 

5.68080 

3.59769 

9037 

-94.763 

-118*049 

-104.88967 

105.09696 

6.59786 

9037 

2.164 

1,942 

2,01994 

2.02064 

,06023 

9037 

6.914 

3.293 

5.14864 

5.19791 

,71407 

9037 

49.366 

2.533 

23.79845 

25.73772 

9.85417 

9037 

23.820 

-27.249 

-8.19980 

11.33378 

7,82462 

9037 

51.412 

3.724 

23.99334 

28.14096 

10.78237 

9037 

330.181 

223.791 

288.65429 

2B9. 01875 

14.51058 

9037 

150.181 

45.791 

108.63435 

109.61690 

14.51059 

9037 

330.181 

225.791 

288.65435 

269.01680 

14.51059 

9037 

330.181 

225,791 

288.65433 

289.01680 

14.51039 

9037 

116.833 

92.210 

102.86079 

102.96626 

4.63936 

9037 

115. 066 

90.199 

101.05040 

101.15233 

4.54017 

9037 

117.624 

91.603 

103.32909 

103.43172 

4.60668 

9037 

109.099 

-113.141 

-35.17826 

69.74513 

60.22663 

9037 

97.056 

-117.254 

-39.27849 

70.64242 

58.71910 

9037 

" ~ 9,631 

-18.056 

.00000 

3.93072 

3.93093 * 

9037 

9,630 

-16.490 

.00000 

3.69204 

3.89225 

9037 

10.556 

-15.009 

•00000 

3.89019 

3.89041 

9037 

11.370 

-12.709 

,07236 

3.17998 

5.17977 

9037 

11.207 

-13.431 

•07506 

3.18344 

3.18318 

9037 

11. 347 

-12.009 

.06633 

3.20477 

5,20461 

9037 

11.309 

-14,749 

.00751 

2.52506 

2.32319 

9037 

10.991 

-10.188 

•01268 

2.36776 

2.36786 

9037 

10.205 

-10.865 

•01106 

Z. 42446 

2.42457 

9037 


?40 



Date: Nov. 29, 1982 

Start Time: 11 : 29 : 30 (PST) 

Duration: 338.0 seconds 
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Figure A. 42, (continued). 


TABLE A. 17. Average Turbulence Param 
Correlation Coefficients 

1. Mean Airspeed (m/s): 

\ Vc Vr 

103.15 100.84 102.84 

2. Standard Deviation of 
Gust Velocities (m/s): 


a uR 

Q 

< 

70 

a wR 

4.10 

5.12 

2.40 

a uC 

a vC 

a wC 

4.10 

5.10 

2.21 

r uL 

a vL 

a wL 


4.18 5.12 2.34 

3. Standard Deviation of Gust 
Velocity Differences (m/s): 

ct auRL a AvRL c awRL 
1.35 1.16 1.45 


iters. Integral Length Scales, and 
of Gust Velocities, Flight 31, Run 9. 


Integral 

Length Scale (m):' 

l uR 

l vR 

l wR 

327.8 

338.0 

93.9 

l uRL 

l vRL 

l wRL 

341.5 

338.0 

83.9 


5. Correlation Coefficient of Gust 
Velocities: 

u R u L/ CT UR CT Uj_ v R v L/ a vp CT V|_ wr w l/ Owr ct W[_ 

0.80 0.90 0.80 

tt R v r/ ct u r ct v r VR^R/q VR q WR WrPr/pwrPur 
0.30 . 0.20 0.19 

URVL/ q UR CT v L VRWL/q VR g WL WRUj:/a WR cT UL 
0.30 0.20 0.18 
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Probability Probability 


Gaussian 

Non-Gaussian 


r = 2.5 r = 1 .0 r = 1 .0 



Gust Velocity Distributions 



Gust Velocity Difference Distributions 


Figure A. 43. Probability density functions for gust velocities and 

gust velocity differences (normalized with the standard 
deviation). Flight 31, Run 9 (r = degree of non-Gaussian) 
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One- and two-point cross-spectra. 
Figure A. 45. (continued). 
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OP POOR QUALITY 


TABLE A. 18. List of All Parameters Measured and Their Range of Values, 
FI i ght 31 , Run 9. 


CHANNEL 

UNITS 

1 tine 

SECONDS 

? PHI DOT 

RAD/SEC 

' AC C l N CS 

G UNITS 

< THETA DOT 

RAD/SEC 

5 THETA 

RAD 

6 PHI 

RAD 

7 P 5 1 1 „ 

RAD 

6 DEL PSI 1 

DEG 

9 PSI 2 

RAD 

10 DEL PSI 2 

DEG 

n ACCl N IT 

G UNITS 

12 ACC l N R r 

G UNITS 

nACCi x ct 

C UNITS 

19 ACC L Y CG 

G UNITS 

15ALPHA ctr 

RAO 

lb BE T A CTR 

RAO 

17'TEHP I 

DEG F 

18TEHP P ' 

DEG F 

19 ACCL 2 INS 

G UNITS 

20 ALPHA RT 

RAD 

21 BETA RT 

RAD 

22 ALPHA LT 

RAD 

23 B ET A LT 

* RAD 

?4 PSI DOT 

RAD/SEC 

25TEHP TOT 

DEG C 

26 QC IT 

PSID 

27 OC CTR 

PSID 

23QC RT 

PSID 

29 p $ 

PSI A 

10TENP IRT 

VOLTS 

31HYGR0N 

DEG C 

3 ? QC 2 LT 

PSID 

330C2 CTR J . 

PSID 

34 QC 2 RT 

PSID 

3S0AR 

DEG 

36 DAI 

DEG 

37 DEL IV 

DEG 

33 OST AB 

DEG 

39 DRUD 

DEG 

♦odthrr * 

PCT HAX 

41DTHRL 

PCT MAX 

42DFLP 

POSITION 

430*8 

POSITION 

44 D TO 6 

METERS 

456 TO 0 

DEGREES 

46 LONG 

OEGREES 

47 

DEGREES 

46TRK ANG 

DEGREES 

49HDG 

RADIANS 

60 VE 

H/SEC 

SI VN 

H/SEC 

s ? altituoe 

KH 

&3TEHPC 

DEGREES C 

54 EV WND SPD 

' KNOTS 

55 N5 WHO 5 PD 

KNOTS 

56WINO SPEED 

KNOTS 

57 WIN D DXREC 

DEGREES 

58WIND DIR2 

DEGREES 

59 WIND DIR3 

DEGREES 

60 nIND 0IR4 

OEGREES 

6) AIRSPEED r 

N/SEC 

62 AIR5PEED C 

M/SEC 

6 3 AIRSPEED L 

‘ N/SEC 

64 DELTA ALT 

HETERS 

bSlNRU DISP 

HETERS 

66 UG RIGHT _ 

' H/SEC 

67 UG CENTER 

H/SEC 

6 BUG IIM 

H/SEC 

69 VC RIGHT 

" N/SEC 

70 VC CENTER 

N/SEC 

71 VG LEFT 

N/SEC . 

7?VG RIGHT 

M/SEC 

73 WG CENTER 

H/SEC 

74 WG LEFT 

N/SEC 


HIGH 

LOW 

41704.552 

41369.602 - 

.318 

-.435 

1.921 

-.517 

.112 

-.089 

.168 

-.036 

.152 

-.260 

220.340 

207.660 

5.255 

-7.264 

581.839 

369.517 

5.696 

-6.888 

3.334 

-1.148 

<i.001 

-1.443 

.178 

-.012 

.102 

-.104 

.087 

-.174 

.164 

-.141 

76.857 

IB. 021 

61.961 

47.322 

1.931 

-.390 

. 085 

-.166 

.170 

-.120 

.144 

-.152 

.143 

-.123 

.097 

-.093 

14.653 

10.203 

1.074 

.651 

• 9 06 

.618 

1.039 

.642 

12.509 

12.119 

6.063 

4.319 

4.123 

-2.329 

.079 

,074 

.179 

• 141 

.169 

.134 

-9.910 

-10.242 

-9.202 

-9.575 

6.066 

5.345 

-.374 

-.422 

11.483 

10.402 

67.969 

66.406 

67.773 

66.406 

.320 

.273 

.324 

.299 

7506823.590 

7464351.806 

72,886 

72.861 

-118.062 

-118.248 

35.355 

35.103 

216.361 

207,622 

3.869 

3.644 

-42.984 

-60.999 

-76.290 

-90.819 

1.597 

1.338 

9.502 

5.133 

39.790 

-18*404 

23.456 

-24.192 

44.252 

• 125 

359.936 

.142 

179.936 

-179.658 

359.936 

.142 

843.475 

-533.217 

114,681 

90.490 

111.755 

88.836 

116.498 

91.024 

251.146 

-7.522 

240.544 

-10.792 

11.496 

-11.949 

11.356 

-11.386 

11.359 

-11.158 

13.853 

-22.090 

12.866 

-17,501 

13.763 

-10.901 

13.700 

-14.128 

13.996 

-10.785 

15.640 

-13.421 


KEAN RH5 

1537,07670 41537.18926 

-.00239 .03126 

I. 00123 1.01562 

.00301 .01688 

.05660 .06526 

-.00262 .04819 

213.51567 213.52714 

-1.33276 2.55017 

575.30539 375.30928 

-.92423 2.37329 

i;oi390 1.06483 

1.01099 1.06704 

.05691 .06566 

-.00127 ,02394 

-.02678 .03247 

.01160 .03349 

76.52245 76.52426 

61.87709 61.87728 

1.00521 1.01953 

-.00918 .02266 

.02013 .03324 

.01113 .02246 

.01432 .03243 

•00253 .02486 

12.68530 12.71776 

.63864 .64041 

•80057 .60229 

.63345 .83524 

12.39236 12.39270 

7.15658 7.17453 

2.63444 2.15294 

.07689 .07691 

,16947 .16976 

.13502 .15543 

-10.00441 10.00472 

-9.40686 9.40722 

5.76361 5.76764 

-.40008 .40033 

10.91013 10.91478 

67.12432 67.12593 

67.02022 67.02126 

.29930 .29939 

.31140 .31147 

72.87323 72.67323 

-116.15379 116.15360 

35.23129 35.23137 

211.17269 211.1840b 

3.75054 3.73074 

-30.53699 50.71863 

-83.45358 83.33606 

1.41511 1.41641 

7.61940 7.67039 

II. 90531 16.43143 

-.25178 3.88026 

15.95477 . 17.45192 

240.11266 231.21631 

60.11270 93.23699 

240.11270 251.21834 

178.96187 281.19140 

102.84471 102.69615 

100.84162 100.89187 

103.13851 103.20907 

69.33534 92.13310 

• 63.39001 89.00817 

.00000 4.10906 

.00000 4.10739 

•06000 4.16736 

-.01121 5.12140 

-.00555 3.10793 

-.00246 5.12131 

• 023M 2.40762 

.02626 2.21561 

.02633 2.34861 


STD 

POINTS 

96.70255 

1 3399 

.05120 

13399 

.17035 

13399 

.01661 

13399 

.03248 

13399 

.04812 

13399 

2.21339 

13399 

2.17427 

13399 

2.11560 

13399 

2.16819 

13399 

,32540 

1 3399 

.34128 

13399 

,03279 

13399 

.02391 

13399 

.01835 

13399 

•03142 

13399 

•52620 

13399 

.15120 

13 399 

.17030 

13399 

.02072 

13399 

.02891 

13399 

.01950 

1 3399 

.02910 

133 99 

.02473 

13399 

.90324 

13399 

.05450 

133 99 

,05240 

1 3399 

.05464 

13399 

.09177 

13399 

.50722 

13399 

1.40119 

13399 

.00150 

13399 

.01004 

13399 

•01115 

1 3399 

.07867 

13399 

•08209 

13399 

.21015 

13399 

. Cl 3 87 

13399 

•31864 

13399 

.46708 

1 3399 

~ ~ .37408 * 

13399 

.01310 

13399 

.00676 

13399 

7011.66724 

1 3399 

•00719 

13399 

•05331 

13399 

.07257 

13399 

2.19175 

13399 

.03910 

13399 

4.04635 

13399 

3.7X141 

1 3399 

•06068 

13399 

.68298 

13399 

11.32521 

13399 

5.67309 

13399 

7,07236 

13399 

73.67116 

13399 

73.87120 

1 3399 

73.87120 

13399 

216.89805 

13399 

3.25334 

13399 

3.16411 

13399 

3.23036 

13399 

60.67753 

13399 

60.38935 

1 3399 

4.10921 

13399 

4.10755 

13399 

4.18752 

1 3399 

5.12156 

13399 

5.10812 

13199 

5.12170 

1 3399 

2.40779 

1 3399 

2.21576 

13399 

2.34875 

13399 
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Figure A. 47. (continued) 





TABLE A. 19. Average Turbulence Parameters, Integral Length Scales, and 

Correlation Coefficients of Gust Velocities, Flight 31, Run 10. 


1. Mean Airspeed (m/s): 

\ Vc Vr 

117.27 115.20 116.70 

2. Standard Deviation of 
Gust Velocities (m/s): 


ff uR 

p vR 

ff wR 

2.04 

4.57 

2.40 

a uC 

ff vC 

a wC 

1.99 

4.58 

2.34 

a uL 

a vL 

a wL 

2.02 

4.61 

2.33 


3. Standard Deviation of Gust 
Velocity Differences (m/s): 

ff AuRL a AvRL a AwRL 
0.41 0.31 0.38 


4 . Integral Length Scale (m):- . 

*-uR L vR l wR 

641.8 729.7 832.2 

l uRL l vRL l wRL 

638.3 742.9 863.8 

5. Correlation Coefficient of Gust 
Velocities: 

lj R u l/ <t ur (I ul 
0.98 

TJ R v R /g UR gv R 
0.08 

¥[/ q U|j g v L 

0.09 


VRVL/a VR a VL 

w R w L/ a w R a W|_ 

0.99 

0.98 

V R W R/°VR a WR 

WRUR/a WR a UR 

0.00 

-0.47 

VRWL/a VR a WL 

w R u l/°w r °u l 

-0.01 

-0.45 
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a. One- and two-point common component spectra. 

Figure A. 50. Comparison of normalized one- and two-point spectral density functions 
for gust velocities with theoretical models, Flight 31, Run 10. 
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One- and two-point cross-spectra. 
Figure A. 50. (continued). 
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TABLE A. 20. List of All Parameters Measured and Their Range of Values, 
Flight 31 , Run 10. 


THAN* EL 


UMTS 


HIGH 


LOW 


MEAN 


STD 


POINTS 
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Date: Nov. 29, 19G2 

Start Time: 11 : 46 : bl (P 

Duration: 334.0 seconds 
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Flight altitude and horizontal wind along flight path. Flight 31, Run 11 
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Figure A. 52. Time histories of gust velocities, gust velocity differences, and aircraft's 
normal accelerations, Flight 31, Run 11. 
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Figure A. 52. (continued). 


s, Integral Length Scales, and 
Gust Velocities, Flight 31, Run 11. 


Integral 

Length Scale (m): 

l uR 

L V R 

<-wR 

370.0 

246.1 

203.3 

l uRL 

l vRL 

l wRL 

375.6 

241.7 

193.1 


Correlation Coefficient of Gust 
Velocities: 


u R u l/ g ur ct ul 

VRVj:/o V Rav L 

w r w l/owr ct W|_ 

0.83 

0.81 

0.78 

urvr/^ur^vr 

VRWR/a VR a WR 

w R u r/°wr ct u R 

-0.21 

-0.09 

0.48 

u R v L/ ct ur^v L 

VRWL/a VR aw L 

w R u L/ a WR CT U[_ 

-0.28 

-0.10 

0.40 


67 



Probability Probability 


Gaussian 

— Non-Gaussian 


r = 2.5 


r = 1.5 


r = 1 .5 




Figure A. 53. Probability density functions for gust velocities and 
gust velocity differences (normalized with the standard 
deviation), Flight 31, Run 11 (r = degree of non-Gaussian) . 
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One- and two-point cross-spectra. 
Figure A. 55. (continued). 
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TABLE A. 22. 


List of All Parameters 
FI i ght 31 , Run 1 1 . 


Measured and Their Range of Values 


CHANNEL 

i yi*f“ 

UNITS 

' seconds 

HIGH 

42744.524' 

LOW MEAN RMS 

42 41 0. '52 4'“4 2577. 524 00 4 2 577.63 319 

STD 

96.42831 

POINTS 

13361 

2 PH 1 t> 0 T 

radTTK 

.243 

-.297 

-.00258 

.04193 

.04185 

13361 

3 A C CT~ N "C G 

“■« ‘UNITS 

~ ” 1.921 

.282 

1 • 0003d 

1.00958 

,13652 

13 361 

4 THETA DOT" 

~ ft A J / S EC 

.C8C 

-.069 

.00323 

.01449 

.01413 

llUl 

5 ThFTa 

RAO 

.231 

-.060 

.04493 

.07771 

.06340 

13361 

6 phT 

RAO 

• 1 91 

-.151 

.00474 

.03775 

.03745 

1 3 361 

?psi“i 

RAD 

""’317.9U7 

300.045 

312.03280 

312.03864 

1.90874 

13361 

8 DEL P7T 1 



1.758 

-5.939 

-2.2MA5 

2.59711 

1.37475 

1 3 361 

9 psr? 

RAD 

317.446"' 

309.708" 

313.41944 

“313.42230 

1.33 612 ___ 

13361 


i o p e l~ ps rr 

nirn u lt 


' DEG~ 


12 ACC l 
13aCCI' 


N 91 " 


imTr 

■ a" 0^41 ti 


2.315 


i4ir:L t zz 

i5AiPHi"r*nr 
16BEF1 r 


i ce — fi^uMiTr" 


G UNITS 


riir 

“RID" 


17TFKT 5 I 

DES F 

~77; 9 — 

lemr - ? 

OEC T 

60.522 

19 aC C l l f NS 

6 UNITS 

2.026 

20 ALPHA "RT 

RAD 

.117 

216ET* AT 

RAO 

.102 

2? AL"P HA LT 

"RAO 

' .103 

23BETA il 

raD 

.095 

24PSI DOT 

ftAD/SEC 

.062 " 


26 OC LT" 

27 QC" CTR 

TSTC 
PSID 

— “T.nr 

. 1.073 

28 QC R T 

PS 10 

l.ni 

29PS 

957* 

12.360 

30TENP IRT 

VDLTS 

6.907 


"DEGT 

5.296 


32QC2 LT 



3ADCF77 

PSID 

,157' 

.04 8“ 

•10622 

.11192 

.03527 

13361 


“DFG ‘ ■ 

-5.FST" 

-6.595 

“ ^5795484 

5.96258 

.30362 

13361 

360 * L 

DtG 

-5.455 

-6. *74 

-5.86257 

5.87255 

.34235 

13361 

370ELEr 

DEC 

" 5.150 

'4.721' 

4.94840 

4,95045 

*14219 

1 3)61 

38DST Afl 

DEG ‘ 

-.iff 

-.401” 

—.30076 

”™“ 73PFBF 

.00900 

13361 

390RUD ' 

DEG 

*" 11.768' 

1 G , 629 

11.14561 

11, 15003 

.31 395 

1 3 Sol 

40DT9RR 

' PC T~HAX ” 

69.629 

67.773 

60. 66135 

68.66343 

.5)496 

— ui4 

41DTHRL 

PCT HA* 

7U.02o 

bTTSTT 

697TWUF 

69.34635" 

.42644 

13 361 

42DFLP 

DaSITTITK^" 

.199" 

.160“ 

• 17571 

• 177T10 

•01173“ 

1 3 36l 

43TJSB 

“PD^rrroH 

.3 85“ 

. 377“ 

. 3 8170 

•3B130 

.7581 95 

i 3T5T 

440 TU 5 

PETERS 

MrAI iftAAlii tig 

iflamaa 



361). 41 976 ” 

1 3 3 6 1 

45B'TD 5 

DEG«teE5 

iMMiWi ft ill 

— T2.68r 

72.H0349 

72.80352 

*06 B44 

1 3 361 

46113 'TG 

DETl^TES 

-118. lA?" 

-U 87408' 

-118,29182 118.29183 

.06476 

13361 

47 A AT 

DEGREES 

35.279 

35.068 

35.17133 

35.171)8 

.06122 

13361 

48 T RK ANi 

■ nun — 

323,207 

315.629 

319,11203 

319.11605 

1,60058 _ 

13)61 



49H0G 
60 Vt 
51 W 

52 ALTTTUO 
53 UHEC" 

54 EW WHtTW 

55 IT 


56 WIMP 5 PEEP 

57 WH O OlftTC 
SflVTfTCTUTR 


... hTScc 

^Eim' 

fcs BilOiig — Mlljl 


£5 UU ftii-IT hTTkC 

67 WTtyre; — jr/ 5 rrzi 

fcfl UG LEF T H/SEC 

69 V w A I GHT STSIT 

70 T(TTT R f Eft H/ Si C 

7i ~vrriFT 


iiuiir— 

94 . $ 7 * 

lu 7. 01063 

1077 T 7339 

5.90453 

1 3361 

4 7^601 

-SF.^oF" 

176.23842 

176.29089 

248.54684 

? 49 ; 2 l 27 T^ 

1 > 5. 26792 

1 33 61 
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Figure A. 57. Time histories of gust velocities, gust velocity differences, and aircraft' 
normal accelerations. Flight 31, Run 12. 
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Figure A. 57. (continued). 
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TABLE A. 23. Average Turbulence Parameters, Integral Length Scales, and 

Correlation Coefficients of Gust Velocities, Flight 31, Run 12. 


1. Mean Airspeed (m/s): 4. Integral Length Scale #1 (m): 


Vl 

VC 

Vr 

l uR 

1-vR 1 -wR 


101.03 

98.55 

100.56 

127.7 252.6 202.4 





l uRL 

L vRL l wRL 


Standard Deviation of 
Gust Velocities (m/s): 

137.5 250.1 190.8 


a uR 

ff vR 

a wR 




1.68 

1.43 

1.66 

5. Correlation Coefficient of Gust 

Velocities: 



a uC 

ff vC 

ff wC 

W/*u R "u L 

VRVL/(J VR a VL 

w R w L/ a w R a W(_ 

1.68 

1.47 

1.54 




0.66 

0.81 

0.78 

a uL 

CT vL 

CT wL 




1.70 

1.42 

1.59 

URVR/a UR a VR 

Vrwr/ a VR a WR 

WROR/a WR CTu R 




0.00 

0.30 

-0.22 

Standard Deviation of Gust 

URVL/a UR a VL 

VRWL/a VR a WL 

WRq>WR°u L 

Velocity Differences (m/s): 

0.01 

0.31 

-0.20 

ct auRL 

^ avRL 

a AwRL 





0.90 0.74 0.91 
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Probability Probability 


Gaussian 

Non-Gaussian 


r = 1 .5 r = 1 .0 r = 1 .0 



Gust Velocity Distributions 


r = 1 .5 r = 1 .5 r = 1 . 5 



Gust Velocity Difference Distributions 


Figure A. 58. Probability density functions for gust velocities and 
gust velocity differences (normalized with the standard 
deviation), Flight 31, Run 12 (r = degree of non-Gaussian) . 
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Aircraft data 
von Karman 



a. One- and two-point common component spectra. 

Figure A. 59. Comparison of normalized one- and two-point correlation functions for 
gust velocities with theoretical models. Flight 31, Run 12. 
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Figure A. 60. Comparison of normalized one- and two-point spectral density functions 
for gust velocities with theoretical models. Flight 31, Run 12. 
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One- and two-point cruss-spectra. 
Figure A. 60. (continued). 
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TABLE A. 24. List of All Parameters Measured and Their Range of Values, 
FI i ght 31 , Run 1 2 . 


CHANNEL 

UNITS 

HIGH 

LOW 

MEAN 

RMS 

STD 

POINTS 

1 Tine 

SECONDS 

43146.508 

43009.608 

43078.05820 

4 307 0.07633 

39.53C45 

*** 1 

2 PHI DOT 

rad/TEG 

.085 

-.12?" 

-.00254 

.03036 

.03025 

5477 

- ACClTinG G UNITS 

4 TilCT* — firtt Sl'ftrt'CP 

1.361 

,1)40 

.656" 
-,<j4 5 

' .99898 

.00352 

1.00374 

.00997 

.09 762 
'i‘00933 

5477 

5477 

ThEYa 

RAD 

.080 

.019 

.04071 

.05036 

_.OI278_ 

5477 

6 PHI 

RAD 

.079 

-.115 

.00130 

.02959 

.02957 

5477 

7 psiT 

RAO 

125.590 

118. 898 

123.04541 

123.04936 

.90 859 

5477 

W DEL PSI 1 

DEG 

- • C64 

-6.618 

-2.50441 

2.69251 

.93878 

5477 

9 PSI - 2" 

RAD 

'487.491 " 
.253 

' ”48 1.154 
— 6 • ?4 1 

4B4. 99570 
' -2 71 4 91 4 

484,99675 

'2.36260 

‘.96953 

.98165 

”54 77 
- - 5477. 

11 aCCL n lt 

G UNITS 

1.695 

.291 

1 .01263 

1.03064 

.19101 

5477 

12 ACCl N'ftT 
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~ 1 • 905 

.iia 
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1. 02931 
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.004 
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. 05i 20 
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B fcn CTP 

V'Un? i 


UNITS 

"raO" 
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19 AC C L l INS 


_DEG_F_ 

ore f 


*017 
• 091 " 


76.317 


-.064 

-.058 


-.02879 
.0 0713 ' 


.03050 

*9153 L 


75.98198 


.01006 

-tails i_ 


3 A 77 

3 fill 


G UNITS 


20 ALPHA RT RAD 


1,429 



• 0 AO 
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5 A 77 


22 ALPHA Lf 
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RAD 
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. - .054” 

.009 

-.024' 

-.057 

.00979 

.010«2 

.0145 3 
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t.01015 

.01701 
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l54 77 

24-p$r&5f” 

rao/sKc 

• 062 

-.039 

10.420 

.00275 
1 i . 3 66 1 5 
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11.37100 
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5477 

5417 

26-yrrr 
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Date: Nov. 29, 1982 

Start Time: 12:01:21 (PST) 

Duration: 269.1 seconds 



Figure A. 61. Flight altitude and horizontal wind along flight path. Flight 31, Run 13 





Figure A. 62. Time histories of gust velocities, gust velocity differences, and aircraft's 
normal accelerations. Flight 31, Run 13. 
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Figure A. 62. (continued). 



TABLE A. 25. Average Turbulence Parameters, Integral Length Scales, and 

Correlation Coefficients of Gust Velocities, Flight 31, Run 13. 


1. Mean Airspeed (m/s): 

\ Vc Vr 

103.30 101.40 103.30 

2. Standard Deviation of 
Gust Velocities (m/s): 


*uR 

a vR 

a wR 

2.49 

5.57 

2.43 

a uC 

a vC 

a wC 

2.48 

5.57 

2.29 

a ul 

°vL 

a wL 

2.59 

5.56 

2.41 


3. Standard Deviation of Gust 
Velocity Differences (m/s): 

ct auRL ff AvRL ct awRL 
1.53 1.39 1.59 


4. Integral Length Scale #1 (m): 


*-uR lyR *-wR 

156.0 428.8 83.7 

l uRL l vRL l wRL 

148.6 424.4 82.6 


5. Correlation Coefficient of Gust 


Velocities: 

u R u L/ a UR a U|_ 

VRVL/a VR a VL 

w R w U/ a w R °w L 

0.80 

0.91 

0.79 

URVR/a U pC V R 

VRWR/a VR CT WR 

w R u r/°w r ct ur 

-0.18 

-0.32 

0.25 

TJRVL/a UR CT VL 

v R w L/ a VR cr W[_ 

w R u L/ a w R a ui_ 

-0.19 

-0.32 

0.22 
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Probability Probability 


Gaussian 

Non-Gaussian 


r = 1.5 r = 2.5 r = 1 .0 



Gust Velocity Distributions 


r = 1.5 r = 1 .5 r = 2.5 



Gust Velocity Difference Distributions 


Figure A. 59. Probability density functions for gust velocities and 
gust velocity differences (normalized with the standard 
deviation). Flight 31, Run 13 (r = degree of non-Gaussian) . 
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Aircraft data 
von Karman 







Aircraft data 
von Karman 
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a. One- and two-point common component spectra. 

Figure A. 65. Comparison of normalized one- and two-point spectral density functions 
for gust velocities with theoretical models. Flight 31, Run 13. 
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One- and two-point cross-spectra. 
Figure A. 65. (continued). 
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TABLE A. 26. List of All Parameters Measured and Their Range of Values, 
FI ight 31 , Run 1 3 . 


CHANNEL UNITS 

\ Tt*E secnuos 

? ph! urn m/5Ec 

3 iCCL U CG FTI T 5'~ 

4 THFTT'OOT — ffAD7Scr~ 


HIGH 

43*44 .$&<r 


LOW MEAN RMS 

' 43 * 63.534 4341$. 64760 <.3413711707 


STD 

7776m r 


1701686 
— - 0 



294 




















CO 



296 


f" 


Figure A. 67. Time histories of gust velocities, gust velocity differences, and aircraft' 
normal accelerations. Flight 31, Run 14. 
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Figure A. 67. (continued). 
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TABLE A. 27. Average Turbulence Parameters, Integral Length Scales, and 

Correlation Coefficients of Gust Velocities, Flight 31, Run 14. 


1. Mean Airspeed (m/s): 

V L V C V R 

103.38 101.07 102.99 

2. Standard Deviation of 
Gust Velocities (m/s): 


a uR 

a vR 

°wR 

2.51 

3.54 

2.37 

G uC 

ff vC 

CT wC 

2.47 

3.50 

2.12 

CT uL 

ff vL 

°wL 

2.52 

3.42 

2.28 


3. Standard Deviation of Gust 
Velocity Differences (m/s): 

a AuRL g avRL ct awRL 
1.29 1.12 1.37 


4. Integral Length Scale (m): 


luR *-vR *-wR 

174.9 204.4 66.8 

*-uRL *-vRL *-wRL 

161.3 205.4 64.5 


5. Correlation Coefficient of Gust 
Velocities: 


TOT/^ur^ul 

VRVL/a VR cTv L 

WrWl/ a WR a WL 

0.79 

0.90 

0.77 

u R v r/^ur q vr 

v R w r/°vr ct wr 

W R U R/ 0 WR a UR 

0.18 

0.19 

0.10 

u R v L/ a UR a V|_ 

Vr w l/<J V rCT W [_ 

w R d l/ g wr ct U[_ 

0.13 

0.27 

0.07 
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a. One- and two-point common component spectra. 

Figure A. 70. Comparison of normalized one- and two-point spectral density functions 
for gust velocities with theoretical models. Flight 31, Run 14. 
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One- and two-point cruss-spectra. 
Figure A. 70. (continued). 



ORIGINAL PAGE IS 

OF POOR Q T J a.tj'py 


TABLE A. 28. List of All Parameters Measured 
FI ight 31 , Run 14 . 


and Their Range of Values 


CHANNEL UNITS 

1 TINE SECONDS.. 

2 PHI DOT . - RAD/ 5CC . 

3 .iLZL. N . C C _ _G U NIT s _ 

* JtjETA . UQI F ADZ5 EX — 

5 THETA RAD 

i £HI_ 


7 P5L- l Rad 

8 DFL._P5.J ; PEG 

9 PSL 2 RAP 

TO DEI PS1 2 DEG 

11 A C CL N LT G UNITS 

1 ? iCCL H RT C UN I TS_ 

13 ACCL X CC C UNITS 

14 ACCL Y CG. G UNITS . . „ 

15 ALPHA CTR __.R AD 

16 PF 1 A CTR *AD 

17 T F HI P I . . tic F . . . 

18 TC HP_P PEG _F. . 

19ACCL 7 INS G UNITS 

20 AL p HA RT. . __RAD 

21 ♦'Ll*. R T PAD 

22 ALPHA LT SAD 

23 P F T A LT kAD.. . 

2 * P.5 I DDT RAU/SFI 

25 TFHP TuT PEL C 

26 OC \j PS ID 

27 QC CLR p $J D 

2BQC* RT P 3 ID 

29 p S P S I A 

30 TFHP. 1 R X V0U.S 

31 HYGROM DC GC 

32 OC2 LT P S 5 6 

33 D C ? CTR 


35 OAR 

36 DAL 


HIGH 

1 06 9. 4Z 3. 

. 195 

1. 88 2.. 

LOW 

.. „ _4 3660 . 57 A 
..-,280. 
*A> 4 7 . 

MEAN 

.43764 .9978U 
. _ .-.00275 
1*0064 1_ 

RMS 

43765. OJ934 
.05322 

1 . 02 2b3 

• 0 1 6b5 . ... 

.00263 . 
» 043 b5 

. 00 7_ 

. 112. 

^ 20.0 . 

t. 107 .00350 

_ Oi7 .05980 

-.197 ,0012 8. 

"206.681 

. 190.753 

200. P 701 *1. 

200.80491 . 

5,0 ?0 . 

r 16.08 L 

-.78389 

2.73068 ... 

.566.013 

5 4 8. / 4 7 

5b 3. 09 99 5_ 

__ 563. 10573 _ 

5. 179. 

r 1 0 « 3 b 7 . 

-.40259 

2.44204 

_ 3.661 

-1.016 

1.01093 

. _ 1.07647 

4 . 00 4 

-1,1/? 

1.01694 

_ 1 .OH022 . 

.1/ 9 _ • 0 2 6 

"'”1-"* 056 34* 

,06047 . 

. 107 

-.166 

.00031 

.02915 

...0/3. 

-.1 63 

-.03216 

" .0 3720... 

• 146 

-.191 

.01173 

,03909 

/«*• 4 97 

76.967 

.7,* 70-30568 

76.30566 



. „59. 11203. 

. .. 1.01012 
-.01699 

.Ol^l 

. • 00 *• 68 . 
.01401 

-0.0 2 .b 2 _. 

1 1 .64364 
... ,831bH. 

. 79435 .. 

_ .62 555. 

12.19375 
6. .90 7 ..J.3U7L 

„-fl .0 76 1699 3 

, 0^<l -06153.. 

. 118 4 


1 . 02 b 73. 
,02763 

.019 76 

.01950 . 
.0 3729 

..-01197.. 

_ 11 . 66 n 79 . 

.83326 

.79663. 

.82666 
12.19397 
7.31627., 
. 3.1^26.. 

,061b! 

. 1 4243_ 



,02168 

.02915 

”__.O1072_ 
.03811 
. 10776 

. 0 60 51 __ 

■ 1 H 395 
.02179 
. .-O3A70_. 
.01893 
.03056 

03.1 80 

,7b57b 

.00760 

- 00510 . 

,00650 .. 
.07275 . 

. 25760.. 

. . 3.1387L 
.00307 
.023 15 


40 C THAR 

41 OTHRL 

42 DFLfL-... 

43D5B 

44 0 TO G 

45 fi _TO__D_ 

46 LONG 

47 L AT 


PS 1 Q 

PFO 

• 147 

.117 

. 12796 ... 

. 12831 . 

.00960 „ 

8355 

- 0 . **4(* 

-9.054 

-0,90240 

8.90402 

. 16993 

8 355 

dec’ 

- 7. 301 

-8.430 

r 7,96140 

..7,96819 

.32666 

8355 

8355 

8355 

DEG 

" DEG 

5.01? 5.033 

” -.372 " -.399 

6.43606 

-. 3033 1.. 

5 .4 4035 . 

. 38339 

_ .00782 .. . __ 

PEG ’ 

1 C. 91 4 

10.606 

10. 754? 4 

10.75 4 47 

.07044 

8355 

8355 

~PCT h AX 

6 7.5 70... 

67. 303 

67.44859 

67.44B62 

.0 59 91 .. _ 

P CT MAX 7.7. 

68. 066 

67.773 

6 7 1 89 629 

67,89634 . 

,082 83 

8355 

" PUS 1 T I ON _ 

.250 

.,2 30 

.24308 __ . 

.24310 

, 0031 2 

.8 355 


POSITION. 

HETFRS 

DECREES 


,3i0. ... .306 

~7* 19 7 719. 11 7 74860(3. 0 t. 


DECREFS ._ 

48 TRK_ANG DE G P E FS 

49 HDG KAD IANS 

50 VE .... n/SEC .. 

51 V N P/S.EC 

52 ALTITUDE ** 

53 TFPPC _ _ DEGREES C. 

54 F h HKD .5 PD. _J<NQJ$ 

55 NS HhD S PD . KNOTS 

56 t* TNP 5PFE0 ANUT5 

57 WIND 01 R EC DEGREES 

5 B HIND DEGPFES 

59 wIND PIR3._. DEGREES 


-ua. 178 

35.725 

.190, 153 

3,6 JU. . 
-22*0/4 
-92.662 
1,(82 . 
9,UJb _ 

_43. 302 

10,3/9 . 

4 3 . 1 7 0 _ 
_35 9,*ij5 . 

1 79,4 35 
359,035 . 
797, 759 


.34722 

I?.p i9 72.4352? _7UJ5J4« 

10.239 -U8.20B02 UB. 20802 


3265.93300 
..OOifiL 


AIR 

speed' r' 

f/ SEC 

* 1 10 , 

-401 . . 

A1P 

SPFEO C 

p/sre 

1 OH, 

, 4i6 

AIR 

5PEEP L_ 

N/src 

1 11 . 

, 0 7 e._. 

PF L 

TA ALT 

METERS 

135, 

,311 

INRTL DI5P 

METFRS 

136, 

>30? . . 

uc_ 

RIGHT M/SEC. 

M, 

*343 

UC 

CENTER 

H/5FC 

7, 

,011 

UC 

LEFT 

M/SEC 

0, 

,272 

1 VG 

R I CtlT 

H/SEC 

1 7. 

.6 10 

1 VC 

CENTER 

H/SEC ... 

15, 

,010 

VG 

If FT 

M/SEC 

15, 

. 304 

, wc 

RIGHT _M/StC 

10 , 

. 700 

1 HC 
I «C 

CENTER 
LEFT T, 

K/5EC 

fr*/SEC 

10 , 

12 , 

.126 
.5 75 


35.23301 
_.X95.2 5335 _ 
3.52647 
- 2 b. 461 7t 

-96. 87623.. 

. 1 » 54 676 
6.55501 

20,14629 

-5.49680. 

21.51716 . 

_ ?H }.b207J . 
.7 101.62077 
_ .283.62077 
AHn. 220*1 1 _ 
Vo 2V999 7?" 
' 101.0/71*7 
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Figure A. 72. Time histories of gust velocities, gust velocity differences, and aircraft 
normal accelerations. Flight 31, Run 15. 
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Figure A. 72. (continued). 


TABLE A. 29. Average Turbulence Parameters, Integral Length Scales, and 

Correlation Coefficients of Gust Velocities, Flight 31, Run 15. 


1. Mean Airspeed (m/s): 

V L V C V R 

107.74 105.40 107.23 

2. Standard Deviation of 
Gust Velocities (m/s): 


CT uR 

ff vR 

a wR 

7.46 

2.84 

3.45 

a uC 

ff vC 

ff wC 

7.31 

2.89 

3.29 

a uL 


a wL 

7.32 

2.87 

3.35 


3. Standard Deviation of Gust 
Velocity Differences (m/s): 

a AuRL ct avRL ff AwRL 
1.45 1.24 1.49 


4. Integral Length Scale (m): 

LuR LyR t-wR 

540.0 225.8 526.1 

luRL l vRL l wRL 

526.5 225.3 494.0 

5. Correlation Coefficient of Gust 
Velocities: 

u R u l/ ct ur°U[_ 

0.90 

u R v r/°ur 0 vr 
0.32 

TiRVL/q UR gv L 

0.30 


VRVL/a VR q VL 

w R w L/a WR °w L 

0.88 

0.88 

VRWR/a VR a WR 

WR u R/a WR a u R 

0.06 

0.02 

v R w L/ a VR a W[_ 

W R U L/ a w R a U[_ 

0.01 

0.00 
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Figure A. 73. Probability density functions for gust velocities and 
gust velocity differences (normalized with the standard 
deviation). Flight 31, Run 15 (r = degree of non-Gaussian) . 













One- and two-point cross-correlations 
Figure A. 74. (continued). 







1E+02 



Figure A. 75. Comparison of normalized one- and two-point spectral density functions 
for gust velocities with theoretical models, Flight 31, Run 15. 
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One- and two-point cruss-spectra. 
Figure A. 75. (continued). 
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TABLE A. 30. 


List of All Parameters Measured and Their Range of Values 
FI i ght 31 , Run 1 5 . 



CHANNEL 

1 T l»r 

2 pm! pnt 

3 *T' rr, 

A tuft 4 not 

5 TMFT4_ 

6 PM|_ 

7 P<V 1 

s 

9 Tjr 

lOlfi 

11 TtrL_w 1 T 

12 *CCi F ST 

13 *cn t re 

14 V 6$ c liwTtS 

11*4* ettf *40 
16 ** T* CTS - 

17 T^p“r~ "Dfrr 


oc? L T 

pmd 


34 OC > ST 

P5Tr> - 

6910 ' 


»T»7 

n?G 

-5 

36 OIL 

PEG 

-5 

5 PPLPV 

DEG 

6 

M O^Ti4 

„ D*ur — 

T)?G 

PEG 

11 

40 hTHS* 

PC T MAX 

69 

41 nfMft ’ 

PCT HAX 

“69 


42 r>us 


‘ pns tTinn 


.nr 


.230 


43 

PnSITTON 

,361 



hMfrS' 

7506561.672 7436279, 

4 c; i *>n n 
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46 1 ^ 

47 ** 4T 

0 FCSF f 9 
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4.011 
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51 VN 

H/5FC 
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K H 
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1 ."50 

Q , 5 

1 
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55>I9“wmp srn 
S6 VT«0 «PTFP 

KNOTS 

KNOTS 
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Figure A. 77. Time histories of gust velocities, gust velocity differences, and aircraft's 
normal accelerations. Flight 31, Run 16. 



TABLE A. 31. Average Turbulence Parameters, Integral Length Scales, and 

Correlation Coefficients of Gust Velocities, Flight 31, Run 16. 


1. Mean Airspeed (m/s): 

VL Vc Vr 

109.41 107.07 108.82 

2. Standard Deviation of 
Gust Velocities (m/s): 


ff uR 

CT vR 

ff wR 

5.68 

3.21 

3.21 

a uC 

ff vC 

a wC 

5.59 

3.44 

3.02 

a uL 

ff vL 

a wL 

5.74 

3.29 

3.14 


3. Standard Deviation of Gust 
Velocity Differences (m/s): 

a AuRL ff AvRL a AwRL 
1.75 1.61 2.00 


4. Integral Length Scale (m): 


luR lyR IwR 

348.1 362.2 95.0 

LuRL *-vRL L w r[_ 

347.5 336.5 115.3 


5. Correlation Coefficient of Gust 


Velocities: 

TI R D L/ CT UR a U|_ 

VRVL/a VR a VL 

w R w l/ a WR cy WL 

0.85 

0.86 

0.85 

^TO/^ur^vr 

^R w r/ ct vr c wr 

WrUr/^wr^ur 

0.49 

0.05 

-0.10 

TJRVL/a UR a VL 

v R w L/ a VR CT w L 

w R u L/ a WR a U[_ 

0.49 

0.02 

-0.10 
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Frequency (Hz) Frequency (Hz) Frequency (Hz) 

a. One- and two-point common component spectra. 

Figure A. 80. Comparison of normalized one- and two-point spectral density functions 
for gust velocities with theoretical models, Flight 31, Run 16. 
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One- and two-point cruss-spectra. 
Figure A. 80. (continued). 
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TABLE A. 32. 


List of All Parameters Measured and Their Range of Values, 
FI i ght 31 , Run 16 . 


CHANNEL 

UNITS 

HIGH 

LOW 

MEAN 

RMS 

STO 

POINTS 

1 me 


*4650*419 

44530,619 

44600.51910 

44600 .526*1 

26.62063 

3993 

2 PHT nnt 

RAD/SFC 

.173 

-.170 

-.00219 

.05469 

.05465 

3993 

3 i^CL 5 cg 

G UNTTS 

1**36 

.184 

.99295 

1.02C51 

,23560 

3993 

' THf nnT 

R AD/SEC 
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-.115 

_ _ .00353. 

.02160. 
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3993 

^ ThFTA 

PAD 
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-.028 
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• 0 2 6 9 7 

3993 

9 PMT 
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. 12 * 

-.141 

.00010 

• 0 5 3 C 3 
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7 Ml ! 

PAD *" 
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3993 

8 DPI P*T 1 

DEG 
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-3.444 

2.57783 

3.36447 

2 • i 6232 

3993 

9 * 

RAD 

456. 5il 

446.301 

452.21091 
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2.11081 

3993 

iDnpi P9 f ? 
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7.40? 
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.00199_ 

•04261 
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APPENDIX B 


DERIVATION OF EQUATIONS 

This appendix contains a more complete derivation of the gust equations to 
compare with the specialized form of those used by the NASA Langley Research 
Center (LaRC) and to show the significance during certain manuevers, of terms 
which are not present in the specialized form. The wind velocity vector 

components at some position r measured from the c.g. of a rigid aircraft are 
designated u g , v g , and w g . These are measured In the coordinate system with 
the x axis pointing north, the y axis pointing east, and the vertical axis 
pointing along the local vertical (gravity vector; positive downward). The 
coordinate system Is called the true north coordinate system and Is taken as 
the Inertial system in this analysis (however, see Rhyne, 1976).* The u g , v g , 
and w g components point north, east, and vertical, respectively, and are given 
by: 

U+Ugl fV N 1 fUR 

V+Vg = + VR (B.l) 

w+Wgj [Vaz] [ w R 

The symbols u, v, and w designate the components of the aircraft velocity 
vector relative to the air mass measured In the true north coordinate system; 
Vn, Ve, and V^z are the inertial velocity vector components of the c.g. of the 
aircraft; and ur, vr, and wr are the velocity components of the position r 

*Gr1d north is true north at the platform alignment location, but as the 
platform moves east or west from Its Initial alignment point, Its north-south 
axis Is not torqued to point at true north but remains parallel to a vertical 
plane through the meridian at which it was aligned. (The north-south and 
east-west axes are torqued to be perpendicular to the local vertical at all 
times, however.) For all practical purposes, the inertial-platform axis 
system can be assumed to be aligned with true north, considering the latitudes 
of operation and the east-west distances flown in a preceding project. 
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relative to the c.g. of the aircraft due to rotation of the frame of reference 
fixed in the airplane, I.e., the body coordinate system. 

The matrix L/\i transforms the velocity components in the true north 
coordinate system to the average flight path coordinates. This transform 
matrix has the following form: bank angle (¥), track angle (¥). and eleva- 

tional angle (" 0 ). 

cos ¥ cos ¥ cos 0 sin ¥ -sin 0 

sin ¥ sin ¥ cos ¥ sin ¥ sin 0 sin ¥ sin <|> cos 0 

l AI = - cos <J> sin ¥ + cos <t> cos ¥ (B . 2) 

cos ¥ sin ¥ cos ip cos ¥ sin ¥ sin ¥ cos ¥ cos 0 

+ sin ¥ sin ¥ - sin <J> cos ¥ 

The velocity components ur, vr, and wr In Equation B.l are derived 
as follows. The velocity of a point r = JL x i + &yi + &zi< measured in the 
airplane frame of reference (I.e., body coordinates) which is rotating 
relative to the fixed frame of reference (I.e., Inertial frame taken as the 
true north coordinates In this report) is given by ^ * f, where $ Is the 
angular velocity of the airplane frame of reference relative to the inertial 
frame of reference. 3 has the components p, q, and r and ft x r has the 
components ur, vr, and wr expressed In body coordinates, I.e.: 

[ 13*1 

ft x r ■ p q r = ur^ + vrJ + wr£ (B.3) 

.Jt x 2-y &z- 

Note H z measured down is positive and JLy measured to the right is positive. 
Expanding Equation B.3 gives: 

ur] [d a z " r V 

vr = r£ x - pJt 2 (B . 4) 

wrJ = IpJly - qS-x 
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In terms of the Euler angles (ip,0,4>) of the body axis relative to the true 
north or Inertial frame of reference: 


p 


'1 

0 

-sin 0 

q 

= 

0 

cos <p 

sin <p cos 0 

r 


0 

-sin <p 

cos <p cos 0. 


hence 



6 


(B.5) 


(B.6) 


Thus, the components of the rotational velocity of the position r about the 
c.g. measured In the body coordinate system are: 


2, z (9 cos <p + sin <p cos 6) + fi.y(0 sin <p - ip cos <p cos 0) 

-[Jt z (4» - sin 0) + Jl x (0 sin <p - ip cos <p cos 0)] (B.7) 

Jt y (<p - ip sin 0) - J, x (0 cos <|> + sin <p cos 0) 

Now since these are velocity components in the body coordinate system they 
must be transformed to the average flight path coordinates. 


UR 


VR 

= 

WR. 



^ I 


[UR] 


UR' 

► < ) 
PJ 

= l AB 

V R 

[WR 


WR 


Lab 


cos 0 cos 


cos 0 sin ip 


-sin 0 


sin 4> sin 0 cos ip 
- cos <p sin ip 
sin $ sin 0 sin ip 


+ cos $ cos ip 


sin $ cos 0 


cos <p sin 6 cos ip 
+ sin $ sin ip 
cos $ sin 0 sin ip 
- sin $ cos 
cos $ cos 0 


(B.8) 


(B.9) 


where ur, vr, and wr are the components of the rotation vector expressed In 
the average flight path coordinates; ip, 0, and $ are the Euler angles of the 
body axis relative to the average flight path axis. 
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The north, east, and vertical Inertial velocity components expressed In 
the average flight path coordinates are denoted with capital letters having a 
C) are: 


V N cos Q cos + V£ cos Q sin 4> - sin 0 

V^(sin 4> sin 9 cos i> - cos 4> sin Tp) + V^(sin $ sin 9 sin 

+ cos $ cos 4?) + sin 4> cos 0 (B.10) 

V^(cos J sin 0 cos + sin $ sin ijJ) + V^Ccos $ sin 0 sin 

- sin $ cos 4>) + V /\2 cos ^ cos 0 

The values of u, v, and w which are the true airspeed velocity components 
In the average frame of reference are not measured directly In the flight 
experiments. Rather the true airspeed of the aircraft, V, is measured. 
Therefore, u, v, and w must be expressed In terms of this variable. The 
velocity components u', v', and w 1 (i.e., measured In body coordinates) are 
related to the true airspeed by the relationship: 


(B.ll) 


u' 


V 

V* 

= Lew 

0 

w 1 . 


.0. 


where 


l BW = 


[cos a cos B 
sin B 


-cos a sin B -sin a 

cos B 0 (B.12) 

[sin a cos B -sin a sin B cos a 

and a and B are the angle of attack (= tan"* w'/u') and sideslip angle (= 

sin - ! v'/V), respectively. Lgw transforms the velocity components measured in 

a frame of reference for which the x axis Is located along the relative 

velocity vector (Etkin (1972) calls this the "wind" coordinate system) to the 

body coordinate system. Thus: 



(B • 13) 


1 

u 


V cos a cos B 

V* 

= 

V sin B 

1 

IW 


V sin a cos B 


The above assumes that the pitot tube measures actual magnitude of the 
relative velocity or true airspeed and not some fractional component. 

The above values must be rotated into the average flight path frame of 
reference with the transform Lab» i.e. f 

I 


u 


u 

A 

V 

= l ab 

1 

V 

w. 


V . 


The wind velocity measured 
by: 


(B. 14) 

in the flight path coordinate system is thus given 


[“gl 


'Ul 

A 


u R 


U' 


[ VN 1 


UR 

1 

- u 

cn 

<> 

= 

V 

A 


VR 

- 

V 

= l ai 

Ve 

+ l ab 

VR 

- V 1 

lw q J 


.w. 


.WR 


w. 


lv a2 l 


WR 

- w\ 


(B.15) 


Consider the transform Lab (Equation B.9). The angles <J>, 0, and iJj are not 
measured In the flight program; therefore, L/^q must be expressed In terms of 
<j>, e, and which are measured and 4>, ?, and i|> which may be determined In 
post-flight analysis. This is achieved as follows: 


Va = LaiVj and Vj = LjrVb (B.16) 

hence 


Va = l ai l ibVb (B.17) 
and thus 

l AB ■ l AI l IB = l AI l BI (B.18) 
where the superscript T denotes the transpose. The terms of Lab f° r the 
general case are very complex; however, assuming wings level flight, i.e., <j> = 
0, which does not lose any generality for the present problem, results in: 
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COS 0 COS 0 COS (tMO 
+ sin 0 sin ¥ 


sin 4> sin 0 cos 0 
cos (4H>) - cos 0 
cos <j> sin (4>-4>) 

- cos 0 sin <{> sin 0 


cos <}) sin 0 cos 0 
cos (4H>) + sin <j> 
cos 0 sin (tM<) 

- cos 0 cos <}> sin 0 


l AB 


cos 0 sin (4>-4>) 


cos 0 sin 0 cos (4>-4>) 
- sin 0 cos 0 


cos 4> cos (4>-4>) 

+ sin <j> sin 0 
sin (4» - ip) 

sin <j> sin 0 sin 0 
cos (4’-4>) - cos 4> 
sin ¥ sin (iM>) 

+ cos ¥ sin <|> cos 0 


-sin <t> cos (4>-4>) 

+ cos <j> sin 0 
sin (4'-40 

cos 4> sin 0 sin ¥ 
cos (4H>) + sin <J> 
sin ¥ sin (4H>) 

+ cos 4 > cos 0 cos ¥. 

(B.19) 


Now assuming 4>, ip-ip, and 0-0 are small angles and neglecting high order terms. 
Equation B.19 reduces to: 


l AB 


1 -(4>-4>) cos ¥ 0-¥ 

(4>-¥) cos 0 1 -4> 

-(0-0) 4> - (4>-40 sin 0 1 


(B . 20) 


Substituting L AB from above into Equation B.15 and similar assuming small 
angles (or angle differences) with second order terms neglected in the 
expressions ur - u'» vr - v 1 , and wr - w 1 (see Equations B.7 and B.13) the 
second term on the right-hand side of Equation B.15 becomes: 


ur - u'] 


l AB 


V R 


I 

V 


wr - w' 


Z z 9 - fcy(-04> + 4> cos 0) - V - (i|>-¥) cos ¥[j, x 4> cos 0 - VB 

- l z (4> - 4) sin 0)] + (0-¥)[Jly(4> - 4> sin 0) - S. x 0 - Va] 
(4>-4i) cos 0[A Z 0 - fi. y 4> cos 0 - V] + Jl x (-04> + 4» cos 0) 

- VB - Jl 2 (4> - 4> sin 0) - 4>[Jly($ - 4> sin 0) - H x 0 - Va] 
-(0-¥)[H z 0 - Jt y 4> cos 0 - V] + C4» - (4>-4») sin ¥][£ x 4; cos 0 

- VB - M4> - ^ sin e )3 + - 4 sin 0) 

- Jl x (0 + 44 cos 0) - Va 


(B.21) 
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The derivation of the equations currently used in the data reduction 
algorithms at the NASA Langley Research Center computer laboratory treats the 
values of <j>, 0, and i> as small. Moreover, assuming the position vector r = 
jt x i + Jtyj + lies in the x- y Pl ane of the body coordinate system, i.e., 
jt z = 0 and introducing these assumptions into Equation B.21 gives upon 
neglecting higher order terms: 



1 1 „ 
UR - u 


-Jyp cos 0 - V + (iM>) V0 COS 9 - (0-0) Va 

l AB 

V R ” v ’ 

= 

- ('Mp) V cos 0 + Jl x ip cos 0 - VB + V<(>a 


.WR - wV 


(0-0) V + [<j> - (lp-^j) sin ¥] VB + ~ £ x 0 - Va. 


Recalling that we have assumed wings level flight, I.e., <j> = 0, the first 
term on the right-hand side of Equation B.15 becomes: 


A 


U 


'Vn cos ¥ cos ip + V£ cos 0 sin ^ - V^2 sin 0' 

A 

V 

= 

-Vn sin ip + V^ cos 

W: 


V^ sin 0 cos ip + V^ sin ¥ sin + \lpj cos 0. 


(B.23) 


Therefore, under the following assumptions: 

1. Wing level flight, I.e., 4> * 0. 

2. (j), 0-¥, i|>-ip, a, and & are small (<10 deg, cos () = 1, sin () = (); 
error <2%) and high order terms of the products of these angles are 
negligible (error <3%). 

3. The wind velocity probe Is measured at a given point in the x-y plane 
of the body coordinate system (I.e., % z = 0). 

4. The values of <j>, 0, and ip are small (<10 deg/sec, error <2%) and high 
order terms of the products of these values are negligible (error 
<3%). 

The wind velocity vector components expressed in the average flight path 
coordinate system Is given by adding Equation B.22 and Equation B.23: 
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(B.24) 


A 



V N cos ¥ cos ¥ + V£ cos ¥ sin ip - Vaz sin 0 - £y ip cos 0 

- V + (\p-ip) V0 cos ¥ - (0-¥)Va 

-Vfs| sin ip + V£ cos ¥ - (lp-ij)) V cos 0 + i> cos 0 - VB + V<pa 

sin ¥ cos 4> + V£ sin ¥ sin ip + Vaz cos 0 + (0-0) V 

- [<p -(ip-ip) sin ¥]VB + Jty4> - Jl x 0 _ Va 


The NASA LaRC algorithm assumes level flight given by 0 = ¥ which 
implies the angle 0 is small. Equation B.24 then becomes: 


u g 


A 


v g 

= 

A 


w g 



Vn cos \p + V£ sin ip - fly^ - V + (ip-ip) VB - V0a 
V|r cos ip - V N sin ip - (ip-\p) V + i X Tp - VB + V<j>a 
Vaz + V0 - V<pB + £ y <P - M " v “ 


(B.25) 


These equations represent the total wind velocity components, however, 
interest is generally In the fluctuations about the mean, hence the terms in 
Equation B.25 are expressed as a mean quantity plus a fluctuation quantity, 
i.e., A = A + A 
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u g 
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u g' 


^ cos ip + V£ sin ip - V ' 


V^j cos ip + Vg sin ip - V - Jlyip 





_L_ _1_ _L_ 

- V 0 a 


+ ipVB - V0a 
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v g 

+ 

A 

v g 

= 

V £ cos ip - V^| sin ip 

+ 

V£ cos ip - Vn sin ip - Vip + 

i 1 






- VB + V4>a 

A 

Wg* 


A 

Wg * 


_L_ Jl L- -i_ JU 

IVaz + V 0 - V a 


Vaz + vo + V 0 - v<}>b + fly*!* 




. - JL x 0 - (Va + Va) 


(B.26) 


Note ip* ip, ip = Tp, 0 = 0, 4> = |, ij> = B = 0 (thus, B = B and 4 > = <p): e and a 
are not necessarily zero. The right-hand most term Is the velocity 

fluctuation about the mean where the mean is given by the expression 
Immediately following the equal sign. 

The equations used In the NASA LaRC algorithm for the fluctuating gust 
velocities are given by: 

333 



(B . 27) 


Ug 


V|\| cos ip + V£ sin - Hyp - V 

A 

v 9 

= 

A# W « f ** 

V £ cos 4> - Vn sin - V^h- £ x 4> - V0 + V<J>a 

A 

w g 


A* A* 1 ** 

Vaz + V0 - V<{>& + S,y4> _ 


In the NASA LaRC algorithm, the signs of B» V A z, and w are defined opposite to 
those used in the previous derivation. Therefore, to obtain the exact form of 

A 

the NASA LaRC equations, one has to change the signs of B, V A z, and w g In 
Equation B.27. Also, it must be noted that the values of a and B in Equation 
B.12 are measured relative to the body axis of the aircraft whereas a and B 
measured in the NASA Gust Gradient Program are relative to the axis of the 
boom. To obtain the angle of attack, a, in Equation B.12 one must add the 
angle between the projection of the boom in the body x-z plane and the body 
x-axis to the measured a. Since for the full equations a must be the value 
relative to the body axis, the angle between the boom and the body x-axis was 
estimated by subtracting the average measured value of a from the average 
value of pitch angle for the total number of straight and level runs of Flight 
31 and Flight 21. This value was determined to be approximately 4.4 degrees 
for the center probe. 

There are some differences in Equation B.26 and B.27 that can be 

explained as follows. The terms ipVB and V0a in Equation B.26 are 

neglected in Equation B.27. This is consistent with the assumption that 
second-order small terms are negligible. However, based on this reasoning, 
the terms V<f>a and V<t>B should also be neglected but It Is not. The reason 

Is that in early studies, V4>a and V<J>B were found not to be small compared 

to the other terms in the equation (Rhyne, 1976) and have therefore been 

** i m ** i 

retained. Also, the expressions Va + Va and V0 + V0 in Equation B.26 are 

M Af 

simply written Va and V0 in Equation B.27. Justification for this is that 
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since a, in Va + Va for example, is a small angle even on the average, then Va 

v* 

Is negligible compare to Vo. This is reasonable in view of the fact that Va 
may be 1 to 2 orders of magnitude larger than Vo because V is typically two 
orders of magnitude larger than V whereas a is probably of the same order of 
magnitude as a. Finally, if second-order terms are strictly neglected, then 
Va should actually be Va; however, there Is no saving in computing Va since it 
is just as easy to compute Va. This is true of V9 and V& as well. 
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